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Abstract 


In this thesis, I study the sign problem appearing in the path-integral approach to strongly- 
coupled systems. This problem appears in many interesting systems, such as finite-density 
quantum chromodynamics (QCD), the repulsive Hubbard model away from the half filling, 
frustrated spin systems, and so on. In these systems, the Euclidean classical action S takes 
complex values, and the oscillatory nature of the Boltzmann weight e7’ causes the serious 
technical problems in ab initio numerical computations. In this study, I use the recently 
developing technique, called the path integral on Lefschetz thimbles, in order to understand 
the essence and difficulties of the sign problem in the path-integral formalism. 

The basic idea of the Lefschetz-thimble method is to deform the integration contours in the 
complexified configuration space using Cauchy integration theorem. Under this deformation, 
the original oscillatory integral is decomposed into the sum of nicely convergent integrals, and 
the sign problem becomes moderate. Each integration cycle in this decomposition is called a 
Lefschetz thimble. In this thesis, I first give a concise review on this method, and demonstrate 
how it works by computing real-time Feynman kernels for free theories. 

Next, I discuss the spontaneous breaking of chiral symmetry in 0-dimensional field theories 
without the sign problem. Chiral symmetry breaking is one of the characteristic features of low- 
energy QCD spectrum, and it is important to revisit this phenomenon using Lefschetz thimbles. 
I developed an efficient way to compute Lefschetz thimbles when a continuous symmetry is 
slightly broken due to the small explicit breaking term. This study also elucidates that the 
Lefschetz-thimble decomposition is suitable to analyze Lee-Yang zeros in the path-integral 
formulation, and the phase structure of the model is discussed for complex coupling constants. 

If one tries to analyze finite-density QCD in the same way, the sign problem appears even 
at the level of the mean-field approximation. I show that the Lefschetz-thimble decomposition 
respects the anti-linearly extended charge conjugation so that the effective potential becomes 
manifestly real, and the saddle-point analysis based on this theorem solves the sign problem 
appearing in the mean-field approximation. The Polyakov-loop effective model of heavy-dense 
QCD is studied with this method as a lucid demonstration. Since all the quarks are heavy in 
this model, fermion dynamics is simplified and the mean-field treatment becomes accurate. 

The sign problem in finite-density QCD with light flavors is known to become too severe 
as the quark chemical potential goes beyond the threshold of pion mass. In order to see what 
happens there, I consider the one-site Fermi Hubbard model, which is exactly solvable but 
has the severe sign problem in the path-integral expression. In this case, the above mean-field 
approximation is not applicable, and I elucidate its reason by studying topological structures 
of Lefschetz thimbles. I also show that interference of complex phases among multiple classical 
solutions play a pivotal role to understand the sign problem and the correct phase structures 
using path integrals. I discuss this interference is also important in the finite-density QCD with 
light flavors if the baryon chemical potential exceeds the pion mass. 

In order to give a feedback of this finding to other approaches to the sign problem, I re- 


ill 


late this newly developing technique, the Lefschetz-thimble path integral, with a conventional 
approach to the sign problem, the complex Langevin method. It is known that the complex 
Langevin method sometimes does not give correct answers, and I analytically show that the 
original complex Langevin method cannot give correct results if interference of complex phases 
among classical solutions becomes important. This gives a simple criterion for incorrectness of 
the complex Langevin method without doing its numerical simulation, and suggests that the 
complex Langevin method is not directly applicable to the finite-density QCD. To resolve this 
problem, I propose to modify the complex Langevin method with some technical working hy- 
pothesis, and it is numerically tested for the one-site Hubbard model. The original formulation 
of the complex Langevin method gives the wrong answer for this model. The modified complex 
Langevin method improves the result, but there exists systematic discrepancy, which requires 
further systematic studies on its properties. 
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Chapter 1 


Introduction 


Path integral is a useful tool to study strongly-correlated systems, however it sometimes suffer 
from the sign problem in its numerical computations. In this thesis, I study the sign problem 
by applying the Picard—Lefschetz theory to unveil its essence and difficulties. First, I would 
like to give a general introduction to study nonperturbative aspects of quantum statistical 
physics using path integrals in Sec. 1.1. In Sec. 1.2, we explain the sign problem appearing 
in the finite-density Quantum ChromoDynamics (QCD), and its relevance to neutron stars is 
mentioned. The basic idea of the Picard—Lefschetz theory is introduced in order to analyze the 
sign problem in Sec. 1.3 via the simple example of the Airy integral. We show the outline of 
this thesis in Sec. 1.4. 


1.1 Path integral quantization 


Path integral is invented as one of the ways to quantize the classical system [1, 2]. The biggest 
difference of quantum mechanics from classical mechanics is the existence of interference of 
transition amplitudes. In order to deal with it, the Feynman path integral tells us that the 
time development of quantum states can be interpreted as a summation of classical amplitudes 
exp(iS{a]/h) over all possible spacetime classical paths x(t). That is, we can compute the time 
development kernel K (x;, ti; £f, tp) of a quantum particle in the potential V(x) by performing 
the infinite-dimensional formal integration, 


Rasa [Pro ; [ dt k (2) -v@) , (1.1.1) 


with the boundary condition z(t;) = x; and x(t) = xp. This reformulation of the Schrödinger 
equation is conceptually interesting, and it becomes easy to observe connection between classical 
and quantum mechanics. 

By replacing the time t with the imaginary time —i7, one can study the thermal equilibrium 
property instead of the quantum dynamical processes [3]. This relates the d-dimensional quan- 
tum systems with (d+ 1)-dimensional classical statistical systems, and we can study statistical 
physics using computational techniques of quantum physics and vice versa. In the path-integral 
representation, the partition function Z(T) is given by 


-fas k (2) vie] 


1 


Z(T) = [cow (1.1.2) 


, 
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with the periodic boundary condition x(0) = x(@), and 6 = 1/T is the inverse temperature. 

From the viewpoint of numerical computation, however, there is an outstanding difference 
between the real-time formalism (1.1.1) and the imaginary-time one (1.1.2). The integrand of 
(1.1.2) is positive definite, which means that we can regard it as a probability density. Instead 
of solving the quantum system exactly, we can create the ensemble of spacetime paths x(T) 
according to the Boltzmann weight and solve the system in a stochastic way. This is called the 
Monte Carlo simulation of the path integral, and its computational cost scales algebraically with 
the system size in this method even when the system consists of many particles. Many thermo- 
dynamic quantities can be computed for various systems using the path-integral Monte Carlo 
method, such as finite-temperature quantum chromodynamics (QCD) in hadron physics [4], 
and liquid helium [5], ultracold atomic gases [6], Bose-Fermi mixtures [7] in condensed matter 
physics. 


1.2 QCD at finite chemical potential, and sign problem 


Quantum ChromoDynamics QCD is a part of the Standard Mode in particle physics, and 
is a non-Abelian gauge theory with the color gauge group SU(3). The gauge bosons are called 
gluons that belong to the adjoint color representation, while the spin-5 particles called quarks 
are in the fundamental color representation. QCD describes strong interaction among quarks 
and gluons, which are confined inside of hadrons. Therefore, this theory plays a pivotal role to 
understand properties of nuclei, compact stars, and also early Universe. 

The QCD Lagrangian has quite a simple form, 


Loco = Lym(A) + 9[(Da+m) — pay] a, (1.2.1) 
1 
Lym = qe Ew), (1.2.2) 


where Fi, = 0,A,—O,A,+i[A,, Av] is the field strength of the SU(3) gauge field, q is the quark 
Dirac field, and I) = y,(0, +iA,) is the Dirac operator (See, e.g., Refs. [8-11] as a textbook). 
Here, fg is the quark chemical potential. In order to study the equation of state of QCD, one 
have to compute the partition function 


| DADP (- fate Loon) 


l DA Det |(P4 +m) — u74] exp (- / dtz Lx) (1.2.3) 


Zocp(T, Ha) 


In the second line, quark and antiquarks are integrated out, and we have to take the statistical 
average in terms of the gauge field A to find the thermodynamic properties. 

Let us first describe the properties of QCD in the vacuum and its spectrum. One of the 
remarkable property of QCD as a four-dimensional gauge theory is its asymptotic freedom |12— 
14]: The gauge coupling constant g depends on the renormalization scale Q, and the one-loop 
renormalization group shows that 


g(Q)? 1 


as(Q) = An ~ Amy a(R? Abon) 


(1.2.4) 
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with o = (11N. — 2N+)/487°. In QCD, the number of color is Ne = 3, and the number of 
fermion flavor is Ny = 6, and thus 69 > 0. Agen is the typical energy scale of strong interactions, 
and its appearance in quantization is called the dimensional transmutation. The world average 
value is Agcp = 213(8) MeV when a, is measured at the Z boson mass using the modified 
minimal subtraction (MS) renormalization scheme based on the next-to-next-to-leading-order 
QCD predictions |15, 16]. This means that the QCD coupling constant becomes logarithmically 
weaker as the energy scale Q? —> oo, and this property is called asymptotic freedom. Thanks to 
the asymptotic freedom of QCD, we can study the deep inelastic scattering experiments [17, 18] 
based on the first-principle computations, and the parton picture inside nucleons can be justified 
in high-energy nuclear collisions. 

At low energies Q S Aqcp, physics of strong interaction is much more difficult. Free quarks 
and gluons have never been observed even though QCD Lagrangian contains them. Because of 
the nature of strong interaction, they cannot live longer than the timescale ~ 1/Aqgcp. All the 
quarks and gluons are confined inside hadrons, and this is called the quark confinement [19]. 
Thanks to the lattice formulation of gauge theories [19] and its Monte Carlo numerical sim- 
ulation [20-22], we have a solid theoretical foundation to study the nonperturbative aspects 
of low-energy QCD. Analytical study of nonperturbative gauge theories is still a big dream 
in theoretical and mathematical physics. Confinement of color degrees of freedom has a close 
relation with topological aspects of gauge theories [23-28]. Indeed, development of supersym- 
metric gauge theories elucidates that confining phases have rich topological structures [29, 30]. 
In order to extend this knowledge for non-supersymmetric cases, we must look more carefully 
at the microscopic dynamics of gauge fields, and there is much effort in that direction [31-44]. 

Let us briefly mention the phenomenological aspects of low-energy QCD spectrum. In the 
quark model, conventional hadrons are classified into mesons and baryons: 


meson (qg) --- m, K, 1, ete. 


baryon (qqq) --- WN, ÈX, E, ete. 25) 


hadron l 
Mesons are bound states of a quark q and an anti-quark @, and thus they have no baryon charge 
Qp = 0. The lightest meson is the pion (7), whose mass m, is about 140MeV. Baryons are 
bound states of three quarks q, and thus they are fermions with the baryon number Qg = 1. 
The lightest baryon is nucleon N, e.g. proton p or neutron n, whose mass my is about 940MeV. 
It is important to notice that there is a huge difference between m,/2 and my/3. In the QCD 
Lagrangian, the light flavor of quarks, u, d, and s, have approximate chiral symmetry under 
q = e} q because their masses are smaller than the strong scale Agcp. At 2 GeV, the current 
quark masses of u, d, and s in the MS scheme are m, ~ 2.15(15) MeV, ma ~ 4.70(20) MeV, and 
Ms ~ 93.5(2.5) MeV, respectively [16]. This chiral symmetry SU (nf) x SU (np) is spontaneously 
broken to SU (np) in the QCD vacuum (ny = 2 or 3), 


(Gq) = (—230 MeV)’, (1.2.6) 


and pseudo-scalar mesons, such as pions, can be regarded as Nambu-Goldstone bosons [45, 46]. 
Because of the small explicit breaking of the chiral symmetry, those Nambu—Goldstone bosons 
become massive. The pion mass, My X \/Agcp™Mux,a, is parametrically much smaller than Agcp, 
which sets the natural mass scale of other hadrons. 


QCD matters The possibility of existence of quark matters inside neutron stars has been 
discussed since the discovery of the asymptotic freedom [47]. Since the coupling constant be- 
comes smaller at high energies or densities, quarks and gluons are expected to be liberated from 
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confinement. Based on such an idea, a prototype of the QCD phase diagram was proposed in 
Ref. [48]. The transition temperature there was given by the limiting temperature of the sta- 
tistical model of hadrons, which comes from the exponential growth of hadron resonances [49]. 
For a modern version of the QCD phase diagram, see Ref. [50]. 

Finite-temperature QCD at ju, = 0 can be studied using Monte Carlo simulations of lattice 
QCD [4]. Around the physical quark masses with (2+ 1)-flavor, lattice QCD suggests that there 
is no phase transition between hadron and quark-gluon plasma phases [51-53]. Thermodynamic 
functions show rapid crossover around the pseudo-critical temperature Tpe ~ 160 MeV, and it 
is a big surprise to have a continuity of two seemingly-disparate systems. Only at sufficiently 
high temperatures T > Tpe, does the system become weakly coupled plasma. Near the pseudo- 
critical temperature, the system is most likely strongly coupled, and the non-perturbative 
treatment is inevitable. This system is also studied experimentally by colliding ultra-relativistic 
nuclei, and its phenomenology and model study are developing [54]. 

However, the situation becomes totally different for the finite-density QCD. We have no 
lattice QCD simulation of finite-density nuclear matters, and this is because of the big obstacle 
called the sign problem. Before its explanation, let us mention physical importance of finite- 
density QCD. One way to realize a high-density QCD matter is to mash up nuclei by external 
forces, like gravity, and the neutron star is a strongly promising candidate of a laboratory for 
its study. We must have rather stiff neutron star cores in order to explain recently observed 
two-solar-mass neutron stars [55]. The origin of large magnetic fields of magnetars is also still 
veiled. Forthcoming observation of gravitational waves from neutron star mergers is expected to 
provide many clues about high-density nuclear matters. Thus, its ab initio calculation becomes 
more and more important to understand dense nuclear matters. 


Sign problem and Silver Blaze phenomenon At finite uq, the quark determinant becomes 
oscillatory complex functions of the gauge field A. In order to see this fact, let us study the 
properties of the Dirac spectrum. Let us set D(A, 4) = D4 — Yaflg +m, and consider 


[Da a Yaltg| Wie Aig = iA,,(A, iq) Pra Bie (1.2.7) 


Here, iàn (A, Hq) is the n-th eigenvalue and Yn;A u is its eigenfunction. Let us multiply ys from 
the left on the both sides of (1.2.7), then we obtain that 


[Da A Yaltq] VV n:A,tq = —iAn(A, Hq) Y5 Vn;A, ua (1.2.8) 


Therefore, det(D(A, ua)) = [],,(An(A, Hq)? + m°). The quark determinant is positive definite 
if uq = 0 because y,D,(A) is anti-Hermitian and A,,(A,0) € R. On the other hand, this does 
not hold for u4 # 0, and An (A, uq) is complex in general. This means that the integrand of 
the QCD partition function does not accept the probabilistic interpretation, and importance 
sampling breaks down [56]. This is called the sign problem, and the same problem appears in 
many other condensed matter systems [57], such as non-Fermi liquids, the repulsive Hubbard 
model away from half filling, frustrated spin systems, etc. 

One way to perform the Monte Carlo simulation is to consider the phase-quenched QCD, 
defined by 


Ziacp\(T; Ha) = J PA [Det(D(A.14))| exp (- fst] ! (1.2.9) 


Since there is no oscillatory nature, we can generate the ensemble of the gauge field for this path 
integral and perform the Monte Carlo simulation. If the phase-quenched QCD is sufficiently 
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“similar” to the finite-density QCD, we can use this ensemble of the gauge fields to study 
finite-density nuclear matters and we have no problems in a practical sense. 

This is not the case, however. In order to explain this lucidly, we consider the two-flavor 
QCD at T = 0. In this case, the phase-quenched QCD is the same with QCD with isospin 
chemical potential because |Det D(A, u,)| is realized by assigning the positive chemical potential 
+44 to the up quark and the negative one —j1, to the down quark. Since J), is anti-Hermitian, 
we find that 


det(D(A, f4g)) = det |—D4 — Yalta +m] 
= [|] Qn(4, =H)? +m’) 


n 


det( D(A, —u4)). (1.2.10) 


Therefore, |det(D(A, u4))|? = det(D(A, ua))det( D(A, —u4)) as we mentioned. At T = 0, the 
thermodynamic partition function picks up the minimum of the free energy. The lightest hadron 
is pion (m, ~ 140 MeV), and it belongs to the three dimensional representation of the isospin. 
Therefore, Ziacp)(0, 44) has a singularity at fg = m,/2, and the pion condensation appears, 
(rt) #0, in the phase-quenched QCD for ug > m,/2. 


For QCD at finite baryon densities, nothing must happen at uq = m,/2 because pion has 
no baryon charges. In the hadron spectrum, the lightest baryon is a nucleon (my ~ 940 MeV). 
Therefore, Zgcp(0, uq) has the first-order phase transition at fg = (my — B)/3, which is known 
as the nuclear liquid-gas phase transition. Here, B is the nuclear binding energy ~ 16MeV. For 
[lg < (my — B)/3, the QCD partition function must pick up the QCD vacuum state at T = 0, 
but this is highly nontrivial from the consideration on the phase-quenched QCD. 


In order to perform the ab initio computation of the finite-density QCD, we have to under- 
stand why QCD at low temperatures experiences nothing for u4 < (my — B)/3. Only after 
comprehension of this “trivial physics” from the viewpoint of path integral, can we study the 
properties of finite-density nuclear matters by using lattice QCD. We must give much effort on 
nothing by saying it the “curious incident”, and, in Ref. [58], this fact is called with the funny, 
fashionable name “Silver Blaze problem” after the famous story [59] of Sherlock Holmes. There, 
“the curious incident of the dog in the night-time” is an essential clue for Sherlock Holmes to 
unveil the mystery. It was quite natural that the dog did nothing at ordinary times, however 
Holmes said “J had grasped the significance of the silence of the dog, for one true inference 
invariably suggests others...” . 


In our case, the phenomenological evidence clearly tells us that the partition function has no 
singularities at low temperatures for ug < (my —B)/3, but this becomes a curious incident from 
the viewpoint of the quark picture. The hadron spectrum contains the quasi-massless pions. 
Since the Dirac spectrum knows it, its real part starts from m,/2. Naively thinking, quarks can 
be excited without any cost for u4 > m,/2, and the thermodynamic state must experience the 
phase transition from the QCD vacuum. We must have much imagination and guess the reason 
why the quark chemical potential does nothing for m,/2 S H4 S (my — B)/3. Interestingly, we 
can argue that all the phase transitions are hidden by the Silver Blaze problem: The orbifold 
equivalence in the large-N, limit suggests that all the critical points in the QCD phase diagram 
must lie inside the pion condensation phase of the phase quenched QCD [60, 61]. 
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1.3 Lefschetz-thimble approach to the path integral 


Now I can concretely state the purpose of this thesis. We would like to understand how the 
partition function becomes constant on the chemical potential at T = 0 in the path-integral 
formulation. After unveiling its origin and analytical structure, I would like to propose a new 
way to solve the sign problem which can explain the Silver Blaze phenomenon. 

For this purpose, I apply the Picard—Lefschetz theory to the path integral. The Picard- 
Lefschetz theory is a basic technique in the complex analysis with many variables. This theory 
is very useful in analyzing the asymptotic behavior of multiple oscillatory integrals [62, 63], 
because it gives a suitable extension of the steepest descent method for the one-dimensional 
integration. In the context of physics, its usefulness was found out in the study of analytic 
properties of the Chern-Simons and Liouville theory [64-66]. 

In the case of one-dimensional integrals, f do exp (—S(¢)/h), the steepest descent path can 
be determined by solving the stationary phase condition, Im [S(¢)| = constant, in the complex 
plane o € C. The stationary phase condition is also important even for multiple oscillatory 
integrals, and it looks quite suitable for solving the sign problem. In Refs. [67—70], the Monte 
Carlo integration on a steepest descent cycle was proposed and numerically tested for scalar 
field theories. 


Example: Airy integral Its formulation will be reviewed in the following chapter, and let 
us quickly remind the steepest descent method in the case of Airy integral. This is the simplest 
model to understand the philosophy of the Lefschetz-thimble method. When x € R, the Airy 
function is given by 


Ait) = f £ ep (-S(6,2)), (1.3.1) 
S(¢,2) = —i ($+). (1.3.2) 


If one tries to evaluate this integration literally, this is a hard numerical task because of its 
oscillatory nature. Such a hard task is not necessary thanks to Cauchy’s integration theorem. 
We can continuously change the integration contour in the complex plane C without changing 
the result of integrals. 

We show the usefulness of changing the integration contour in Fig. 1.1 for the Airy integral 
with x = 1. If we perform the integration on the real axis, we must perform the oscillatory 
integral, which integrand is shown with red dashed line in Fig. 1.1 (b). In order to apply the 
steepest descent method, we complexify the integration variable, o € C and find the complex 
saddle points: 


045(¢,1) = -i(¢? +1) =0, = ¢= ti. (1.3.3) 


Let us change the integration contour to the steepest descent path J as shown in Fig. 1.1 (a), 
which passes the saddle point ¢ = +i. In the case of one-variable integrals, the path of steepest 
descent can be found by solving 


Im [S(¢, 1)] = Im [S(i, 1)]. (1.3.4) 


The integrand shows the nice convergence after the deformation of the integration contour, and 
it is shown with the blue solid line in Fig. 1.1 (b). Furthermore, the condition (1.3.4) tells us 
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Figure 1.1: Effectiveness of the change of integration contours in the Airy integral (x = 1). 
(a) Gray region represents the nicely convergent direction Re(S) — oo, and the blue curve is 
the steepest descent path J that connects two distinct convergent directions. (b) Dashed red 
and solid blue curves represents behaviors of the integrand Re(e~*) on R and J, respectively. 


Figure 1.2: Steepest descent and ascent paths, J+iyz (solid black lines) and K4iys (dashed 


orange lines), in the complex ¢-plane. Blue blobs show complex saddle points, ¢@ = e*'*/?,/z. 


that the oscillatory complex phase does not appear from the exponential part in the integrand 
of (1.3.1). As we already mentioned, this property is useful to study the sign problem, and 
Monte Carlo integrals on steepest descent paths were proposed in Refs. [67-73]. 


Let us consider a generic case when x takes complex values and perform analytic continua- 
tion of Ai(x). Once the integration contour in (1.3.1) is deformed to the steepest descent path 
J for x > 0, the integral expression is convergent for generic x. Let us rewrite that path J 
in Fig. 1.1 using the steepest descent paths at x. We show steepest descent paths at x = e°" 
and x = e™ in Fig. 1.2. If the phase of x is sufficiently small, or, more precisely speaking, if 
|arg(x)| < 27/3, only the steepest descent path 7,;/ contributes as we can see in Fig. 1.2 (a). 
That is, for jarg(x)| < 27/3, 


Ai(x) = f hi eke) (1.3.5) 


ive 


However, when the phase of x is close to 7, the topological structure of steepest descent paths 
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changes drastically, as is shown in Fig. Fig. 1.2 (b). For jarg(”) — m| < 7/3, 


Ai(z) = f GP o-s) 4 f d? 802), (1.3.6) 
T 2T T 2T 


tiva -ivs 


The jump between these two expression occurs at arg(x) = 27/3, which is called the Stokes 
jump. At arg(x) = 27/3, the complex phases of the classical actions at two saddle points 
become identical, 

Im [S(+i/z, x)| = Im [S(-iVz, 2)] , (1.3.7) 


since S(+iyz, x) = +22°/?. Therefore, the stationary phase condition Im(S(¢)) = constant 
connects two saddle points ¢ = +i,/x. This is called the Stokes ray, at which the topological 
structure of 74; ,/¢ changes. In order to apply the steepest descent method, tracing this change 
of the topological structure of J> is essential for o € {+i\/x}. Using the steepest ascent paths 
Ko shown with orange dashed lines in Fig. 1.2, this topological information can be encoded as 
intersection numbers. The integral on Js must be taken into account if the intersection number 
(Ko, R) between Ko and the original contour R is nonzero, and 


Aiiz)= > CoR) | pa, (1.3.8) 
o€{ti/zr} Jo 


This single equation summarizes (1.3.5) and (1.3.6). Its extension to multiple oscillatory inte- 
grals will be considered in Sec. 2.1, and the steepest descent contours J, are called Lefschetz 
thimbles. 

Based on the steepest descent method with (1.3.5) and (1.3.6), let us study the asymptotic 
behavior of Ai(x) with |z| => oo. According to (1.3.5), for jarg(a)| < 27/3, 


exp (—22°/?) 


(1.3.9) 
Since the second term in (1.3.6) is exponentially smaller than the first term, this expression 
turns out to be valid not only for |arg(x)| < 27/3 but also for |arg(xz)| < m. On the other hand, 
for jarg(x) — r| < 7/3, we obtain 


cos (2(—2)*/? — £) 


Ai(a) ~ dan 


(1.3.10) 


This expression is also valid for |arg(x) — m| < 27/3. Remarkably, there are two different 
asymptotic expressions of Ai(x) for 7/3 < arg(x) < a, although Ai(z) is a single-valued function 
on C. This is one realization of the Stokes phenomenon. The boundary of applicability of each 
asymptotic approximation is called the anti-Stokes ray, at which the dominance of different 
Lefschetz thimbles is interchanged. 


Stokes jumps, Borel summation For completeness of our analysis, let us briefly review 


the above analysis from the viewpoint of the differential equation. Let us recall that the Airy 
function is a solution of the Airy equation: 


(a - r) Zj =ð. (1.3.11) 
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One can easily see that (1.3.1) satisfies the Airy equation (1.3.11) by using the integration by 
parts. By regarding Ai(x) as a partition function Z(x), the Airy equation is nothing but the 
Dyson—Schwinger equation, 


dod. /¢ . f de q 
o= [ Seewi (5 +20) =i [Eeto (Er). (1.3.12) 


Since ¢° in the last expression can be replaced by —d?/dz?, we can obtain (1.3.11). This 
is a second-order ordinary differential equation and has two linearly independent solutions, 
commonly called Ai(x) and Bi(x). This corresponds to the fact that this model has two complex 
saddle points and thus two steepest descent contours. One can also easily see that using a new 
variable ¢ = 2r3/ 2 the Airy equation has a regular singularity at ¢ = 0 and an irregular 
singularity of rank 1 at = oo. Because of these singularities, the perturbative expansion of 
these partition functions Ai(x) and Bi(x) at infinities can only have asymptotic expansions, i.e., 
the formal power series in terms of 1/¢ are divergent. This fact leads an exotic but interesting 
fact, which is called the Stokes phenomenon. 

Let us observe an example of the Stokes phenomenon. In the limit x — oo with arg(x) < r, 
the Airy function has a formal series expansion 


exp (—223/*) = „T (k+2) 0 (k+ 2) 1 
SOn Se Era Game P 


This series is asymptotic because the coefficients show factorial growth. In order to make sense 
out of this divergent series, we perform the Borel transformation; 


seS e e tae K+ a) (1.3.14) 


co TOT a) AAAI 


Here, the series coefficient is divided by k! compared with that of the corresponding divergent 
series, and we can easily check that this series has the convergent radius |ņ| < 2. By regarding 
B(n) as an analytic function by using the analytic continuation, it has a branch point at n = —2. 
In this easy case, one can perform this summation explicitly to find that 


15 n 
B = oF, | ~,-=;1;-=}. 1.3.15 
(n) 2 (G3 , z) ( ) 


In order to relate B(ņ) with the Airy function, we perform the Laplace transformation and get 


Ai(e) = 8" [ane B(n) (1.3.16) 
25/631/6 71/2 J n n), oe 


with Ç = $2°/? and arg(¢) < 7/2. This is a well-known expression of Ai(x) using the confluent 
hypergeometric function. In this way, we can reconstruct the original analytic function using 
its asymptotic series expression, and this process is called the Borel summation. 

Using the Borel summation (1.3.16), let us observe how Stokes jump plays into the game. 
Let us add a phase C +> ¢e'®, then we also rotate the phase of the integration contour of 7 so 
that 7 € e~9[0, 00). Therefore, for =r < 0 < 7, 

i0)5/6,—Ce? —pe*oo 
Ai(xel?9/3) = Sauna f dn eS" B(n). (1.3.17) 
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As 0 — 7, the integration contour comes close to the negative real axis, while B(¢) has a branch 
cut (—oo,—-2]. For 0 > a, the additional contribution comes from the contour integration 
around the branch cut, and this phase clearly corresponds to arg(x) = ar Therefore, the 
appearance of this additional contribution is the Stokes jumping, and 0 = ~ is the Stokes ray. 
More interestingly, asymptotic expansions at two complex saddle points have some relations; 
for example, large order perturbations at a saddle point determines the small order behaviors of 
perturbations at another saddle point. This is an example of resurgence, proposed by Ecalle [74]. 
For a readable review of resurgence, see Ref. [75] for example. 


1.4 Outline of this thesis 


In Chapter 2, we briefly review the Picard—Lefschetz theory and study its basic properties in 
the context of oscillatory integrations. Picard—Lefschetz theory is established only for finite- 
dimensional oscillatory integrals, and thus we confirm that the Picard—Lefschetz theory works 
well also for free quantum mechanical systems. 

In order to understand how this technique describes the symmetry breaking phenomena, we 
discuss, in Chapter 3, the Z and O(4) chiral symmetry breaking in simple models without the 
sign problem. By relating the Lefschetz-thimble analysis with the Lee-Yang zeros of partition 
functions, we can find how the phase transition is described in the Lefschetz-thimble path 
integral for these models. 

In order to understand the origin and structures of the sign problem, we consider two 
examples with the sign problem in Chapter 4. The first one is the Polyakov-loop effective model 
of the finite-density quantum chromodynamics (QCD) with heavy quarks. In order to study it, 
we developed the mean-field approximation which is applicable even if the sign problem exists, 
and we can judge whether this works nicely by computing topologies of Lefschetz thimbles. This 
gives the mathematical foundation to study the Polyakov-loop models at finite quark densities. 
The second one is the one-site Fermi Hubbard model, which can be exactly solved and has 
the severe sign problem in the path-integral expression. In this case, the above mean-field 
approximation is not applicable, and we elucidate its reason by studying topological structures 
of Lefschetz thimbles. We show that interference of complex phases among complex classical 
solutions play a pivotal role to understand the sign problem and the correct phase structures 
using path integrals. We discuss this interference is also important in the finite-density QCD 
with light flavors if uq exceeds the half of the pion mass. 

We relate the Lefschetz-thimble path integral with the complex Langevin method, in Chap- 
ter 5. Since the complex Langevin method has longer history for the study of the sign problem, 
this establishment enables us to translate the new findings in the Lefschetz-thimble approach 
to other methods, and vice versa. We show that the original complex Langevin method cannot 
give correct results if interference of complex phases among classical solutions is important. 
This gives the simplest criterion for incorrectness of the complex Langevin method without 
doing its numerical simulation. This also means that the complex Langevin method is not 
directly applicable to the baryon Silver Blaze phenomenon in finite-density QCD. To resolve 
this problem, we propose to modify the complex Langevin method by introducing a technical 
working assumption, and it is numerically tested with the one-site Hubbard model. The original 
complex Langevin method does not give the correct answer, while the modified one shows rapid 
transitions of number densities. Although the result is not yet perfect, this study shows the 
importance of relating the Lefschetz-thimble method to other approaches to the sign problem. 


Chapter 2 


Path integral on Lefschetz thimbles 


In this chapter, we describe the application of the Picard—Lefschetz theory to the oscillatory 
multi-dimensional integrals, and its formal application to the functional integration. We study 
its generic properties, and compute it for solvable examples in quantum mechanics. 


2.1 Picard—Lefschetz theory for oscillatory integrations 


In this section, we first explain how to apply the Picard—Lefschetz theory to oscillatory in- 
tegrals [62, 64, 65]. This gives a mathematically rigorous extension of the steepest descent 
analysis given in Sec. 1.3. In this formulation, the complex analogue of the Morse homology 
theory turns out to be useful, which basics is briefly summarized in Appendix A. The procedure 
for its formal application to the path integral is summarized in Sec. 2.2. 

We consider the following integral, 


Za = f dé exp (—S(0)/ħ), (2.1.1) 
yX 


where fA is a real parameter. Here, the integration region Y is taken as a real affine variety 
with a natural volume form d0. One can think of Y = R” and dð = d"zx as the simplest 
example, but one can also think of classical real Lie groups such as Y = O(n), U (n), etc. with 
the Haar measure. Let S : Y — C be a complex-valued polynomial on VY, and Re(S) > co as 
|0| — co if Y is non-compact. In the language of physics, V is the configuration space, d is the 
path-integral measure, and S' is the (Euclidean) classical action of the system to be quantized. 
Since S(@) is complex in general, this is an oscillatory integral, which often impedes the study 
of analytic properties of the partition function Z}. Picard—Lefschetz theory provides a beautiful 
framework, which converts such oscillatory integration into a sum of nicely convergent integrals. 
This A in quantum physics appears with different symbols in different contexts; e.g., it can be 
replaced by kgT in statistical mechanics, 1/N in N x N matrix models, etc. 

We assume the existence of complexification V of Y. This means that there is a complex 
manifold X which contains JY as its subspace, and that Y is fixed under the complex conjugation 
-of X; i.e. VOX and y = y for any y € Y. For Y = R”, O(n), U(n), their complexifications 
are ¥ = C”, O(n,C), GL(n, C), respectively. By assumption, ¥ has a Kahler metric 


1 : — — . 
ds? = 598 (az @ dz + dz & dz’) f (2.1.2) 
and we denote its Kähler form as w = iggd7 A dzi. Here, (2t,...,2”) is a local holomorphic 


coordinate on ¥. 
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One of the primary purposes of this section is to construct all the possible n-dimensional 
integration cycles J in XY, on which Im(S/h) is constant so that the integrand is not oscillatory. 
After that, we represent (2.1.1) using those integrations. Let us denote J = —S/h and regard 
it as a holomorphic function on Y (A can now be complex in general). The Morse function h 
is defined by its real part: 

T+T 


h=ReD=——. (2.1.3) 


Since the integrand is given by e”~'!™(6/), this function h represents the importance of field 
configurations in the complexified space. 

Let p E€ X be a critical point of h; i.e., ih), = ihl, = 0 using a local holomorphic 
coordinate (z!,..., z”) around p. Equivalently, p € ¥ is a critical point if and only if 


ƏT], = 0, (2.1.4) 


because 0,Z = 0 follows automatically from the Cauchy-Riemann condition. This is nothing 
but the classical equation of motion. It is important to study the property of the Hessian 
matrix around the critical point. Using a complex version of the Morse lemma, we obtain that 
there exists a local coordinate around a non-degenerate critical point p such that 


T(z) = (0) + (2)? +--+ + (2). (2.1.5) 


By taking its real part, we find that the Morse function behaves around the non-degenerate 
critical point as 


h(z) = h(0) + @')? +++ a") — GP me — GT)’, (2.1.6) 


with zf = xf + iy’. Therefore, the Morse index of non-degenerate critical points of h is always 
n=dimcr. 

In order to find the integration cycles 7, it turns out to be important to compute a relative 
homology of ¥ using the Morse function h. In order to explain it, we introduce the Morse’s 
gradient flow equation by 


dz’ gOh zə da jôh OT 


==) — J — ==) J 7 = — K aris 2.1. 
dt 4 Oz I zi dt 1 Ja 1 Jg a) 


Along this flow equation, the Morse function h decreases monotonically: 


dh ie pte 

a —2g" ð;hðzh (=: —|Vh]?) < 0. (2.1.8) 
It is clear that the equality holds only at a critical point of h on X. It is of great importance 
for our purpose to notice that Im(S/h) is constant on each flow line given by the gradient flow. 
Let us denote H := Im Z = (T — T)/2i, and then 


dH dz? dzi_ g7 _ E 
de de ror OF = — gy OAE (2.1.9) 


This phenomenon has a clear interpretation based on the classical mechanics [64], and it be- 
comes useful when we extend this formalism to the case with continuous symmetries. Such an 
extension will be discussed in Sec. 3.2, and for a moment we continue to assume that all the 
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critical points are non-degenerate. In order to explain the classical-mechanical interpretation, 
we pay attention to the symplectic structure w of 4, which defines the Poisson bracket {-,-}p 
by 

{f, g}p = —2ig™ (O. f079 — Orf Ong) - (2.1.10) 


In the language of classical mechanics, this says that x’ and yê describe coordinates and canon- 
ical momenta, respectively. Now, the conservation of H becomes crystal-clear because the 
gradient flow equation is nothing but the Hamilton equation of motion: 
dz’ ; dz — 
={Hztp, —={H,z tp. 2.1.11 
di { je dt {H, zi} p ( ) 


The conservation of H is very important to define 7, because we would like to find integration 
cycles satisfying that property. 
Let us denote the set of critical points by X. It labels critical points pz, i.e., 


OZ |p, = 0, (2.1.12) 

for any ø € X. In order for J to be an integration cycle of the integration measure ei) q"z, 

it must be an n-dimensional object in ¥ and h + +00 as |z| > co on J. Therefore, we can 

identify all the possible integration cycles as an element of the relative homology H,,(¥, X_7; Z) 
for T >> 1, where 

A_p:={x € X | h(x) < -T}. (2.1.13) 


In order words, we identify all the points in Æ if e” is sufficiently close to zero, and compute the 
n-th homology of ¥/X_r. The homology can be computed in the following way (see Appendix A 
for its detail). Assume that any two critical points are not connected via the gradient flow, 
then the Morse function h diverges to —oo as the flow time goes to oo. The Lefschetz thimble 
J, associated to the critical point p, is defined as the moduli space for endpoints of solutions 
c: R — & of the downward flow 


C 


a {Im T, c}p, (2.1.14) 
with the initial condition c(—0oo) = p,. That is, 
Is = { (0) EX | lim c(t) = Po}. (2.1.15) 
——00 


Since all the critical points of h has the Morse index n, the Lefschetz thimbles J, generates 
the n-th homology H,,(4, 4_7;Z). By construction, the integrand becomes non-oscillatory and 
exponentially convergent on Js. 

Let us rewrite the original oscillatory integral (2.1.1) by using integrals on J. For that 
purpose, we regard the real cycle Y as an element of H,(4,X_7;Z) and express it in terms of 
Lefschetz thimbles as 

Vea > feds (2.1.16) 


o€y 


where no E€ Z. With this expression, the oscillatory integral becomes 


Zr=X no [ dé exp(—S/h). (2.1.17) 


oEX 
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Since Im(S/ħ) is constant on each Lefschetz thimble Jp, each integration on the right hand 
side is expected not to suffer from the sign problem. 

The coefficients n, is still unknown, and thus we must give the formula for them. Those 
coefficients can be obtained by introducing the intersection pairing. Consider the relative 
homology H,,(¥, &T; Z), where 


AT = {x EX | h(x) >T} (2.1.18) 
As we have shown for H,(4, ¥_7;Z) for T > 1, the n-th homology H,(4, 7; Z) is generated 
by 
C= { (0) ex | lim c(t) = Po}. (2.1.19) 
t+00 
By definition, h(z) < h(p,) for z € J, and h(z) > h(p,) for z € K,, and the equality holds only 
at z = po. Therefore, Jo and Ko intersects only at pọ. Since we have assumed that any two 


distinct critical points are not connected via the gradient flow, J, and K, does not intersect 
with each other if o # 7. We find that the intersection pairing is naturally given as 


(Jo, Kr) = dor (2.1.20) 
under some appropriate orientations. Therefore, the coefficients n, can be calculated as 
no = (Y, Ko). (2.1.21) 


That is, no is given by the intersection number between the original integration cycle Y and 
Ko, which gives an extension of the formula (1.3.8) for the example of the Airy integral. We 
now find that that the integration cycle can be decomposed into the sum of Lefschetz thimbles 
without any ambiguities. 


2.2 Path integral on Lefschetz thimbles 


We apply the Picard—Lefschetz theory to path integrals, which method is called the Lefschetz- 
thimble path integral. In this section, we first summarize the way to compute the Lefschetz- 
thimble path integral by formally using the argument in Sec. 2.1. This decomposes an infinite- 
dimensional oscillatory integral into the sum of well-defined, i.e. convergent, path integrals on 
Lefschetz thimbles. On each Lefschetz thimble, we show that the Schwinger—Dyson equation 
holds. This section is adapted from Ref. [76]. 


Application of the Picard—Lefschetz theory to path integrals Let us consider the 
path-integral quantization of a classical action S[x(t)] = f dt L(«,da/dt)(For a Minkowski 
path integral, —S should be regarded as iS in the usual convention). In the path-integral quan- 
tization, the transition amplitude K(a;,ts;x;,t;) is computed by summing up the probability 
amplitude exp (—S|z]/h) over the set of all the possible space-time paths VY = {x : fti, ts] > 
R[x (ti) = Xi, u(t) a Tf}, 


K (x,t, £i ti) = [ Drexp (-32). (2.2.1) 


In general, the action S[z] takes complex values for generic spacetime paths x : [t;, ts] > R, 
and (2.2.1) must be regarded as a limit of multiple oscillatory integrals. 

The formal application of the Picard—Lefschetz theory to (2.2.1) can be summarized as 
follows [64—66]: 
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1. First, we must find all the solutions of the “complex” classical equation of motion, ôS [z] = 
0, in X = {2 : fti ti] > Clz(t) = xi, z(t) = xf}. We denote those solutions by zo (o € X 
is the label). 


2. For each z,, we compute the Lefschetz thimble J, by solving the gradient flow equation 


Oz(t; u) 65 |z(t; u)| 
Ou = ôz(t; u) ), Vea) 


with the boundary conditions z(t;—0oo) = z,(t), and z(t; u) = x; and z(ty,u) = ay. As 
we have shown in Sec. 2.1, Re S — +00 and Im S is a constant along flows. Therefore, 
the path integral on each Lefschetz thimble, 


Ip = faco ‘ete — C| lim z(t,u)= zld (2.2.3) 


uU— — Oo 
shows good convergence. 


3. The path integral (2.2.1) can be decomposed into the sum of the path integral on Lefschetz 
thimbles Js; 


K(£f,tf; £i ti) = S| (Ko, R) Dz exp (2) l (2.2.4) 


oE To 


Here, (Ko, R) is the intersection number of Y and the dual thimbles, 


Ks f(; 0) : ftat >C | z(t, +00) > Oy (2.2.5) 


Since we have only shown that this method is valid for finite-dimensional integrals in Sec. 2.1, the 
above procedure is somewhat formal and must be performed with appropriate regularization. 
We assume that some regularization is chosen, and that Lefschetz thimbles are well-defined as 
an integration cycle. If two distinct saddle points are connected via (2.2.2), the system is on 
the Stokes ray and the second assumption breaks down. However, this breakdown does not 
happen for generic cases because Im S is a conserved quantity. In the following of this section, 
we use these assumptions implicitly. 


Quantum equation of motions on Lefschetz thimbles Let us consider the consequence 
of the quantum equation of motion from the viewpoint of Lefschetz thimbles. The quantum 
equation of motions, or the Dyson—Schwinger equation, is given by 


f Dr a eS] = 9 (2.2.6) 


for any boundary conditions at t = t; and ty. We first show that this equality (2.2.6) holds on 
each Lefschetz thimble, that is, for any o € X 


7 ôS [2] e-s] — 
[ D T 0. (22.7) 
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Let us show a finite-dimensional analogue of (2.2.7) since the path integral is assumed to 


be regularized. $(z) a holomorphic function of z = (z',...,2"), which has critical points zo. 
Then, we have to prove that 
OS 
I arz P92) os _ 9, (2.2.8) 
Is oz’ 
The following notations are useful to prove (2.2.8): let e = (e1,...,e") € C”, then 
Se(z) := S(z +£), Zoe = 1 d™ ge Pel), (2.2.9) 

Since we can find that ay ae 

ae ass f ge ee) (2.2.10) 

Oe? e=0 J. oz’ 


it suffices to prove that Zs is independent of ¢ for sufficiently small e. 
All the saddle points z,,- of Są are clearly given by 


Xone = Zo =E; (2.2.11) 
We can also readily find that the Lefschetz thimble Js e of Ss is given by 
Joe = Jo — € := {z — E|z E€ Jo} (2.2.12) 


as a subset of C”. When ¢ = 0, Js does not intersect K, if o # 7 and it intersects K, 
transversally only at one point. This is an open condition, and thus valid for an open region of 
£ including £ = 0. For |e| < 1, 

(Toes Kr) = for, (2.2.13) 


which implies that J,. = J, as integration cycles. This completes the proof of Zoe = Zoo, 
and we obtain (2.2.8). 

Let us reverse the story. Assume that we have a generating functional Z|j] satisfying the 
Dyson-Schwinger equation, or the quantum equation of motion, 


(-S'[ð/ðj] + hj) Z[j] = 0. (2.2.14) 
Using Lefschetz thimbles 7%, we can construct the general solution of this functional differential 


equation as 
Zijl= > de] Dae Sie, (2.2.15) 
a Te 


with some dọ E€ C. These coefficients d, can be anything if one just requires the Dyson- 
Schwinger equation. In order to relate those generic solutions with the path-integral quantiza- 
tion, we must choose dọ = (Ko, V). 


Flow equations around complex classical solutions It would be useful to consider the 
structure of the perturbative fluctuations around a complex classical solution z,. Let us consider 
a classical and real-time action Z = i f dt [32° — V(z)] as an example; then the gradient flow 


2.2.2) is given b 
(2.2.2) is given by det) (eo LV: )) (2.2.16) 
P DE oe n 
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In order to restrict our consideration to the perturbative computations, we consider its lin- 
earization, z = Zo + Az, around the classical solution zo. It gives a parabolic partial differential 
equation, 


Ot? 


In order to find analytic properties of this linearized equation, we try to construct its solutions 
by separating variables. Using an Ansatz Az(t; u) = ei*/4 exp(Au) f(t) (A > 0), (2.2.17) reduces 
to an eigenvalue equation 


Sanh u) = (2 + V" (Z(t ») Az(t; u). (2.2.17) 


o? = 
-(5 + V” (Z(t ») FOHAT UE): (2.2.18) 
By separating the real and imaginary parts, f(t) = f(t) +ifo(t), (2.2.18) reads 


(- ae i OF RV GD] ) ( E ) =A ( f ) l (2.2.19) 


On the space of smooth functions with the Dirichlet boundary condition, this differential op- 
erator L on the left hand side of (2.2.19) is self-adjoint. Its eigenvalues \; are real and eigen- 
functions associated with different eigenvalues are orthogonal to each other. By assumption, 
A; # 0 for any j due to the non-degeneracy of critical points. Using this information, we can 
show that eigenvalues of L must be paired with opposite sign, +A;. To see this, the self-adjoint 
operator L is represented as L = —[0?+ Re V"(Z,(t))] 03 — Im V"(Z,(t)) cı using Pauli matrices 
ci. Introducing the “chirality” matrix € = ig, we can easily show that 


elLe = —L, (2.2.20) 


where cel = 1. 

It may be useful to remark here that the above property on the complexified fluctuation 
operator holds in general, although we are considering the Minkowski path integral as an 
example. It follows from the fact that the fluctuation operator is given by taking the continuum 
limit of the complex symmetric matrix. 

The orthonormal property of L is not enough in order to perform the perturbative calcula- 
tions around a complex saddle point. Let A, and Àp be eigenvalues of L, and their eigenfunctions 
are denoted by (fa, fa,2) and (fo1, fo.2), respectively. Assuming that |Aa| Æ [Ab], we must show 


[ett =0, (2.2.21) 


where fa = faa + ifa and fp = foi + ife for our purpose. Since the eigenfunctions are 
paired by multiplying €t, (fs2,—fo1) = €'(foa, fo2) is an eigenfunction with the eigenvalue 
—,. Orthogonality of L shows that 


care + fa2fo2) = care — fa2fo1) = 0. (2.2.22) 


This is nothing but (2.2.21). The tangent space of J, at z, is spanned by eigenfunctions with 
positive eigenvalues. In contrast, those with negative eigenvalues span the tangent space of Ko. 
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Now, we construct the tangent space T, 7, of J, at zo explicitly. Let {An}; denote the 
positive eigenvalues of L, and {f,}, the corresponding normalized eigenfunctions. Then, the 
tangent space Tz, J, is given by 


Tada = ` age falt) 


n=1 


an € e}. (2.2.23) 
The behavior of Z[z] on J, up to O(Az?) is given by (Az(t) = >, anei”/*f,,(t)) 


Z[z,+Az] = Tz] +Y = / dt f, (t) (5 i v"(ex(t)) fin(t) 


2 
lee) An ; 
= Tz.]—- >> =a%. (2.2.24) 


This means that the real-time path integral on the Lefschetz thimble Jo is realized as a Wiener 
integration at least for small quantum fluctuations around zo. 

Let me give one notice before closing this section. In quantum and many body physics, 
several low-lying eigenvalues À, may be sufficiently close to zero, e.g. A, ~ O(Ah%™) with 
Qn > 0. For these quasi-zero modes, the Gaussian approximation is not necessarily enough and 
we have to take into account the higher order corrections for those modes. 


2.3 Simple examples in quantum mechanics 


In this section, we compute real-time Feynman kernels using Lefschetz-thimble path integral 
for its demonstration. Three simple examples of quantum mechanics will be considered; free 
particles on a line and on a circle, and a harmonic oscillator. Since the our argument in previous 
sections is general and abstract, it would be useful to see the concrete way to compute Lefschetz 
thimbles of simple examples. It is also interesting to see that the real-time path integral for 
those systems turns out to be constructed as Lebesgue integrals on Lefschetz thimbles. This 
section is adapted from Ref. [76]. 


2.3.1 Free non-relativistic particles 


Free particle on a line First, we consider path integral for a free particle on the real line 
R. Let us solve the fundamental solution Kfree(xp, tp : £i, ti) of the Schrödinger equation, 


2 492 


o 
thy v(z, t) = -y gaVe D) (2.3.1) 


by using the Lefschetz-thimble path integral. 
The classical action of this system in (2.2.1) is given by, 


E i fae : OE (2.3.2) 


This can formally be regarded as a holomorphic functional on ¥ = {2 : [tats] > Clz(ti) = 

£i, z(t) = £f}. The classical solution za of this system is given by solving the Euler-Lagrange 
equation: 

d?z4(t) _ 

dt2 


(2.3.3) 
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Under the Dirichlet boundary condition z(t;) = x; and z(t;) = wy, it is solved as 


t-ti 
za(t) = (zf — rome + Ti, (2.3.4) 


and the classical action is Z[z.] = i(xp — 2;)°/2(ty —t;). This is a real-valued function, and 
then it belongs to the original integration cycle Y of (2.2.1). 

Let us compute Lefschetz thimbles, and, for that purpose, we have to write down the 
gradient flow. The Kahler metric of ¥ can be chosen as 


ds? = I dt 5 (52(t) Q Z(t) + dz(t) @ dz(t)), (2.3.5) 


and then the cotangent space is spanned by ôz(t) and 6z(t). The Kahler form of this metric is 
w = f dtiðz(t) A 6z(t). Let us recall the definition of the Lefschetz thimble J around z(t): It 
is defined as a set of gradient flows emanating from z(t): 


z(t; —00) = za(t), 2) = - (Seno) (2.3.6) 


J= {260 ew 


This is a middle-dimensional cycle in ¥ and can be used as a domain of integration. Let us 
rewrite the path-integral expression of the Feynman kernel of the free particle Kfree as 


T 
KralEntriznti) = | Dzexp _ (2.3.7) 
T 


In the following, we will explicitly find that Za has the intersection number 1. In the real-time 
path integrals, it can generally be shown that the real-valued saddle point always has the unit 
intersection number [64, 65]. 

The Lefschetz thimble J are computed as follows. Let us write z(t; u) = zalt) + Az(t; u), 
where Az —> 0 as u + —oo. The gradient flow equation in (2.3.6) becomes a partial differential 
equation of the parabolic type, 


—Az(t;u) = -i— Az(t; u). (2.3.8) 


The boundary condition is Az(t; oo) = 0, and the Dirichlet boundary condition Az(t;;u) = 
Az(tp;u) = 0 is also required. The basis of solutions Az(t; u) for (2.3.8) is given by 


= t-t 
Az(t;u) = /* exp (ar) sin @ ) , (2.3.9) 


with £ € Zo. As a result, any element of Lefschetz thimbles J can be written as 


z(t) = zalt) + e'"/* S asin (~ = ) , (2.3.10) 


where the coefficients ap are real values. 
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Since the Fourier transformation is unitary, the Jacobian associated with the change of 
integration variables does not appear. Indeed, the induced metric on the Lefschetz thimble 7 
from (2.3.5) is given by 


tf 
ds*| , = `. dapday dt sin (e 
The path-integral measure on the Lefschetz thimble J is proportional to 


f P=N | TIVda (2.3.12) 


120 
with some normalization factor M. Since Z[z] — Z[za] is a real-valued Gaussian functional of 
{a¢}e, the real-time path integral becomes the Wiener integration on the Lefschetz thimble. It 


reads 
VACA 
Dz expT[z| =N |] | FE | exp (2.3.13) 
[ ae 4 7a 2(t¢ — ti)? h 
Since the computation of the normalization factor M is standard, we do not give its details. 
The result is 


t) sin n (ae E) = t=) Dao, (2.3.11) 


=A C40 


1 ra [ze] 
[> exp Z|z] = Piht ta) exp ——. (2.3.14) 

It would be useful to point out the relationship of real- and imaginary-time path integrals 
from the viewpoint of Lefschetz thimbles. The imaginary-time formalism is also a mathemati- 
cally rigorous approach to the path integral [77-79]. By formally replacing the real time t by 
the imaginary time —i et, with 0 < y < }, the classical action Z[z] becomes positive definite 
at y = 0. Therefore, the oath integral at Ge = 0 is well-defined as the Wiener integration. The 
complex time t, matches the real time t at y = 7/2, and we compute Lefschetz thimbles at 
each y to show its connection to the real-time Lefschetz thimble. The path integral (2.2.1) 


becomes 
-ip 1/dr\’ 
K free(@ ¢, —i etr; £i —i et) = [ Drexp -o fates 2: : (2.3.15) 
y h 2 \dty 


It is now clear that this path integral with the imaginary time, y = 0, is well-defined as a 
Wiener integral, and that the Lefschetz thimble at y = 0 is nothing but the original integration 
cycle. In order to compute Lefschetz thimbles of (2.3.15) for general y, we introduce 


Z,[2(ty)] = 2 f dtes Ca (2.3.16) 


as a holomorphic functional of z : [t;, ts] + C with z(t;) = x; and z(t) = £p. Since the classical 
solution is obtained by solving the same Euler-Lagrange equation (2.3.3), it is given by zalto) 
of (2.3.4). By repeating the same procedure in (2.3.8-2.3.10), the elements of the Lefschetz 
thimble J around za(tp) are obtained as 


ti 
z(ty) = za(t,) + el? So ae sin (i). (2.3.17) 
m=i 
£>0 
At p = 7/2, as we remarked above, this is equal to the real-time Lefschetz thimble (2.3.10). 
Therefore, our results (2.3.7) and (2.3.14) have a clear relationship with the imaginary-time 
path integral of free theories by deforming Lefschetz thimbles smoothly as given in (2.3.17). 
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Free particle on a circle Using these results, the Lefschetz-thimble path integral can be 
computed also for quantum mechanics of a free particle on a circle St. This is an interesting 
and simple example to learn about the path integral of quantum mechanics on a non-trivial 
topological space. We solve the Schrödinger equation (2.3.1) under the different boundary 
condition for x € St = R/2nZ. The wave functions obey 


ylz + 2n,t) = e y(x, t). 


Here, 0 is a real parameter. 
The classical action in (2.2.1) of this system has a topological 6-term, 


Tz] =i fat; (%) + sf (2.3.18) 


This is a holomorphic functional of z : R + C/2aZ with the boundary condition z(t;) = x; and 
z(tr) = ay. The topological term distinguishes the first homotopy class 7(R/27Z) of possible 
real-valued paths in the usual path integral. This is the same even after the complexification. 
The real-valued target space is embedded as R/27Z C C/2nZ, and thus the topological term 
distinguishes 7(C/27Z) (= Z) of complex paths. 

Solving the same Euler-Lagrange equation (2.3.3) under the boundary condition z(t;) = x; 
and z(t) = x, in R/2aZ C C/2nZ, the set of solutions {Zc,w(t)}wez are parametrized by the 
winding number w €E Z; 


t—t 
Zalw(t) = (xy + 2rw — Ti) 


7 2.3.19 
E a (P 


Let Ju be the Lefschetz thimble for Zea,w(t), and its elements are calculable as (2.3.10). We can 
compute the path integral (2.2.1) of this system as we have done in (2.3.13), and the result is 


z] oo 
aa D [Prey h T E Le (2.3.20) 


wEZ 


Each classical solution has the intersection number 1, and thus all of them labeled by the 
topological number w must be summed up. 


2.3.2 Harmonic oscillator 


So far, path integrals of free particles are reconsidered by using the path integral on Lefschetz 
thimbles. Next, we compute quantum mechanics of a harmonic oscillator. The classical action 


is given by 
. 1f/dz\* 1, 


and we regard it as a holomorphic functional on ¥. The Euler-Lagrange equation is 


d?z(t) 
a = Ale): (2.3.22) 
Thus, the classical solution is obtained as 
T vi t ‘ 
Zalt) = ope acosiy ty) sin(t — t;) + z; cos(t — t;). (2.3.23) 


sin(tp — t,) 
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Let us assume that ts — t; 4 na for any n € Z>o in order to make (2.3.23) meaningful. Since 
this condition can be satisfied for generic initial and final times, t; and tf, the assumption is 
harmless in the following computations. The classical action is given as 


i 


T[za] = 2sin(t, — t,) [ 


(op +, cost — 1) ea] (2.3.24) 
In order to calculate the Lefschetz thimble, the gradient flow equation is considered, 


ð a eo — 
get w) =i ($ + i) Az(t; u). (2.3.25) 
Here, the complex path z is decomposed as z(t; u) = zalt) + Az(t; u) with boundary conditions 


Az(t; —œ0) = 0 and Az(t;u) = Az(ts,u) = 0. The Lefschetz thimble is spanned by the 
following solutions Az, 


eiT/4 exp (=)= 1 
ee mae 3 z in a (t— ti), (nm < (ty ti)). D 


Let v denote the maximal non-negative integer smaller than (ts — t;)/r. Then, the Lefschetz 
thimble J consists of complex paths, 


z(t) = zalt eu 2 Qe sinz 


Tht- ti), (nT > (ts —t)), 
Azn(t; u) = 


Tl 
en + e ag sin 
el S aen” 
ł=v+1 


A (t — ti), (2.3.27) 
with real-valued coefficients ay. Therefore, the path-integral measure on the Lefschetz thimble 
J is given by 


[>= [Te T dan II et dam = corn f T] vide (2.3.28) 


m=v+1 
Compared with that of the free particle case (2.3.12), the path-integral measure has an extra 
factor exp (-i%). This integer v is called the Maslov—Morse index of the classical trajectory 
(see, e.g., Appendix 11 of Ref. [80]), and this represents the number of turning points during 
the time interval [t;, tf]. 
We can readily compute the Feynman kernel Kh.o. of the harmonic oscillator by using the 
Lefschetz thimble in the following way: 
Ilza] - E 
exp (Fal = K “TL Vines (- (= £) —1 e) 
Zell : 
exp (7 ax 2rih( TEA — ((tp — t;)/m€)? 
Heal (2.3.29) 
ex —i— ]. oO. 
2nih| Siu — t;)| PR 2 
The normalization factor M in this calculation is the same with that for the free particle. In the 


case of the harmonic oscillator, the Maslov-Morse index describes how the Lefschetz thimble 
intersects with the original space of paths. 


Kio: (tz, tf; Ti, ti) 
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2.4 Physical relevance of complex classical solutions 


Since the original path integral is defined as the summation over real spacetime paths, it may 
not be clear what is the physical meaning of complex classical solutions. In this section, we 
demonstrate their possible usefulness by considering the tunneling or decaying process in the 
Minkowski path integral at the semi-classical level. This section is adapted from Ref. [76]. 


2.4.1 Complex classical solutions of the double-well potential 


Tunneling is a characteristic process of quantum mechanics. Many dynamical phenomena are 
related to quantum tunneling; e.g., false vacua decays, bubble nucleation at first order phase 
transitions, domain wall fusions, etc. [81, 82]. Although quantum tunneling can be described 
very well by the WKB analysis of the Schrödinger equation, its understanding based on path 
integrals is not fully developed. This is mainly because any classical solutions cannot connect 
two classically stable regions when the potential barrier exceeds the energy. The purpose of this 
section is to open a new possibility to approach this problem by using the Lefschetz-thimble 
path integral, which shows the physical relevance of complex classical solutions. 


Double-well potential As a model to discuss the tunneling, we consider the double-well 
quantum mechanics. In the real-time formalism, the classical action is 


gj: ifa E (4) e s 2 J (2.4.1) 


The Euler-Lagrange equation of (2.4.1) is 


d?z 5 
with the boundary condition z(t;) = x; and z(t) = wy. Instead of solving (2.4.2) directly, we 
use the conservation law of energy, 


() + (22-1)? = p?, (2.4.3) 


where p(€ C) is a complex parameter. Using the Jacobian elliptic functions, the general 
solutions of (2.4.3) can be obtained as follows: 


z(t) = Jg (v2 c, =) l (2.4.4) 


Here, c is an integration constant to satisfy the boundary condition z(t;) = x;, and k = 
(1 + p)/2p is called an elliptic modulus. Our notation of elliptic functions is based on Chap. 22 
of Ref. [83]. Later, we will find all the possible parameters p so as to satisfy z(t) = xy. 
Before calculating them, let us comment on properties of the double-well potential by ob- 
serving qualitative behaviors of certain classical solutions. The potential term in (2.4.1) has two 
classical minima at z = +1, and thus there are two different kinds of real solutions depending on 
their classical energies p?/2. At high energies p? > 1, the energy exceeds the potential barrier at 
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Figure 2.1: Typical behaviors of real classical solutions. The left and right panels show them 
for p > 1 and p < 1, respectively. (Figures are taken from Ref. [76].) 


the origin z = 0, and thus the oscillation of the particle is given by the left panel of Fig. 2.1. In 
contrast, at low energies p? < 1, the potential traps the particle in one of the potential minima 
and its behavior is shown in the right panel of Fig. 2.1. At low energies, it is impossible to 
connect two potential minima if only the real-valued solutions are considered. There are also 
complex classical solutions, which properties will be discussed in the next subsection. 


Classification of classical solutions By solving the boundary condition z(t;) = x; and 
z(tr) = xy, we calculate the set of elliptic modulus k for classical solutions. Let us recall 
that the Jacobian elliptic function in (2.4.4) has the double half periodicities (see Sec. 22. 4 of 
Ref. [83]): 

sd(z+ 2K(k),k) = sd(z + 2iK(v1 — k?), k) = —sd(z,k). (2.4.5) 


Here, K(k) represents the complete elliptic integral given by K(k) = we dO//1 — k? sin? 0. 
The inverse elliptic function has ambiguities because of the above half-periodicities, and a label 
(n,m) must be introduced in order to remove them. The boundary condition is (k’ = v1 — k?) 


o E af ut ya (VEDE, 
o JCR- a( CD22 i ( kk se) a), (2.4.6) 


and it can be solved as 


2 (n K(k) +im K(vI—®)) = Me 


+ {st = r) (=i "ed? (Ea r) l i (2.4.7) 


It is possible to show that an appropriate subset of Z? gives the set of classical solutions by 
using this relation (2.4.7)'. 


Intersection numbers n, and short-time asymptotic behaviors Since we have identi- 
fied the set of complex classical solutions, we can calculate classical actions of those classical 
solutions. By taking the short-time limit t — t; + +0, the leading asymptotic behaviors of 
the classical action is easily calculable, and they provide useful information for intersection 
numbers n,. 


'The proof is given in Ref. [76]. 
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Let us set x; = xf = O just for simplicity. Indeed, the leading short-time asymptotic 
behavior does not depend on x;, xf. For (n,m) € Z, we must solve the boundary condition 
(2.4.7) to find the elliptic modulus K: 


— : (nk (k) + imK(VT—#)) = so. (2.4.8) 


In the limit (ts — ti) + +0, the right-hand-side of (2.4.8) goes to zero. Therefore, the elliptic 
modulus k converges to 1/./2 to satisfy this equality. We find from this observation that 


2 1 ty — 1; 2 l aAa 
o 2 (aE) l O ((ty ti) E (2.4.9) 


If we adopt the standard choice of the branch cut, n must be a non-negative integer. The 
approximate expression for the classical solution is given by 


(n + im)K(1/V2) . 


tf—ti 


ona) (t)= 


1 
d{ 2(n+i 1 2.4.10 
(2(0 + imp avai, =): (2.4.10) 
According to this formula of short-time behaviors, classical solutions are real-valued if and only 
if nm = 0. Other solutions take complex values. The total classical energy of this solution is 
given by 


2 2(2k2 — 1)? ty —t; 


Pinm) _ 1 29 (ama) (2.4.11) 


Here, Pin) refers the parameter p of the label (n,m), where k = \/(1+p)/(2p). We can 
compute the leading term of the classical action Z in this limit only by evaluating the terms of 
O(1/(t; —t;)*) in the Lagrangian. As a result, we find 


Tinm] = i fae (Zoe - Com wae er a E (2.4.12) 


á 3 (t= t) 


The last expression can be obtained by using a formula in Sec. 23.14 of Ref. [83]. We can use 
these formulas also for the imaginary time by replacing t, t;, and ty by —it, —it;, and —it,. 

By using this expression (2.4.12), we can derive a constraint on the intersection number 
no = (Ko, YV). The real-time classical action Z is purely imaginary for real-valued paths, and 
thus complex classical solutions with Re(Z) > 0 does not contribute to the Lefschetz-thimble 
decomposition of the path integral [64, 65]. Moreover, there are real-valued classical solutions 
when nm = 0, and we can conclude that their intersection numbers are unity, nz = 1. It 
means that there are infinitely many real-valued solutions, which contribute to the short-time 
asymptotic behavior of the path integral. In order to consider about complex classical solutions, 
we assume that n > 0 and m # 0. Using (2.4.12), we find the real part of T as 


8K (1/3)! 
3t —t)3 


Since Re(Z) is positive for complex solutions with nm(n?—m7?) < 0, the intersection numbers no 
of those classical solutions must be zero. On the other hand, complex solutions with n >m > 0 


Rel eee | * — m(n? — m’). (2.4.13) 
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or 0 < n < —m may be able to contribute, although the intersection numbers for them are 
difficult to be determined. In order to find the intersection numbers for them, we must solve 
the gradient flow equations. 

It is an open problem to compute undetermined coefficients n,. Moreover, the intersection 
number n, must be defined with an appropriate regularization, since we do not know the notion 
of the intersection between two infinite-dimensional manifolds. This problem always happens 
for interacting quantum field theories. Let N be the number of lattice sites when the lattice 
regularization is used, then the number of complex classical solutions is 


(deg V)” = (degV — 1)”, (2.4.14) 


where deg V denote the degree of the potential term V. Except for free field theories with 
degV = 2, the number of classical solutions grows exponentially with the lattice size N. 
Therefore, it is very difficult to solve the theory completely using the Lefschetz-thimble method. 
Nevertheless, this method sheds a new light on properties of the path integral. In this thesis, 
we will see that the structure of the sign problem is elucidated thanks to this technique. 


2.4.2 Quantum tunneling in the real-time path integral 


Let us try to describe the real-time tunneling phenomena by using the path integral on Lefschetz 
thimbles. Although our discussion in this part is incomplete because of limitation of our current 
knowledge on Lefschetz thimbles, this trial will make clear the physical relevance of complex 
classical solutions. Let us first review the consequence of quantum tunneling in the double-well 
potential. The ground and first excited states, Wo and Yı, are parity even and odd, respectively. 
In order to realize localized wave-functions in the semiclassical limit, we must consider the 
following superpositions, Y4 := (Yo + 1)/V2. In the limit A > 0, 


|ys(x)|? > O(a F 1). (2.4.15) 


These states w+ fail to be eigenstates of the Hamiltonian: A particle localized on a classi- 
cal minimum transits to another classical minimum. This is the quantum tunneling, and its 
oscillation is characterized by the difference of two energies, E1 — Ey ~ e~5°/" (So > 0). 

In order to construct a classical solution describing the quantum tunneling, let us evaluate 
the long-time asymptotic behaviors of complex classical solutions. This would be important 
as a first step to develop our physical insight to describe the quantum tunneling using the 
real-time path integral. 

For our purpose, it would be useful to start from instantons of the imaginary-time path inte- 
gral, which play an important role to study the quantum tunneling. In Fig. 2.2 (a), instanton- 


like solutions in the imaginary-time formalism are shown for the case x; = —1, x7 = 1 at the 
finite-time interval tp —t; = —10i. The two left figures, (n,m) = (0,0) and (1,0), show classical 
solutions with the one-instanton classical action Zọ = —4/3 in the limit ty — t; + —ioo. The 


right figure, (n,m) = (3,0), shows a classical solution with Z ~ 3Z9(= —4), and thus it seems 
to be naturally interpreted as an instanton—anti-instanton—instanton process. Let us study an 
instanton by setting (n,m) = (0,0); then the boundary condition (2.4.7) becomes 


i 22 — 1 2k? — 1 I T JeF 
-`T = -1 a | sa . (2.4.1 
2 7 ( DK2(1 — k2)’ \ tag (Vi = Cee 
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(a) Instanton-like solutions in the imaginary time 


(n,m)=(1,0) (n,m)=(3,0) 
(n,m)=(0,0) 


(b) Corresponding classical solutions in the real time 


Figure 2.2: (a) Instanton-like solutions in the imaginary-time formalism. They correspond to 
(n,m) = (0,0), (1,0), and (3,0), respectively, for x; = —1, xş = 1, and tf — ti = —10i in 
(2.4.7). (b) Corresponding solutions in the real-time formalism. Especially, the solution with 
(n,m) = (0,0) is very close to the unstable classical solution, and it seems to be natural to 
interpret them as sphalerons. (Figures are taken from Ref. [76].) 


We rewrite the right hand side into a suitable form to study an asymptotic behavior in p > 0 
from the side Im p > 0 by using the formula given in Sec. 22. 17 of Ref. [83]. In the limit 
T — +00, we can solve the asymptotic behavior of pio) = 8ie~”, and obtain the well-known 
one-instanton solution (t = —ir): 


l—p 
2(0,0)(T) = y1 — poo) sn (v 1+ p(0,0)T; | ~ tanh r. (2.4.17) 


1+ poo) 


By doing the same analysis for (n,m) = (1,0), we find that another instanton-like solution 
behaves as 


T 1 8 T, 1 — pao) 
~ 1— 1 =l EA 
21,0)(7) P(1,0) S0 ( + Pio) (- rao" P10) a V 1+ pao 


~ tanh (7 n Ti) l (2.4.18) 


2 


with pao) = —8e7T. This explains that both solutions has the same action, Zọ = —4/3, in the 
limit T > oo. 

In Fig. 2.2 (b), we show the corresponding classical solutions in the real-time formalism. 
These processes have sufficient energies, p?/2 ~ 1/2, to overcome the potential barrier clas- 
sically. This always happen when the label (n,m) is fixed and the limit tp — ti — +00 is 
considered: Since the label (n,m) designates the number of oscillation within the time interval 
(ts —t,), the classical solution must stop almost at the top of the potential barrier. They might 
be related to unstable classical solutions, called sphalerons, in the long time limit tf — t; — 00, 
and then these processes are not suppressed in the semi-classical limit [84-86]. 
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Figure 2.3: Highly-oscillatory complex solutions. In spite of their oscillatory nature, actions of 
these solutions can be comparable with multi-instanton actions. Indeed, we obtained p(31,30) = 
0.427 + 0.155i and Z(31,30) = —1.072 + 0.007i for (n,m) = (31,30) with (tp — t:i) = 100, and 
(52,50) = 0.528 +0.185i and Z(52,50) = —2.892 + 0.092i for (n,m) = (52,50) with (tp — t;) = 172, 
in these examples. (Figures are taken from Ref. [76].) 


We expect that there exist complex classical solutions with p ~ 0 because the quantum 
tunneling is a low-energy phenomenon [76, 87]. Indeed, if we perform the Wick rotation T + it 
in (2.4.17) or (2.4.18), we can find classical solutions with p? = 0. Although they do not satisfy 
the boundary condition x; = —1 and xy = 1, it is intriguing to notice that those solutions 
oscillate with the period m. This observation may suggest that the quantum tunneling in our 
formulation is described by complex solutions with large labels (n,m). Let us set the labels 
n, m of the order of tf — t; (see Fig. 2.3). Such solutions can have finite classical actions 
with non-vanishing negative real parts, since they oscillate in the complexified configuration 
space (see Ref. [87] for detailed discussion on this property). Indeed, if n, of such solutions are 
nonzero, then their classical actions have negative real parts and those transition amplitudes 
are exponentially suppressed in the limit A — 0. 

Therefore, it is very promising to conclude that infinitely many complex solutions including 
large (n,m) are significant in the real-time path-integral description of the quantum oscillation. 
This insight can be checked by computing intersection numbers n, of those complex solutions z,, 
and then we expect that those intersection numbers n, must be non-zero. Similar observation 
has been done also in the context of the particle production [88-90], and we will see its great 
importance to study the sign problem in Sec. 4. 


2.5 Brief summary 


In the present Chapter, we reviewed the foundation of the Picard—Lefschetz theory in order to 
apply it to the evaluation of oscillatory multiple integrals in Sec. 2.1 [62-66]. This naturally 
extends the idea of steepest descent methods (see Sec. 1.3 in the introduction). If the number 
of integration variables is one, we can easily designate the steepest descent path by using 
the stationary phase condition. When the number becomes larger, however, that condition is 
insufficient to pick up half-dimensional integration cycles. Mathematically, we can overcome 
this difficulty by using the Morse theory, and the steepest descent cycles can be found by solving 
the gradient flow in the complex phase space. Those steepest descent integration cycles are 
called Lefschetz thimbles, and they form a basis of a linear space. 

In Sec. 2.2, we concretely describe the procedure to apply the Picard—Lefschetz theory to 
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the path integral in quantum mechanics following Refs. [64, 65]. After that, we confirm that the 
quantum equation of motion on each Lefschetz thimble. This explicitly shows that Lefschetz 
thimbles form a basis of complex vector space formed by the set of solutions of the Dyson— 
Schwinger equation. We also study the properties of the Gaussian fluctuation around complex 
classical solutions in a direct analytic way. 

It is important and helpful to understand this method by applying it to the simplest exam- 
ples of quantum mechanics. For that purpose, we compute the real-time Feynman kernels of 
free particles and the harmonic oscillator with this method in Sec. 2.3. The Lefschetz-thimble 
method reproduces the well-known results for these simple models, and we can confirm that 
this model works well. It is still interesting to see that we can compute Maslov—Morse index for 
the Feynman kernel of the harmonic oscillator without introducing explicit regulator explicitly. 

In the following chapters, we will see the usefulness of the Lefschetz-thimble method to ob- 
serve phase transitions and to understand the sign problem appearing in the statistical systems. 
In Sec. 2.4, we mention another possible direction for this method to be helpful. The real-time 
tunneling process can be studied using the WKB analysis, but its path-integral derivation is 
not fully understood. If we can understand it using path integral, then that computation is 
available not only in quantum mechanics but also in quantum field theory. As a prototype of 
this study, we consider the real-time tunneling of the double-well system based on the path 
integral on Lefschetz thimbles. We elucidate that the complex solutions play an important role 
there. 
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Chapter 3 


Lefschetz-thimble approach to the 
spontaneous symmetry breaking 


Chiral symmetry breaking is one of the characteristic nonperturbative phenomena of low-energy 
QCD. In this chapter, we study symmetry breaking phenomena for simple zero-dimensional toy 
models from the viewpoint of Lefschetz-thimble path integrals. In Sec. 3.1, the 0-dimensional 
Gross—Neveu model is considered for studying Zə chiral symmetry breaking, and we directly 
relate the Lefschetz-thimble decomposition and the Lee-Yang zeros in the study of phase tran- 
sitions. We develop a computational method to analyze the O(n)-symmetric systems using 
Lefschetz thimbles in Sec. 3.2. This formalism is effectively used in Sec. 3.3 to study the O(4) 
chiral symmetry breaking of the 0-dimensional Nambu—Jona-Lasinio model. 


3.1 Spontaneous breaking of Zə chiral symmetry 


In this section, we discuss the spontaneous breaking of the discrete symmetry in the theory with 
chiral fermions based on the Lefschetz-thimble approach to the path integral. The discussion 
is adapted from Ref. [91]. 


3.1.1 O-dimensional Gross—Neveu model 


We discuss a prototype of the spontaneous breaking of a discrete symmetry. For this purpose, 
we introduce the 0-dimensional analogue of the Gross-Neveu (GN) model [92]. This is the 
theory of 2-component Grassmann variables 7, and Y, with N colors (a = 1,...,.N), and its 
partition function is defined by 


Zn(G,m) = [aoa exp (a3. (ip +m) va + T ha) ‘). (3.1.1) 


Here, G > 0 is a coupling constant that describes the four-fermion interaction, m is a bare 
fermion mass, and p := 52 pio; is a 2 x 2 matrix-valued constant that mimics the effect of 
nonzero-momentum modes in higher dimensions. In the massless case m = 0, the classical 
action in (3.1.1) is invariant under the Zə chiral transformation defined by Y +> ys. Then, the 
theory has a discrete symmetry, and we will discuss its spontaneous breaking in the large-N 
limit. 
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In order to consider the path integral on Lefschetz thimbles for this system, we perform the 
Hubbard-Stratonovich transformation, 


N 
Zn(G,m) =4/ Z | Tavao exp (>: Palip +m + o)ha — 5] (3.1.2) 
a=1 


= UZ [ew [det (ip + m +0)]™ exp (-3") (3.1.3) 


Here, ø is an auxiliary bosonic field. The equation of motion tells us that ø is nothing but the 
fermion condensate, 


(0) = aa Do (3.1.4) 


The partition function is given by 


N 
Zn(G,m) = y5 fare, (3.1.5) 


where the classical action is 


Slo) = — — log[p? + (o + m)?] (3.1.6) 


#2 
G 
with p? = p}+p3 > 0. Because of the fermion loop, the classical action has a logarithmic term. 
One of important purposes of this section is to demonstrate that the Lefschetz-thimble path 
integral is applicable also for such systems. According to the expression (3.1.5), we can see that 
the saddle-point analysis becomes accurate in the large-N limit. We separate the cases with 
m = 0 and m > 0 in the following, and the chiral limit m = 0 is discussed at first. 


3.1.2 Structures of Lefschetz thimbles in massless case 


Complex saddle points in the massless limit In order to understand the physical prop- 
erties of the model, we assume that G > 0 for a moment and consider the case m = 0: 


2 

S(q) = a — log(p? + 0”). (3.1.7) 
The classical action in the chiral limit takes its minimum at ø = 0 for 0 < G < p? and at o £0 
for G > p?. Because of this difference, the second-order phase transition associated with the 
chiral symmetry breaking occurs at G = p? in the limit N — oo. 

In the following, we discuss this phase transition of the 0-dimensional GN model by using the 
path integral on Lefschetz thimbles. Because of the logarithmic singularity, the complexification 
of o € R is given by z € C \ {tip}. The Morse function Re(S(z)) and the integrand e~%3) 
are well-defined as functions on C \ {ip}, and thus we can compute Lefschetz thimbles as 
one-dimensional manifolds in C \ {tip}. The complex saddle points are obtained as 


OS(z) 2z 2z 
= = — => =0, t/G-—p?. A 
j Oz G p+ am iC a 
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Figure 3.1: Lefschetz thimbles J (black lines) and duals K (dashed red lines) on the z-plane 
for the GN-like model in the chiral limit. The three orange blobs at z = 0 and z4 = iy 1 — G 
are the critical points of S(z), while the two red blobs at z = +i are the points where S(z) 
diverges. The background color scale describes Re S(z). (Figures are taken from Ref. [91].) 


Let us set z4 := +,/G — p?. The three critical points coalesce at G = p?. As a result, Lefschetz 
thimbles become ill-defined at G = p*, and thus we hereafter put G £ p°. 
When G is positive, all the classical actions at the saddle points z = 0, z4 are real, and 


2 2 

p 
oi A. 
aoe (3.1.9) 


$0)\=S6)==14 a — log 


The equality holds if and only if G = p?. Therefore, if z} = +,/G—p? contribute, the 
symmetry is spontaneously broken in the limit N — oo. This can be judged by computing 
Lefschetz thimbles and their duals. 

Before starting its computation, we must notice that Im S(0) = Im $(zz) = 0 for real G 
and p. When this happens, the gradient flow may connect different complex saddle points, and 
the Lefschetz thimbles become ill-defined. This is a typical behavior when the theory is on the 
Stokes ray, as we have explained in Sec. 1.3 [64, 93]. In order to avoid this, we endow G with a 
phase factor e? (0 < |0| < 1)'. This solves the degeneracy of Im S(z) among the critical points. 
We postpone more detailed analysis of Stokes jumps in the complex G-plane until Sec. 3.1.3. 


Lefschetz thimbles in symmetric phase Let us consider the case where the chiral sym- 
metry is unbroken, i.e., G < p?. We compute Lefschetz thimbles and their duals for p = 1 
and G = 0.7e with 6 = +0.1, and the result is shown in Fig. 3.1. There are three Lefschetz 
thimbles J, J.,, and J,_ (solid black lines in Fig. 3.1) and three accompanying duals Ko, K.,, 
and K,_ (dashed red lines in Fig. 3.1). Lefschetz thimbles 7,, and J,_ terminate at branch 
points of the logarithmic term (red blobs), and this behavior is characteristic when fermions 
are included. Even though S(z+) < S(0), only the Lefschetz thimble Jọ contributes since it 
has a nonzero intersection number. Thus, the symmetry is unbroken. 


It is not appropriate to rotate the entire action as S(z) > et? S(z) because it allows the 27-ambiguity in 
the imaginary part of $(z) to induce an ambiguity in the real part of e? S(z). 
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Figure 3.2: Same as Figure 3.1 but with |G| = 1.1. (Figures are taken from Ref. [91].) 

Since 0 = 0 is a Stokes ray, J; and Ko jump as 0 crosses zero. This is a Stokes jumping, 


which is explained analytically in Sec. 1.3. When |G| < p°, Jo does not jump at 6 = 0. The 
Stokes jump does not affect the computation of Zy(G, 0) itself at 0 = 0. 


Lefschetz thimbles in symmetry-broken phase Next, let us set |G] = 1.1 and p? = 1 in 
order to study the symmetry broken phase. A small phase factor is attached to the four-fermion 
coupling G in order to circumvent the Stokes jump. Figure 3.2 shows the behavior of Lefschetz 
thimbles and their duals in this case. The structure of Lefschetz thimbles in Fig. 3.2 totally 
changes in comparison with Fig. 3.1. This difference plays an important role to describe the 
different phases, but the connection is not so straightforward, as we will see in Sec. 3.1.3. 

In Fig. 3.2, all the dual thimbles Ko, K,,, and K,_ (dashed red lines) intersect with R. 
This is the biggest difference between Figs. 3.1 and 3.2, and it indicates that all three Lefschetz 
thimbles J, Jz, and J,_ contribute to Zy(G,0). Indeed it is visually clear that the union of 
the three thimbles is homologically equivalent to R. Since S(z+) < S(0) in the limit 6 = 0 
according to (3.1.9), the behavior of Zy(G,0) in the large-N limit is completely dominated 
by the nontrivial saddle points z+. As a result of the spontaneous breaking of the Zə chiral 
symmetry, fermions obtain a dynamical mass. 

Let us analyze the Stokes jump in Figure 3.2 in more details. The orientation of 7, is fixed 
so that Rez increases along the positive direction, and that of J is defined by the direction 
from z = —ip to z = +ip (p > 0). Then, the Lefschetz-thimble decompositions of the real 
integration cycle R at 0 = OF can be expressed as follows: 


= for 0 = 07, 


Ie tht.  forĝ=0*. (3.1.10) 


There are two remarks on (3.1.10). 


e As we can see in Figs. 3.1 and 3.2, Jo does not jump across the Stokes ray. This fact can 
be intuitively understood: By definition of the gradient flow and (3.1.9), Re S(z)|z, > 
Re S(0) > Re S(z+). Any flows along Jo do not approach z+. 
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e The integral coefficients of J, in (3.1.10) do not jump across 0 = 0. This fact can be 
readily understood by applying the same argument to show that K., do not jump across 
the Stakes ray. This is convincing also from the viewpoint of the asymptotic analysis. 
The asymptotic behavior of Zy(G,0) for G > p? at N > 1 is given by the nontrivial 
saddle points z = z4. Since Zy(G,0) is a continuous function on G, the contributions 
from J,, cannot jump discontinuously in the large-N limit. 


Both arguments have been presented in Sec. 3 of Ref. [64] in the context of bosonic integrals. 
Here, their usefulness is highlighted in a fermionic model. 


3.1.3 Stokes phenomena and Lee—Yang zeros 


Stokes lines So far, we have considered the phase transition of the 0-dimensional GN model, 
and the Stokes jumps for G > 0 are studied. Through those explicit calculations, the Lefschetz 
thimbles turn out to be well-defined if a small complex phase e? (0 < |0| < 1) is attached to the 
four-fermion coupling G. Let us now study the Stokes phenomenon for a generic four-fermion 
coupling G € C.? 

Before starting our analysis, let us remark a technical comment. In the original path-integral 
expression (3.1.5), the integral is convergent only when Re(G) > 0. Only after its computation, 
can we perform the analytic continuation to study Zy(G,0) as an entirely analytic function in 
terms of G. Here, the Lefschetz-thimble decomposition again plays an essential role. Lefschetz 
thimbles connect nicely convergent regions to provide an integration cycle, and thus the path 
integral on Lefschetz thimbles always converges even if the coupling G is complexified. Thus, 
we can study analytic properties of Z(G, 0) while keeping its path-integral expression by using 
Lefschetz thimbles. 

Let us at first study the Stokes jumping in the complex G-plane. A Stokes jump can occur 
if multiple saddle points have the same imaginary part of the classical action. This condition 


reads 
2 2 


0 = Im[S(0) — S(z4)] = Im| =i G -log & , 
Recall that p? > 0 by assumption. This condition can be solved by G € Ryo or by G € 
{ pre’? |r = 2 and — 7 < ġ < r}. We show the union of these curves by blue lines in 
Fig. 3.3, on which the Stokes jump happens. In the shaded area of Fig. 3.3, Jo is the unique 
Lefschetz thimble that contributes to Zy(G,0) as in Fig. 3.1. Outside the shaded area, all the 
Lefschetz thimbles Jo and J, contribute to Zy(G,0) as in Fig. 3.2. 

It must be pointed out that the Stokes jump at the boundary of the shaded area does not 
describe the phase transition directly. In the large-N limit, we find 


p (3.1.11) 


e7090) inside the curve, 
Zn(G,0) l He NSO) 4 He-NSG+) + He-NS52-) outside the curve, i) 
with some coefficients symbolically represented by #. In both cases, Z(G, 0) is entirely dom- 
inated by the same term e~“5) since Re S(0) < Re S(z4) on the boundary. Because of the 
continuity of Zy(G, 0), subdominant contributions can emerge across the Stokes ray but the 
dominant one cannot change discontinuously [93, 96]. 


2A complex four-fermion coupling appears in studies of the 0 vacuum in QCD [94, 95]. 
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Figure 3.3: Stokes lines for the GN-like model with p = 1 and m = 0 (blue lines). The 
global topology of Lefschetz thimbles and their duals changes across the Stokes lines. In the 


shaded area, only Jọ contributes to Zy(G,0). Outside the shaded area, all the three thimbles 
contribute. (The Figure is taken from Ref. [91].) 


Anti-Stokes lines and Lee-Yang zeros Let us now study the phase transition of the chiral 
symmetry breaking in the complex-G plane by applying the path integral on Lefschetz thimbles. 
In order to find the nonzero chiral condensate in the large-N limit, the following two conditions 
must be satisfied simultaneously: 


1. J.,, contribute to Z(G, 0), and 


2. Re S(z+) < Re S(O). 


We have already seen that J}, contribute outside the shaded region in Fig. 3.3. The second 
condition says that 7-,, must give the dominant contributions in the large-N limit in order to 
break the chiral symmetry spontaneously. This is characterized by an anti-Stokes line, at which 
saddle points exchange their dominance. In the present model, the anti-Stokes line is specified 
by 


2 


0=Re [S(0) — S(z4)] = Re ee — log © (3.1.13) 


The anti-Stokes line is shown with the green curve in Fig. 3.4 for p = 1, together with the Stokes 
rays from Fig. 3.3.3 For this simple model, we can neglect the effect of fluctuations around 
saddle points z,, but in general we must compare the real part of “free energies” instead of 
those of classical actions. Outside the green anti-Stokes curve, chiral symmetry is spontaneously 
broken in the limit N — oo. The phase transition along the anti-Stokes curve is generally of 
first order because several saddle points exchange their dominance. The only exception is the 
point G = p?, at which the phase transition is of second order. This happens since the Stokes 
and anti-Stokes curves intersect there, and the physical phase transition occurs at G = p?. 


3The anti-Stokes line actually extends into the interior of the Stokes curve, but this part is not shown in 
Fig, 3.4 because Jz, and Jz_ do not contribute to Zy(G, 0) there. 
+ 
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Figure 3.4: Anti-Stokes line for the GN-like model with p = 1 and m = 0 (green curve), 
overlaid with Lee-Yang zeros for N = 40 (red bullets) and the Stokes line in Figure 3.3 (blue 
curve). (The figure is taken from Ref. [91].) 


Next, we would like to explore a connection between the anti-Stokes line and zeros of the 
partition function. According to the seminal work by Lee and Yang [97, 98], zeros of a finite- 
volume partition function in a complex parameter space, called Lee-Yang zeros, provides rich 
information on the phase transition in the thermodynamic limit. We here give only a brief 
review on this topic (See [99] for a review and [100-107]). Assume that we put a field theory on 
a compact torus so that the partition function Z(j,1) is analytic in terms of a parameter u, such 
as the chemical potential. If the physical parameter is real, the partition function is positive 
definite, and thus it has no zeros for u € R. However, if we allow the chemical potential u 
being complex after computing the partition function, it can have zeros in C \ R, and let us 
denote them by {fia}, i.e., 


Z(u) ~ |] = to). (3.1.14) 


a 


Therefore, the partition function satisfies the following second-order differential equation: 
3 ð ln Z(u) = X 16? (ut — Ha). (3.1.15) 


This is nothing but the Poisson equation in the complex pi-plane, and we can exploit our 
knowledge on electrostatics. In this interpretation, the number density n = Soon In Z is nothing 
but the electric field in the real direction. In the decompactification limit, V —> oo, the set of 
zeros {Ha} typically forms lines in the complex p-plane. If the line of electric charges cuts the 
real y-axis smoothly at ue € R, the electric field „Z has a discontinuous singularity at that 
point and thus the system experiences the first-order phase transition there. If the line ends or 
has a singularity, such as a kink, at some ue € R, 0,,Z() is continuous but singular at u = He. 
This is nothing but the higher-order phase transition. 


In order to see the connection between the anti-Stokes curve and Lee-Yang zeros, we can 
compare them directly. We can evaluate the partition function (3.1.1) or (3.1.5) explicitly to 
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find that 


Zy(G,0) = ey (3) a k! G) . (3.1.16) 


We show the numerical result of Lee-Yang zeros at N = 40 and p = 1 in Fig. 3.4. All Lee-Yang 
zeros are located in the vicinity of the anti-Stokes curve (green curve). As N become larger, 
we can observe that Lee-Yang zeros align on the anti-Stokes curve more and more densely. 
In the large-N limit, they form the anti-Stokes curve. Since the anti-Stokes curve in Fig. 3.4 
has a kink at G = p°, this phase transition is of higher order. By using (3.1.13), we can show 
that the anti-Stokes curve pinches the real axis at angle +7/4. According to a general theory 
of Lee-Yang zeros [99], this behavior characterizes the second-order phase transition with the 
mean-field critical exponent. This is exactly what happens in this model at G = p°. 


3.1.4 Structures of Lefschetz thibles in a massive case 


At the nonzero fermion mass m # 0, it breaks the Zo chiral symmetry explicitly and the “chiral 
condensate” (0) becomes nonzero even for small G > 0. Although a phase transition is absent 
at m Æ 0, it is still interesting to see the behavior of Lefschetz thimbles because of a crossover. 
Let us assume m > 0 in the following. 

We complexify o to z € C in order to use the path integral on Lefschetz thimbles. The 
saddle points of S(z) are obtained as solutions of 
— OS(z) 2z 2(z +m) 


: Oz G p+(z+m)° 


(3.1.17) 


There are three solutions: One of them is always real, and other two solutions are both real or 
form a complex conjugate pair. 

Figure 3.5 (left panel) is the phase diagram of the 0-dimensional Gross-Neveu model. The 
typical shape of S(c) is also shown in each region. Across the blue dashed line in Fig. 3.5 (left 
panel), the number of real saddle points changes. This line does not describe a phase tran- 
sition, but a metastable state appears/disappears across this line. We schematically show in 
Fig. 3.5 (right panel) how complex saddle points move when the line of metastability is crossed 
from below by making the fermion mass m larger. 

We can now elucidate the behavior of Lefschetz thimbles. Let us take two points A and B 
in Fig. 3.5 (left panel) that represent white and blue regions, respectively. In Fig. 3.6 (a), we 
draw the Lefschetz thimbles with p = 1 at the point A. It is interesting to see that no gradient 
flows connect distinct saddle points at m # 0 and G < p?. We do not need to introduce a 
complex factor to G in this case, and the chiral condensate vanishes in the limit m = 0. Just 
as we saw in Fig. 3.1, we find that only one Lefschetz thimble (the real axis, R) contributes 
to the partition function Zy(G,m). In the large-N limit, the saddle point associated with this 
Lefschetz thimble gives the expectation value of the chiral condensate (o) Æ 0. 

Figure 3.6 (b) displays the Lefschetz thimbles at the point B in Fig. 3.5 (left panel) at 
G = 1.5e7?" p = 1, and m = 0.05. When G > p’, the three saddle points lie on the real axis 
with Im S = 0. Therefore, G € R is right on the Stokes ray and we must attach a complex phase 
to G.4 The overall structure of the thimbles is the same as in Fig. 3.2 at m = 0. Among the 


4More generally, a Stokes phenomenon occurs everywhere in the shaded region of Fig. 3.5 (left panel), since 
all critical points are real there. 
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Figure 3.5: Left: Phase diagram of the GN model. The solid red line represents a first-order 
phase transition line. The dashed blue line is a limit of metastability (not a phase transition). 
Concerning the points A and B, see Figs. 3.6 (a) and (b). Right: The behavior of saddle 
points of $(z) when m is increased at fixed G > p°. As we move out of the blue region in the 
left panel upward, two of the three saddles on the real axis merge and then migrate into the 
complex plane. (Figures are taken from Ref. [91].) 
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Figure 3.6: Same as Fig. 3.1 but with G = 0.7 and m = 0.05 (a), and at G = 1.5e-°" and 
m = 0.05 (b). (Figures are taken from Ref. [91].) 


three saddle points, the right-most one has the smallest Re S and hence governs the partition 
function and condensate at N > 1. 


Between G = 0.7 and G = 1.5, the number of contributing Lefschetz thimbles changes from 
one to three. This change happens when the four-fermion coupling G traverses the blue dashed 
boundary in Fig. 3.5 (left panel). This is not a phase transition since the symmetry is always 
broken. From the viewpoint of asymptotic analysis, this is because the dominant saddle point 
always lies on Rx» and moves smoothly in terms of G. Although the Stokes line crosses the 
real axis in the complex G plane, the anti-Stokes line does not if m # 0. 
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3.2 Lefschetz thimbles with continuous symmetry 


In order to study the phase transition associated with continuous symmetry breaking, we must 
develop a formalism to compute the Lefschetz thimble if continuous symmetry exists. If there 
exists a continuous symmetry, the classical solutions of the equation of motion are degener- 
ate in general, which forbids us naive application of the Picard—Lefschetz theory described in 
Sec. 2.1. In Sec. 3.2.1, we give a brief review about the way to compute Lefschetz thimbles 
with continuous symmetry following Refs. [64, 65]. Applying this knowledge, a general way to 
compute Lefschetz thimbles with small explicit breaking terms is developed in Sec. 3.2.2. This 
extension is necessary to study the phase transition associated with the spontaneous continu- 
ous symmetry breaking, and we study its application to the Nambu—Jona-Lasinio-like model 
in Sec. 3.3. This section is adapted from Ref. [108]. 


3.2.1 Exact continuous symmetry 


In this subsection, we review applications of Picard—Lefschetz theory to oscillatory integrations 
with unbroken continuous symmetries [64, 65]. For simplicity, we concretely consider a zero- 
dimensional O(n) sigma model. 


Let o = (0a)a € R”, and then consider the Lie group action O(n) œ R” given by the 
defining representation. Using real anti-symmetric matrices €, its infinitesimal transformation 
reads 


Ons = Eabb (3.2.1) 


with €gp = —Eba € R. We consider a classical action function Sp : R” — R, which is invariant 
under O(n), and study the property of integration, 


Zo(h) = f do exp(—Solo]/h), (3.2.2) 


where d"o = do, A --- A don is the Lebesgue measure. 

Before trying Lefschetz-thimble methods to evaluate (3.2.2), let us consider the result of 
O(n) symmetry directly. The O(n) symmetry of the integrand and integration measure tells 
us that we can integrate radial and angular directions separately, R” \ {0} ~ Rso x S”71. By 
integrating out the angular directions S”~', we just need to evaluate the following integration 


2g”/2 


Zo(h) = ROA) a dao" exp(—So(c)/h), (3.2.3) 


with So(a) = Soloê] (ê is an arbitrary unit vector), instead of the original one (3.2.2). This 
result can be interpreted from the viewpoint of the reduction of integration cycles: The quotient 
space (R” \ {0})/O(n) is identified with the positive real axis Ryo(= (0,00)). In the rest of this 
section, let us review a general way to derive (3.2.3) without taking quotients of the original 
integration cycle at first [64]. This consideration also provides a basic setup for taking into 
account the effect of explicit breaking terms. 


“=> 


We take the complexification C” of R”, and introduce the complex coordinate Ea = Ca +ina. 
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The complexified space C” canonically admits the symplectic structurež, 
w= — 5A, A dé, = dna A doa . (3.2.4) 
We exploit the fact that a gradient flow for the Lefschetz thimble can be viewed as a Hamiltonian 


flow of classical mechanics [64] (see also Sec. 2.1). Namely, we regard {o,}, as canonical 
coordinates and {7,}, as conjugate momenta, and consider the Poisson bracket 


Of Og Of Og 
= — . 3.2.5 
{figs OTa ONa ANa Oe ( ) 
Let us consider the Hamilton mechanics with the Hamiltonian 


then its equation of motion with respect to “time” t, e = { Ho, ĉa}, is nothing but Morse’s 
downward flow equation, 


déa o 
T (s); (3.2.7) 
Along the flow, the Hamiltonian Họ = Im Sof] as well as other constants of motion are 
conserved, which plays a pivotal role in the following discussions. 

The group action (3.2.1) must also be extended. The real orthogonal group O(n) is complex- 
ified to the complex orthogonal group O(n, C) = {A € GL(n, C)|AAT = idp}. The complexified 
group action O(n, C) œ C” is given by 


Oza == Eabêb, (3.2.8) 


with Eab = —Eba € C. This extension makes clear that the O(n) symmetry of the original action 
Slo] is extended to the O(n,C) symmetry of S|]. Under this O(n, C) transformation, the 
symplectic structure (3.2.4) transforms as 


zw = Im(Equ) déa A dép, (3.2.9) 


and thus it is invariant if and only if the infinitesimal parameter Zap is real. Therefore, the 
symmetry of this Hamilton system is O(n), although the Hamiltonian Ho is invariant under 
O(n, C). Noether charges of the O(n) symmetry are 


JablÊ, £) = —Na% + 1bOa = Leg = bbq); (3.2.10) 


and a collection of these Noether charges is called the momentum map, denoted by Ho(n) (See 
Appendix 5 of Ref. [80] for example). It is important to notice that join) = 0 on the original 
integration cycle R”. Since Lefschetz thimbles contributing to the original integration must be 
connected with R” via a flow [64], such Lefschetz thimbles must be connected to the region 
Ho(n) = 0 through O(n, C) transformations. 


5 As we have already seen in Sec. 2.1, the complexified space must be chosen to admit a Kahler metric, which 
is a closed, non-degenerate two-form and compatible with the complex structure. The important point is that 
the Kahler metric uniquely provides the symplectic form w in (3.2.4), which defines the Poisson bracket. 
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Figure 3.7: Downward flows of the reduced Hamilton system (pọ = 0) in the rp, plane with 
a = 0.1. Two points with (r,p,) = (0,0) and (1,0) are critical points of this system. Red 
solid lines represent Lefschetz thimbles J, i.e., downward flows emanating from critical points. 
Green dashed lines are their homological duals KX, which are characterized by downward flows 
getting sucked into critical points. Blue arrows show Hamiltonian vector fields. (The figure is 
taken from Ref. [108].) 


Let us consider the simplest toy model: 


Solé] = FE- 1), (3.2.11) 


with €? = €?+---+€?. Here, a should be formally regarded as an infinitesimally small positive 
constant, œ = 0*, which makes Lefschetz thimbles well-defined as integration cycles. Critical 
points of Sof] consist of two sets: One of them is the origin O = {€ = 0}, and the other one 
is a complex quadric, 


Qr* = {€ EC" | €? = 1}. (92.19) 


Complexified critical orbits can be nicely parametrized with slow variables. For each symmetry, 
there exists two degenerate directions around a critical point: One of them is given by the 
symmetry transformation, and another one by the corresponding Noether charge. This provides 
the easiest way to find the set of zero modes, or slow variables, of the flow equation (3.2.7). 
Mathematically, this implies the fact that Q”! can be identified with the cotangent bundle 
T*S"—' of the hypersphere S"~! = {ø € R” | o? = 1}. We will see this explicitly in Sec. 3.2.2 
for n = 2. 

Since the critical point € = 0 is nondegenerate, we can compute its Lefschetz thimble in 
a usual way. Thus, we only explain how to construct the Lefschetz thimble of €? = 1 in the 
following. Let us introduce the polar coordinate to emphasize the rotational symmetry by 
using the (point) canonical transformation: r = |ø| is a radial coordinate, p, = (a - 7)/|o| 
is its conjugate momentum, and pz; = >> „<p J& represents the total angular momentum. The 
O(n, C)-invariant variable €? can be written using r, pr, and pg as 


& = (a +in) = (07? — n’) + 2io- 


p2 
. (r+ | + 2irp,. (3.2.13) 


r2 


I 
4 
| 
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The Hamiltonian Ho of this O(n)-invariant system becomes 


: 2 2 
H = Z9 G == (0? + 2) -arz 
4 i 
p2 
+ COS @ rpr (e —1)- (r + z) l (3.2.14) 
r 


Because of the conservation law of Ho(n), dynamics of the rp, direction is completely determined 
once på is fixed. The result for the case pj = 0 is shown in Fig. 3.7, and this is nothing but 
the downward flow for the integral (3.2.3) by regarding r + ip, as a complexified variable of ø: 
Hoļpp=0 = Im(So(r + ip,)) with So(o) in (3.2.3). The red and green dashed lines in Fig. 3.7 
show downward flows emanating from and getting sucked into a critical point at r+ ip, = 1, 
respectively, which are given by 


Ho(r, Pr) |pe—0 = Ho(r = 1, pr = 0) |p)=0- (3.2.15) 


Once the radial motion is totally determined in this way, the angular motion is automatically 
determined, and thus we can embed Fig. 3.7 into C”. In order to obtain an n-dimensional 
convergent integration cycle 7,24" and its dual KT, we must pick up other (n — 1) directions 
from the 2(n — 1)-dimensional critical orbit Q’~!(= T*S"~). 

The convenient procedure to construct the Lefschetz thimble, which is proposed in [64], 
is to rotate a red solid line of Fig. 3.7 by using O(n) symmetry. As a result, the Lefschetz 
thimble JST for r = 1 can be determined as the direct product of $"~' = {o? = 1} and the 
red solid line at r = 1 of Fig. 3.7: F224" ~ R x S”! as manifolds. Its dual K224" is given by 
rotating the green dashed line at r = 1 in Fig. 3.7 using the imaginary direction of O(n, ©). 
This procedure ensures that J and K intersect only at a single point. The reader may suspect 
that this construction contains some ambiguities because we pick up (n—1) directions for JET 
by hand, but we can argue that this is the unique choice up to continuous deformations [64]. 
In Sec. 3.2.2, we will see that this special choice of directions in T*S”~' is consistent with the 
downward flows with symmetry-breaking perturbations. 


3.2.2 Continuous symmetry with small explicit breaking 


In this section, we add a small symmetry-breaking term ¢AS, which lifts degeneracies of critical 
points and allows us to use the familiar formulation of Picard—Lefschetz theory. Global behav- 
iors of downward flows are scrutinized, and use of approximate symmetries plays an essential 
role for computation of Lefschetz thimbles. This study can become important in studying 
spontaneous symmetry breaking using Lefschetz thimbles. 


O(2) sigma model We again consider the O(2) linear sigma model with the action (3.2.11), 
and add a symmetry-breaking term, 


cAS[o] = ceo, (3.2.16) 


with 0 < £ <1. We formally set a = 0* to simplify our argument. We study properties of the 
path integral 


Zh) = A doidoe exp [— (So[o] + cAS|a]) /hj , (3.2.17) 
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Figure 3.8: Flows projected on the critical orbit T*S! induced by the perturbation cAH at 
a = 0.1. Black circles at 6 = 0,7 correspond to critical points in the perturbed system. Red 
solid and green dashed curves are downward flows projected onto T*S', which emanate from 
and get sucked into critical points, respectively. Blue arrows show Hamiltonian vector fields. 
(The figure is taken from Ref. [108).) 


using Lefschetz-thimble techniques. We also discuss the relation to the previous formulation in 
the limit € — +0 at the end of this section. 

Since there is an approximate O(2) symmetry, the canonical transformation to the polar 
coordinate is again useful: 


—1 92 
r= 4/0 +08, 0 = tant , 
ai 


onı + GaN 
Pr = fae Po 5 N20 a 
oi + OD 


(3.2.18) 


Here, pg is nothing but the Noether charge Jj in (3.2.10), and w = dp, A dr + dpo A dé. Since 
E? = r? + 2irp, — (p2 + p3/r?), the critical condition €? = 1 can be solved as p, = 0 and 
r? — (pg/r)? = 1. Now, the O(2,C) critical orbit (3.2.12) is explicitly represented as 


_ ( rcosé — i(pe/r) sind r= 4/14/14 4p)/2, 
fe 7 ( rsin 0 + i(pe/r) cos 8 ) y \. (3.2.19) 


po € R,0 € (0, 27) 
As we mentioned, this has a nice parametrization via slow variables 0 and pg, which implies 
that Q! ~ T* S1. 
The Hamiltonian of this system is H = Hyp +eAH = Im(Solé] + cAS|€]). The unperturbed 
Hamiltonian is given in (3.2.14), and the perturbation is given by 


AH = sing r cosl + cosa (pr cos 6 — sin 0) (3.2.20) 
r 


This solves degeneracies of critical orbits, and validates Morse theory. However, naive compu- 
tation of flows can easily fail, because of the remnant of symmetry. After giving a proposal for 
circumventing this technical difficulty, we show its origin by scrutinizing properties of downward 
flows through concrete computations. 

Before starting explicit computations, we mention an efficient way to compute Lefschetz 
thimbles and comment on behaviors of downward flows. We advocate that the downward 
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flow on the nearly critical orbit T*S! should be calculated at first: The perturbation AH 
breaks the O(2, C) symmetry of S[€], and it induces slow downward flows. The slow downward 
flow can be projected onto T*S', which is shown in Fig. 3.8. In order to get Fig. 3.8, the 
Hamilton mechanics in terms of 0 and pg is solved under the critical condition p, = 0 and 


p= ya + ,/1 + 4p? )/2. This step turns out to be important because Lefschetz thimbles in 
this system are two-dimensional cycles with one nearly flat direction: To create such integration 
cycles, downward flows must first travel slowly along red solid curves on T*St in Fig. 3.8 [up 
to O(e)]. After that, flows branch into radial directions, as will be confirmed below. 

The critical condition 0S/0€, = 0 gives 


@+a-0( 2 )+(5)=(5). (3.2.21) 


Critical points of this system consist only of three points: (&,) = (€,0), (£1 — ¢/2,0) up 
to O(c), and two of them (+1 — ¢/2,0) correspond to black blobs in Fig. 3.8 (œa = 0.1 in this 
figure). In the polar coordinate (3.2.18), these solutions are (r,@) = (e,0), (1 — €/2,0), and 
(1 + €/2,7) with p, = pp = 0. Among these critical points, Re(So + ¢AS) takes the smallest 
value at (r,@) = (1+ ¢/2,7), and the largest value at (r,@) = (€,0). Because of the symmetry- 
breaking term, the value at (r,@) = (1 — ¢/2,0) is lifted up from the smallest one roughly by 
2e cos a(= 2e). 

The Hamiltonian H takes 0, sina, and — sina at (r,@) = (e, 0), (1—e/2,0), and (1+¢/2,7), 
respectively, up to O(e°), and thus flows cannot connect distinct critical points. At a = 07, 
however, these differences are infinitesimally small: The Lefschetz thimble J1+-/2,.) around 
(r,0) = (1+ ¢/2,7) can pass infinitesimally close to other critical points. 

The set of Hamilton equations (exactly at a = 0) is 


A =r G — a — ) — 3rp? + ecos8, (3.2.22a) 
Pr = (1 — 3r?)p, + p? — = e sin 0, (3.2.22b) 
“ 2 A — eat (3.2.22c) 
“be =€ (p, sin ô + - cos 0) . (3.2.22d) 


One can find a set of special solutions of this Hamiltonian system such that p, = 0 and 
po = 0, which indeed automatically solves (3.2.22b) and (3.2.22d). Moreover, Figs. 3.7 and 3.8 
imply that these are nothing but the conditions for a part of J142/2,x). Solutions of (3.2.22) are 
shown in Fig. 3.9 (a) in the original coordinate system (01, 02) = r(cos 0, sin 0). Let us consider 
flows starting from (r,@) = (1+ ¢/2,7). Near the critical point, one solution flows along the 
radial direction in the same way as in the symmetric case, but there is another solution along 
the @ direction. Indeed, setting r ~ 1 in (3.2.22c) gives 


O(t) = +2tan7'(exp(—et)). (3.2.23) 


According to (3.2.23), there is a downward flow rotating along a nearly flat direction with a 
time scale 1/e. For consistency with (3.2.22a), the radial direction should be fine-tuned so that 
r(t) ~ 1—(e/2) tanh(et). After traveling along the critical orbit T*.S" in this way, flows branch 
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Figure 3.9: (a) Behaviors of downward flows (3.2.22) in the o,09-plane when p, = pọ = 0 and 
a = 0. Black blobs are critical points of the symmetry-broken system. Red arrowed curves are 
typical solutions starting from the global minimum (r,@) = (1+€/2, 7r) withO < € « 1 (e = 0.1 
in this figure). Before branching into the r direction, every flow from the global minimum moves 
slowly along the critical orbit r = 1. (b) Behaviors of solutions to the differential equation 
(3.2.22) in the o,-plane with 0 = p, = 0 when a = 0. Red arrowed curves show typical 
solutions starting from (r,@) = (1 — ¢/2,0) with € = 0.1. (Figures are taken from Ref. [108].) 


into radial directions, as shown in Fig. 3.9 (a): Drawing flow lines starting from a critical point 
is thus a difficult task, because we must designate the way of branching after traveling the nearly 
critical orbit just by tuning the initial condition. Since we already know the approximate flow 
on the critical orbit in Fig. 3.8; however, we can start and revert a flow from a point of branch 
in drawing Fig. 3.9 (a). This procedure does not require fine-tuning of initial conditions. All 
we have to confirm is that the reverted flow is drawn into the critical point. 

Next, let us calculate the Lefschetz thimble Ja-e/2,0) for (r,0) = (1 — ¢/2,0). For that 
purpose, we consider another set of special solutions of (3.2.22) by putting 0 = 0 and p, = 0, 
which gives the red solid line around 0 = 0 in Fig. 3.8 (but at a = 0). This condition solves 
(3.2.22b) and (3.2.22c) automatically. Behaviors of the downward flow in this plane are shown 
in Fig. 3.9, and, as we expected, downward flows go along the hyperbola o? — 73 = 1 at first, 


/ 2 
which is characterized by pg: Put r(pọ) = ie lp and then (3.2.22d) 


2r(p9(t)) — tanh *(1/r(pe(t))) = et. (3.2.24) 


After this slow motion, they branch into the r direction in order to span a two-dimensional 
cycle J1—</2,0). In drawing red flow lines in Fig. 3.9 (b), we start and revert the flow from a 
point of branch instead of a point near the critical point. It would be really hard to control 
global behaviors of the flow if we compute the flow starting from a point near the critical point. 

Now that we have identified local behaviors of Lefschetz thimbles around critical points 
and properties of approximate downward flows on the nearly critical orbit, we can comment 
on global behaviors of Lefschetz thimbles. This enables us to comment also on the relation 
of Lefschetz thimbles in symmetric and non-symmetric cases. The Lefschetz thimble Ja4+e/2,7) 
around the global minimum goes away to infinities in the ņ direction at 6 = 0, 27 because there 
is another Lefschetz thimble Ja-e/2,0) at 0 = 0. See Fig. 3.10 for its schematic description. In 
order to recover the symmetric Lefschetz thimble JF in the limit € — +0, we must sum 
TJate/2n) and Ja-e/2,0) for the cancellation of the 72 direction. 


gives 
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Figure 3.10: Schematic figure for a global structure of Lefschetz thimbles J1+</2,7) and J(1—</2,n) 
at a = 07, which are shown with green and blue surfaces, respectively. In this limit, flows into 
Pr direction of these thimbles are quite small unless r ~ 0, and we totally neglect that direction. 
Three black blobs represent critical points of this system, and arrowed lines show slow motion 
with time scale 1/e along the critical orbit T*S'. (Figures are taken from Ref. [108].) 


What is the fate of each Lefschetz thimble J(14</2,7), J(—</2,0) in the limit € + +0? Let us 
observe the behavior of integration on Ja-e/2,0) in the limit € —> +0: 


f d?z exp [— (So + £A S) /ħ] 
Ja—e/2,0) 


oo Di 2 oo 
~ -if AdX exp (-“G*) iz do exp (-2 cosh 6) . (3.2.25) 


Here, we have set o1 = A cosh ġ and m = Asinh@. The integration of À is convergent, but the 
integration of ¢ is logarithmically divergent in the limit ¢/h —> +0: 


[eee (-2 cosh 6) ~ 2In À (3.2.26) 


This divergent integration is nothing but the integration along the critical orbit r?—(pg/r)? = 1. 
Exactly the same divergence with the opposite sign appears in the integration over J(+e/2,7), 
and we must sum them up to construct a convergent integration cycle in the limit € > 0. As 


already mentioned, this is identical with 7,27" in the homological sense: 


IEY = Ja+e/2,7) F Ja-—e/2,0)- (3.2.27) 


This might be surprising because the number of possible integration cycles changes under 
continuous change of parameters. Physically, this means that the number of solutions of the 
Dyson-Schwinger equation in the symmetric system is smaller than that of symmetry-broken 
systems. 

We need to pick up half dimensions of the critical orbit T*S"~' in order to construct 
the middle-dimensional integration cycle 7,%;". In the previous section, we picked them up 
in a specific way by rotating the Lefschetz thimble of the reduced system using the original 
symmetry O(n). However, we did not explain why we cannot use imaginary directions of 
O(n, C). The above analysis clearly shows the reason we have to give up imaginary directions 
of the complexified symmetry group: Otherwise, integration does not converge. 

We have stuck to the case a = 0* in order to draw Figs. 3.9, but our proposal is expected 
to be useful for generic a. Since there is no Stokes jump in 0 < a < 7, Lefschetz thimbles at 
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generic a are continuous deformations of Fig. 3.10 (in C?). In order to span such cycles, flows 
should still go along T*St up to O(e) as shown in Fig. 3.8, and branch. The dynamics after 
the branch can be solved in a similar way to get Fig. 3.7: €AH is expected to be negligible 
since critical conditions are broken at O(e°). 


Generalization to O(n) sigma models Let us consider the case with general n. Then, the 
original action has O(n) symmetry, and the perturbation AH breaks it into O(n—1) symmetry: 


OG <i ( : aaa 7 ) C O(n). (3.2.28) 


In order to study the effect of AH, we need nothing new. All the computations are reduced 
into the ones in the case of O(2) symmetry as follows. 

Since the O(n — 1) symmetry is not broken at all, let us integrate it out. The O(n — 1)- 
invariant variables are 


01, 09 = 1/03 +--- +02. (3.2.29) 


The original O(n) invariance ensures that the unperturbed action can only depend on o? + 
ol’. Therefore, even after O(n — 1) symmetry is reduced, there still exists a continuous O(2) 
symmetry. Thus, all of the previous argument for the system with explicit symmetry breaking 
holds. After computing downward flows in the reduced system, we can construct Lefschetz 
thimbles by rotating them under O(n — 1). 


3.3 Spontaneous breaking of U(1), chiral symmetry 


In this section, we discuss the spontaneous breaking of the continuous symmetry in the theory 
with chiral fermions based on the Lefschetz-thimble approach. The computational technique 
developed in Sec. 3.2 is applied in our analysis for that purpose. The discussion is adapted 
from Ref. [91]. 


3.3.1 O-dimensional Nambu—Jona-Lasinio model 


Let us consider two zero-dimensional toy models which have continuous chiral symmetries. 
They are zero-dimensional analogues of the Nambu—Jona-Lasinio (NJL) model [45, 46]. The 
first model is defined by 


N 


Zu p= f dda exp (27. (ip + m)ba + AET Ve) + (L Finwe) Y) (3.3.1) 


We use the same symbols and variables as in (3.1.1). In the chiral limit m = 0, the classical 
action is invariant under the U(1), chiral rotation y —> ey and Y => we”. By introducing 
auxiliary fields ø and 7 via the Hubbard-Stratonovich transformation, we get 


7G 


N 
= 7a dodr eNSlon) | (3.3.2) 


N N 
Zu (1) = o [det (ip +m+ot+ iysT) | N exp (Fe it z?) 
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with 
otn? 
S(o,m) = — log [p? + (o + m)? + 7°] 4 G` (3.3.3) 
The U(1)a S O(2) transformation rotates (o,r) as a two-dimensional vector, and it is a sym- 
metry when m = 0. If we set 7 = 0 by using an O(2) rotation, then the present model reduces 
to the GN-like model in Sec. 3.1 up to an O(2)-invariant measure. This perspective will become 
important in analyzing the structure of Lefschetz thimbles later. 
The second model we consider is defined by the partition function 


Zsa = [ eaves (°F, (ip + mja + AET Ue), +E (Lant) 


A=1 a=1 
(3.3.4) 


where ~ and 7 are two-component Grassmann variables with SU (2) flavors and N colors, {74} 
are the Pauli matrices for SU (2) flavors, and the summation over flavor indices is implicitly 
assumed. In the chiral limit m = 0, this model has an exact SU (2)r x SU (2) symmetry, and 
it is explicitly broken to its diagonal component, SU(2)L+r, by nonzero m. Again using the 
Hubbard-Stratonovich transformation, 


N \? N 
Zsu) = (=) a dodr, [det(ip +m+o+ inst aT“)| exp (-Fe + r) ) (3.3.5) 


z (5) f ioira tonto eo -Fo e30 


Under the chiral symmetry SU (2) xX SU (2)r = O(4), (0, T1, T2, T3) rotates as a four-dimensional 
vector. One can rotate any such vectors to (ø, 0,0, 73) by using an unbroken SU (2)L+r = SO(3) 
rotation. The resulting integral over o and 73 is essentially equivalent to the former model 
s(3.3.2) and does not entail a new feature. For this reason we focus on the first model in the 
following. 


3.3.2 Structures of Lefschetz thimbles in the massless case 


Let us consider the chiral limit m = 0 at first. In this case, there exists an exact continuous 
symmetry. Therefore, in order to apply the Lefschetz-thimble method for Zyga), we have to 
decompose the original integration cycle R? according to Sec. 3.2.1. Upon a complexification 
of variables, the action becomes 


Di a2 
S(z,w) = —log(p? + 2? + w°) + a = ‘ (3.3.7) 
whose domain is {(z,w) € C? |2? + w? Æ —p°}. The set of singular points of logarithm, 
{(z,w) € C? |z? +w? = —p?}, forms a two-dimensional surface in C?. It is equal to the 

O(2, C)-orbit of the singular points (z, w) = (+p, 0) of the massless GN-like model. 

The flow equation reads? 
dz 2z 2z dw 2w 2w 

= — d = — : 3.3.8 
dr pe+z2#+u? G and ar p+2+w? G ( ) 


In this section, the definition of the fictitious time 7 has an extra minus sign. 
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It is important to point out that the symmetry of this flow equation is O(2) instead of O(2, C), 
although the classical action is invariant under O(2, C) (see Sec. 3.2.1). 
The set of saddle points is obtained by solving (3.3.8) with 


dz dw 
— = — =Q. 3.3.9 
dr dr ( ) 
Let us put G Æ p? in order to circumvent accidental degeneracy of complex saddle points. The 
set of saddle points in C consists of two connected components: 


Co := {(0,0)} and Ci:= {(z7,w) E C| +w = G-p}. (3.3.10) 


C; is equal to the O(2, C) critical orbit of the saddle points (+ G — p?, 0) in the 0-dimensional 
massless GN model. It crosses the real plane R? if G—p? > 0 but has no intersection otherwise. 

We now construct Lefschetz thimbles J, and their duals Ko (o = 0,1) for saddle points 
Co and Ci. Since Co = {(0,0)} is a non-degenerate saddle point, we can compute J and Ko 
in the standard manner; Re S(z, w) increases in two directions and decreases in the other two 
directions around Cp. Since the gradient flow preserves the O(2) symmetry, we can obtain Jo 
by rotating a Lefschetz thimble in the GN-like model with the O(2) action: 


z cos? —sinĝ\/z' l 
h = 6 = ae cos 0 6) —m<@O<7 and z E Blan} ; (3.3.11) 


The dual thimble Ko can be computed in the same way. 

Since C is the set of degenerate saddle points, we need a special treatment to construct J1. 
The general method given in Ref. [64] is already explained in Sec. 3.2.1, but we try to compute 
it in a more concrete manner. Let us pick up a point (,/G — p?,0) € C1. Two directions along 
Cı give two null eigenvectors of the Hessian matrix of S, and other two directions are non- 
degenerate. For non-degenerate directions, we can construct the flow lines, and the flow line 
emanating from (,/G' — p?,0) can be identified as a Lefschetz thimble J., of the 0-dimensional 
GN model. Recalling that the flow respects the O(2) symmetry, Jı can simply be obtained as 
an O(2) revolution of J,, in Sec. 3.1. We thus conclude that 


z cos@ —sin@\( 2" , 
x =| (‘) = a e -T <0<rT and z E Jalen) (3.3.12) 


Since the Noether charge u is conserved along a flow, u vanishes everywhere on J. 
By our construction, it is clear that behaviors of the Lefschetz thimbles, (3.3.11) and (3.3.12), 
in terms of G € C can be learned from Sec. 3.1 without additional calculation: 


e J and Jı jump on the Stokes lines (blue lines) in Fig. 3.3. 
e For G inside the shaded area of Fig. 3.3, only J contributes to Zy 1). 
e For G outside the shaded area of Fig. 3.3, both Jo and Jı contribute to Zya). 


e For G inside the anti-Stokes line (green line) in Fig. 3.4, % dominates Jı at large N. 
Their dominance is exchanged as G moves out across the anti-Stokes line. 
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Figure 3.11: Snapshots of the action S(o,7) in (3.3.3) for varying m at G = 3 and p= 1. The 
number of saddle points jumps from 3 to 1 at m œ 0.724. (Figures are taken from Ref. [91].) 


e Spontaneous breaking of the O(2) chiral symmetry is a first-order phase transition for 
generic complex couplings G. For physical G € Ryo, this is the second-order phase 
transition at G = p°. 


As a conclusion, the Lefschetz thimbles for the 0-dimensional massless NJL model (3.3.2) can 
be obtained by O(2) revolutions of the Lefschetz thimbles for the 0-dimensional massless GN 
model (3.1.5). Therefore, we can immediately derive many properties of the chiral symmetry 
breaking for the NJL-like model thanks to the knowledge of the GN-like model. 


3.3.3 Massive case 


With nonzero fermion mass m Æ 0, the O(2) chiral symmetry is explicitly broken. The behaviors 
of the NJL-like model becomes qualitatively similar to those of the massive GN-like model. 
Because of the absence of the continuous symmetry, there are only three saddle points {(z,w) = 
(z5,0)|o = 0,+}, where {z,},~0,4 are the three solutions to (3.1.17), and a critical manifold 
(3.3.10) has different meaning. Associated with three non-degenerate saddle points, there are 
three Lefschetz thimbles J, and their dual thimbles Ko. 

We have already discussed the dependence of {z,}, on G and m in Sec. 3.1.4 (see, e.g., 
Fig. 3.5). For sufficiently large fermion mass m, there is only one saddle point in the real space 
(o,r) € R?, and only one Lefschetz thimble contributes to Zy(1) (see Fig. 3.11 (c)). At small 
m or large G, all the three Lefschetz thimbles contribute (see Figs. 3.11 (a) and (b)). In the 
latter case, the Stokes phenomenon occurs among the three thimbles, but it is hard to visualize 
as it occurs in C?. 

Let us try to calculate Lefschetz thimbles for sufficiently small m. This helps us to visualize 
the Stokes jump in the NJL-like model. So far, we have emphasized similarities of the NJL-like 
model to the GN-like model. However, the flow equation of the NJL-like model has a unique 
character: There are quasi-fixed points of the flow equation because of the approximate O(2, C) 
symmetry of S (see Sec. 3.2.2). That is, the flow along the critical orbit C1 in the massless case 
is extremely slow compared with that of other directions. 

For simplicity, we put p = 0 in the following (This simplification in the GN-like model is 
discussed in Ref. [91]). Here we summarize main features of the NJL-like model at p = 0. 


e m=0: The Lefschetz thimble Jọ in (3.3.11) disappears at p = 0 as Talat no longer 
exists. This is because the singularity of logarithm exists at the origin at p = 0. Therefore, 
J, is the unique Lefschetz thimble. From (3.3.12), we find Jı = R? \ {(0,0)}. 


—m+tV/m?+4G 
9 p} 0 =: (z+,0), 


e m >00: There are only two saddle points at (z, w) = ( 
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Figure 3.12: The downward flow (3.3.13) of the NJL-like model in a plane with Im z = Im w = 
0. The background color scale represents Re $(z,w). The red blob at (—m, 0) is a singularity of 
the flow. The orange circle in the left panel and the orange blobs in the right panel are saddle 
points. (Figures are taken from Ref. [91].) 


and C is no longer the set of saddle points. The associated Lefschetz thimbles are 
denoted by J, and 7. Both J, and J contributes to (3.3.2), and J, gives a dominant 
contribution at N > 1 as Re(S(z+,0)) < Re(S(z_,0)). 


Thus, the Lefschetz thimble Jı is split as J} U J by an arbitrarily small m ¥ 0. In order to 
visualize the Stokes jump, we compute Lefschetz thimbles J+ on the Stokes ray G > 0, and 
the quasi-stationary flows are important there. The flow equations at p = 0 read’ 


dz 2(z +m) 2z 
= — 3.3.13 
dr (z+m)?+w? G’ ( a 
w 2 2 
IP nny i (3.3.13b) 


dr (z+m)?+u? G 


It is easy to check that the condition Im z = Im w = 0 is conserved along the flow. That is, if 
we start a flow from a point (z,w) € R?, then the flow stays in R?. In Fig. 3.12, we show a 
sketch of the flow for (z,w) € R? at m = 0 (left panel) and m > 0 (right panel). At m = 0, 
there is a ring of fixed points z? + w? = G (orange line) that is Cı N R? and attracts all flows 
on the plane Jı = R? \ {(0,0)}. 

For m > 0, there are only two saddle points shown with orange blobs in Fig. 3.12 (b). This 
figure corresponds to the left panel of Fig. 3.9. Now there is a slow but non-vanishing flow 
along the orbit z?+ w? ~ G. It is clear from the figure that (z,,0) attracts all the flows on this 
plane, except for the real axis to the left of (—m,0) on which the flow is attracted to (z_,0). 
Since J, is defined as the union of all flows that sink in (24,0), we find that the set 


Ja := R? \ {(z,0) |£ < =m} c C? (3.3.14) 


is contained in J4. 


TIn this section, the definition of the fictitious time 7 has an extra minus sign. 


3.3. Spontaneous breaking of U(1),4 chiral symmetry 53 


(a) m=0(G=2, p=0) (b) m= 0.4 (G =2, p=0) 


Figure 3.13: Same as Fig. 3.12 but with Imz = Rew = 0. The X-shaped white lines are 
points of logarithmic singularity of the action. (Figures are taken from Ref. [91].) 


Let us go back to (3.3.13), then the condition Imz = Rew = 0 is also conserved along 
the flow. Figure 3.13 shows the flow pattern in the (Rez, Im w) plane with Im z = Rew = 0. 
At m = 0 the fixed points of the flow form a hyperbola (Rez)? — (Im w)? = G, which is the 
non-compact orbit of C4 in Sec. 3.3.2. 

At m > 0, a slow but non-vanishing flow appears along the hyperbola (see Fig. 3.13 (b), 
which corresponds to the right panel of Fig. 3.9). According to Fig. 3.13 (b), we can find that 
the Lefschetz thimble J— is given by the wedge-shaped region 


Jg := {(z,w) € C? | Imz=0, Rew =0, Rez +m <Imw < —Rez — m}. (3.3.15) 


In Fig. 3.14, we combines flow plots in Figs. 3.12 (b) and 3.13 (b) at m = 0.4 within a 
three-dimensional subspace of C? specified by Im z = 0, as we did in Fig. 3.10. Since G > 0 
is on the Stokes ray, two saddle points (z+, 0) are connected via a flow and Lefschetz thimbles 
are ill-defined as cycles. Although the singularity of Lefschetz thimbles can be smoothed out 
through complexification of G as Get, it makes the Lefschetz thimbles extend into the whole 
C? space. 

Since the complexification of G impedes our visualization of the Stokes jump, we shall deduce 
the behavior of Lefschetz thimbles as 0 + 0* by exploiting an indirect argument. The sum of 
J, and J- is identical to the original integration cycle, R’, since both J} and J- contribute 
to Zua) with unit coefficients. Since J_ gives a subdominant contribution compared with J4, 
J, changes by an integer multiple of 7 at 0 = 0*. These considerations indicate that we have 


oe eye S i 


` ` at 0=0+. (3.3.16) 
J- = Jg J- =p 

Here, we assume an appropriate orientation of Jg. We find that the Lefschetz-thimble decom- 
position is given by R? = J} — J- at 0 = 07 and by R? = J} + J- at 0 = 0+. It is important 
to notice that the contribution of Jz is exactly canceled between J, and J. This cancellation 
keeps the partition function real-valued and continuous in G, because the integral over Jz is 
purely imaginary (see below). This is one example that highlights the reason why we need sum 
up multiple Lefschetz thimbles in order to obtain a physically correct result. 
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Figure 3.14: The downward flow (3.3.13) in a 3-dimensional space with Imz = 0. The 
parameters are m = 0.4, G = 2 and p = 0 (same as for the right panels in Figs. 3.12 and 3.13). 
(Figures are taken from Ref. [91].) 


It is intriguing to ask what happens near the chiral limit. As m — 0*, we can regard 
that J4 converges to J, = {(z,w) € R*|(z,w) 4 (0,0)} up to a null set. However, according 
to (3.3.16), both J} and J contain an imaginary component Jp, and thus neither of them 
converges to Jı in the chiral limit. The same thing happens for the bosonic model discussed 
in Sec. 3.3.2. A unique feature of the fermionic model in comparison to bosonic models is that 
one of the thimbles (.7_) is the wedge-shaped region because of the logarithmic singularity. 


Near the chiral limit (m « VG), we have observed that a quasi-stationary flow along the 
orbit z?+w? ~ G emerges (see Figs. 3.12 and 3.13). This originates from the O(2, C) symmetry 
of S(z,w) at m = 0. Since O(2, C) is non compact, the non-compact orbit on Jp (see Fig. 3.13) 
extends to infinity and it gives a divergent contribution in the limit m — 0, as we discussed in 
Sec. 3.2.2. More explicitly, we find for the contribution from Jp 


N N _N 2a? 
Z\;, = a dzdw {(z +m) +w?) eee te) (3.3.17) 
- iP arpn f S JD en Blin? +2m cosh +R?) (3.3.18) 
T 0 —oo 
N Np [S _N 2NR 
= aie G | dR R?NHe oR Ko(—m) ; (3.3.19) 
In the second and third lines, we changed the integration variables as z = —m — Rcosh® and 


w = —iRsinh®. Roughly speaking, © parametrizes the hyperbola in Fig. 3.13. Using the 
asymptotic formula Ko(x) = — log 5 + 7) I(x) + O(a”), we conclude that Z Pa diverges as 
~ logm in the chiral limit, m — 0. Since both J} and J include Jg with the opposite sign 
(see (3.3.16)), this infrared divergence on each thimble exactly cancels out when we sum up 
both contributions. We suspect that this is a generic phenomenon that occurs when we handle 
the spontaneous breaking of continuous symmetries by the path integral on Lefschetz thimbles. 
It is an important check of future studies to compute the behavior of Lefschetz thimbles for 
field-theoretical models in higher dimensions. 
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3.4 Brief summary 


In this chapter, we discussed the chiral symmetry breaking from the viewpoint of the Lefschetz- 
thimble approach. Since the spontaneous breaking of the chiral symmetry is an essential in- 
gredient to understand low-energy dynamics of QCD, this study gives a part of the theoretical 
foundation to study low-temperature and finite-density QCD using Lefschetz thimbles. 

In Sec. 3.1, we first study the case where the chiral symmetry is given by a discrete group 
Zə using a 0-dimensional Gross—Neveu model. After introducing a auxiliary bosonic field for 
the chiral condensate, we apply the Lefschetz-thimble method for evaluating the bosonic in- 
tegral. This is the first application of the Lefschetz-thimble method when the effective action 
has a contribution from the fermionic determinant. The Lefschetz-thimble decomposition of 
the integration cycle provides a nice setup to study the analytic structure of the partition func- 
tion because we can complexify the four-fermion coupling before performing the path integral. 
Thanks to this property, we can discuss the phase transition of this model in the large-N limit 
from the viewpoint of anti-Stokes rays and Lee—Yang zeros. 

In QCD, the chiral symmetry is a continuous and approximate symmetry. Indeed, pions 
are pseudo Nambu-Goldstone bosons of the symmetry breaking from SU(2); x SU(2)p to 
SU(2)r+r. In order to extend our analysis in Sec. 3.1 to spontaneous breaking of continuous 
symmetries, we developed a computational way of Lefschetz thimbles with approximate con- 
tinuous symmetries in Sec. 3.2. If there exists an exact continuous symmetry, the set of saddle 
points form a critical orbit of the Lie group, and there exists a remnant of that structure even 
after adding a small explicit symmetry breaking. In this case, we can approximately decouple 
the quasi-zero mode of the gradient flow to compute Lefschetz thimbles. 

By applying this technique, we study the spontaneous breaking of the axial U(1) symme- 
try of the 0-dimensional Nambu—Jona-Lasinio model in Sec. 3.3. After overcoming technical 
difficulties of computing Lefschetz thimbles with continuous symmetries, we can analyze this 
model in the almost same way as what we have done for the Gross-Neveu model. The only 
difference is that picking up one Lefschetz thimble has a singularity in taking the chiral limit, 
and thus we must sum up Lefschetz thimbles related by the gradient flow along the quasi-zero 
mode. 

It is an important future study to extend those analysis of 0-dimensional models to that 
of lattice field theories or other matrix models. When many degrees of freedom are correlated 
with one another in a complex manner, then the perturbative fluctuation around the saddle 
points becomes important. It is interesting to see whether our techniques and results in this 
chapter are useful under such circumstances. 
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Chapter 4 


Lefschetz-thimble approach to the 
Silver Blaze problem 


In this chapter, the Lefschetz-thimble path integral is applied to the sign problem of several 
models. In this circumstance, the mean-field approximation in a usual sense cannot be applied 
because the self-consistent equation has no real solution. When the classical action is complex, 
we need complexify integration variables in order to search relevant Lefschetz thimbles. Physical 
observables must have real expectation values; then how can this condition be satisfied? This 
question is answered in Sec. 4.1. 

In Sec. 4.2, the one-site Hubbard model is studied in detail, which is a good playground 
to study the sign problem in the path-integral expression. This model is easily and exactly 
solvable, but it has the severe sign problem in the path integral, whose mathematical structure 
is quite similar to that of finite-density QCD at uw, > m,./2. Using the Lefschetz-thimble 
method, we unveil the secret of the severe sign problem from various points of view. 


4.1 Sign problem in the mean-field approximation 


When the sign problem exists, there is no solution of the self-consistent equation of the mean- 
field approximation. When we try to solve it, the mean field usually takes complex values, which 
contradicts the real-valuedness of physical observables. This is the sign problem appearing in the 
mean-field approximation, and we consider its origin from the viewpoint of Lefschetz thimbles. 
We first show the general theorem to ensure the real-valuedness of physical observables in 
the Lefschetz-thimble decomposition, and applies that theorem to the complex saddle-point 
approximation. Using this method, the sign problem of the heavy-dense QCD is shown to be 
completely understood. The discussion in this section is adapted from Ref. [109]. 


4.1.1 Real-valuedness of the partition function and observables 


Let us consider a multiple integration that gives the partition function, 
Z= f d”z es, (4.1.1) 


where S(x) is the action functional of the real field x = (x1,...,2,). Since the partition function 
(4.1.1) for physical systems is given by the trace of a positive operator, 


Z=tr exp -B| ; (4.1.2) 
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it must be a positive quantity; however the Boltzmann weight exp (—S(a)) can be complex in 
general. In this section, we elucidate how the real-valuedness of the partition function holds 
manifestly in the Lefschetz-thimble decomposition 


Z= Y o e | d”z eS), (4.1.3) 


ced 


of (4.1.1). 
One of the sufficient conditions for ensuring Z € R is the existence of a reflection symmetry 
L: (xi) + (Lijx;), which satisfies Li; = Lj; € R, L? = 1 and 


S(x) = S(L- x). (4.1.4) 


In the case of finite-density QCD, we will see later that the charge conjugation gives this 
correspondence. The main purpose of this section is to derive the useful formula for the mean- 
field analysis and its systematic improvement while keeping Z real under Eq. (4.1.4). 

Let us recall Morse’s flow equation which defines Lefschetz thimbles, 


= = (S, (4.1.5) 


In order to show invariance of the Lefschetz-thimble decomposition under the complex conju- 
gation, some invariance of the flow equation (4.1.5) associated with the reflection symmetry 
(4.1.4) is sufficient. We extend the linear map L on R” to an antilinear map on C” by [109] 


The flow equation (4.1.5) shows the covariance under this antilinear reflection: Using Eq. (4.1.4), 
we find = 
dz OS(L -Z) 
= 4.1. 
= ( s ), (4.1.7) 


and thus the antilinearly transformed function z’(t) := K(z(t)) satisfies the same gradient flow 
equation (4.1.5): 


dete r (Pe) = Lig Lay: (32) 


= (3): (4.1.8) 


Many important properties can be deduced from this conclusion [109]: 


1. When 27 is a saddle point, so is K(z’). Let Jọ be the Lefschetz thimble around 27; 
then JE := {K(z) | z € Je} coincides with the Lefschetz thimble around K(z?) up to 
orientation. 


2. Since the flow is covariant under the antilinear map, (Ko, R") = (KE, R”) under an 
appropriate choice of the orientation. This means that the contributing thimbles always 
form an invariant pair J, U J, K : 
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Let us decompose © into three disjoint parts. For simplicity, we assume that S(z7) € R only if 
the saddle point satisfies z7 = K(z’) 1; then, © = No UUs U Ne, where 


Xo = {0o |2 =L-27}, Me = {0 | ImS(27) 2 0}. (4.1.9) 


The transformation K induces a one-to-one and onto correspondence between Uy and Ne. 
Equation (4.1.3) becomes 


Z= ` ne | d”z e75) 4 > a zee), (4.1.10) 


aeX0 TEexs t+ IF 


Each integral on the r.h.s. of the formula (4.1.10) is real or purely imaginary depending on 
whether K changes orientation of Js U Rs . Since the l.h.s. is real, n, must be zero unless the 
integral on J, + J* is real ?. This conclusion can also be applied to expectation values of any 
physical observables that satisfy the symmetry (4.1.4). The decomposition formula (4.1.10) 
takes a suitable form for the saddle-point analysis. 

Let us emphasize here that the real-valuedness of the partition function is ensured only by 
the invariance under the antilinear map K. The point of our discussion is that the Lefschetz- 
thimble decomposition of the integration cycle manifestly respects the antilinear reflection K, 
and so does the saddle-point analysis based on it. 


4.1.2 Mean-field approximation and complex classical actions 
Charge conjugation in finite-density QCD The QCD partition function at temperature 
T = 6"! and quark chemical potential uq is given by 

ZQocp = [PA det D(A, fg) eM 4), (4.1.11) 


where Sym = $tr y da* f d°x|F,.|? (> 0) is the Yang-Mills action, and 
det D(A, p4) = det |y” (3, + igA,) + fha + M] (4.1.12) 


is the quark determinant. When u4 # 0, the quark determinant becomes an oscillatory func- 
tional of the gauge field A, and the sign problem emerges. Even when uq # 0, the charge 
conjugation A œ> —A‘ with the ys hermiticity implies that the fermion determinant still satis- 
fies the identity [113, 114], 


det D(A, fq) = det D(A, — Hq) 
= det D(—A, u). (4.1.13) 


The charge C and complex K conjugation, or the CK transformation, serves as the antilinear 
map (4.1.6) for finite-density QCD [115, 116]. The Lefschetz-thimble decomposition (4.1.10) 


1This technical assumption can be removed by carefully defining X> so that they are conjugate to each other 
under the antilinearly extended map K(z) of L. 

?We here assume that the coupling constant is generic so that Stokes phenomenon does not occur. If it does, 
several saddle points are connected via some downward flows, and subdominant saddle points may give a purely 
imaginary contribution so as to cancel ambiguities of the large-order perturbation series at the dominant saddle. 
This is discussed in the context of the resurgence trans-series of quantum field theories [35, 37, 40, 41, 75, 110- 
12). 
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manifestly respects the CK symmetry so that Zacp € R. Since the discussion is robust, this 
conclusion applies to any effective descriptions of QCD, including lattice QCD simulations 
[67—70, 117] and also matrix models. 

We apply our insight on the Lefschetz-thimble method to the sign problem in Polyakov- 
loop effective models. The Polyakov loop £ of the fundamental representation 3 is an order 
parameter for the confinement-deconfinement transition [24]; 


B 
ly = Str P exp (of Aya") ; (4.1.14) 
0 


where P refers to the path ordering. To understand its properties, the constrained effective 
action Sog(a4) of the Polyakov loop is useful [118, 119]. It describes the (complex) weight of 
the partition function under the background Polyakov-loop phase a4: 


exp (—Seg(a4)) = [PA ell 8 (A, — a4), (4.1.15) 


where A; is the temporal gauge field and S[A] = Syy[A] — In det D(A, uq). We implicitly take 
the Polyakov gauge fixing with 


exp G) = diag [E O e e g e, (4.1.16) 


The Weyl group action (01,02) ++ (—01, 02), (1, 02) > ((01 + 302)/2, (01 — 02) /2) just permutes 
eigenvalues of the Polyakov loop (4.1.16), and thus the parameter region can be restricted to € = 
{0 = (01, 02) | 302| < 01 < a}. This background field method, or the mean-field approximation 
(MFA), gives a single integration over a4(01, 02) to compute the partition function: 


Zeacp = f anwo, 02) exp [- Seg (01, 02)], (4.1.17) 
€ 


with H(0) = sin? 6; sin?((01 + 302)/2) sin?((01 — 362)/2) coming from the SU (3) Haar measure 
[120, 121]. When the quark chemical potential qx is turned on under the nontrivial Polyakov- 
loop background, the effective action S.@(@) takes complex values in general due to the quark 
determinant. This makes the integration (4.1.17) oscillatory, and the sign problem remains in 
the mean-field approximation [113, 114]. 


Lefschetz-thimble analysis of heavy-dense Polyakov-loop models In order to evade 
the sign problem of the MFA, we should rewrite Eq.(4.1.17) based on the decomposition formula 
(4.1.10). Since H(@) identically vanishes at the edge of €, the Lefschetz-thimble method can 
be applied by regarding S(0) = Ses(0) — In H(@). In terms of the Polyakov-loop (4.1.16), the 
CK transformation K becomes 

K(a, 22) = (zZ, —%), (4.1.18) 


where z; is the complexified variable of 6;, up to some Weyl transformations. The CK symmetry 
(4.1.13) leads 
Ser (21, 22) = Sem (Z1, —22). (4.1.19) 


It is important to remark that the Polyakov loop, 


l3(0) = ; (2e' cos 04 + 67%”), (4.1.20) 
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Figure 4.1: Flows at h = 0.1 and u = 2 around the CK-symmetric saddle point (the black blob) 
z* in the Re(z,)-Im(z2) plane. The red solid and green dashed lines are the Lefschetz thimble 
J, and its dual K, of z*, respectively. Since X, intersects with Im(z2) = 0, the integral over 
J, is equal to that over €. (The figure is taken from Ref. [109]) 


satisfies the CK symmetry (4.1.19). For our demonstration, we take a simplified heavy-quark 
model [113, 122-127]: 


(3° 


— 1 
Sef = T (ets(0, 02) + e ”lz(0, 0») 


8h 
er (2 cos 0; cos(@2 — iu) + cos(202 + iu) ). (4.1.21) 


Since the Polyakov loop £3 can be regarded as a free energy for single quark excitation, this 
effective action represents the free energy for static quarks. Therefore, this approximation can 
be justified when quarks are sufficiently heavy. When h 4 0 and u = Hak # 0, the integration 
(4.1.17) is oscillatory because Seg takes complex values. The gluon dynamics is neglected just 
for simplicity. 

If h = 0, the system has the Z symmetry generated by (61, 02) > (01, 02+ 27/3) up to Weyl 
group actions. The Haar measure factor H(0,,62) takes the unique maxima at 6* = (27/3,0) 
in €, and this is Z invariant. For later purposes, it is important to notice that this saddle 
point is CK symmetric in the sense that 6* = K(6*). The eigenvalues of the temporal Wilson 
line (4.1.16) become uniformly separated, and thus (43) = 0. 


When h Æ 0, the Z symmetry is explicitly broken, and the saddle point z* moves away 
from the Z3-symmetric point (27/3,0), but it still must be CK symmetric because there are no 
other saddles to form a pair. By continuity of intersection numbers, z* will contribute to Zacp 
even for finite h and u. In order to demonstrate it lucidly, we explicitly solve the flow equation 
(4.1.5) around the CK-symmetric saddle point (see Fig. 4.1). The dual thimble K, of z* is 
shown with the green dashed curve, and it indeed intersects with the original integration cycle 
€. Therefore, the complex saddle point contributes, but physical quantities take real values 
since the CK symmetry is respected under the Lefschetz-thimble decomposition. 


In the limit u —> +00, the saddle-point approximation becomes accurate. In this limit, the 
effect of the quark potential Sep becomes quite large. By solving the saddle-point equation of 
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S in the limit u > oo, we find that 


3e7H/2 e # 
~ o a 4.1.22 
1 Wh 2 Sh ( ) 
and it indeed approaches the perturbative vacuum 2%... = (0,0) of Seg. Thus, (£3) and (£3) 


converge to 1 in this limit, as was numerically observed in Fig.2 of Ref. [113]. Using the 
saddle-point approximation, we can see that these Polyakov loops have different expectation 
values: 


2 
(fz) — (€3) >~ 3 (sinh 2iz} — 2 cos z{ sinhiz>) > 0. (4.1.23) 


These results are consistent with those of exact computations of expectation values [113], and 
the CK-symmetric saddle point gives real expectation values [115, 116]. The formula (4.1.10) 
shows that this approximation can be systematically improved perturbatively to satisfy these 
physical requirements. 

In this model, the flow equation tells us that only the single saddle point contributes to the 
path integral. In order to make our mean-field analysis valid, the one-thimble approximation 
must be sufficiently accurate and this condition can be checked by solving the flow equation in 
principle. 


4.2 One-site Hubbard model 


In this section, we apply the Lefschetz-thimble approach to the one-site model of electron 
systems. This toy model can be regarded as an extreme limit of strong couplings, because it 
can be obtained by neglecting hopping in the Hubbard model, and its Hamiltonian can be easily 
diagonalized. If we rewrite its partition function using the path integral with the Hubbard— 
Stratonovich transformation, the sign problem emerges and thus this toy model provides us 
a good playground to study theoretical structures of the Lefschetz-thimble approach when we 
apply it to the sign problem. In a previous study [117], the two-dimensional Hubbard model 
is studied using the Lefschetz-thimble Monte Carlo method, but one must use the so-called 
“one-thimble approximation” due to the current limitation of the numerical algorithm. Since 
our system is exactly solvable, we can understand the complete structure of Lefschetz thimbles, 
and suggest appropriate approximation schemes. The discussion in this section is adapted from 
Ref. [128]. 


4.2.1 Physical properties of the model 


The Hubbard model [129-131] is a lattice model of many-body electron systems, and it can 
describe the transition between Mott insulators and metals. Its dynamics is governed by 


A =t X Xis +U > Arityi— ud (fra + yy). (4.2.1) 


(ij) 


The summation is taken over all the lattice sites 7, and the notation (i, j} denotes that only the 
nearest neighbor hopping is considered. él and €,,, are creation and annihilation operators of 


A 


fermions of the spin o(=f, |) at the site i obeying {ei jj 


d = 057043, and fig, = om is the 
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number operator. The parameter t is called the hopping parameter, and we will consider the 
case t = 0 in this section. The parameter U(> 0) describes the on-site repulsive interaction, and 
the energy increases if more than one particles get together at the same site. The Hamiltonian is 
invariant under the global U(1) rotation alt) = exp(+tia)él), and the total number of fermions is 
conserved. In the Hamiltonian (4.2.1), the chemical potential u is introduced for this conserved 
quantity. 

At zero temperature, the Hubbard model shows a Mott insulating phase ift < U. Especially 
in the strong coupling limit t = 0, one can solve the model (4.2.1) analytically and let us see 
the result on the total number density (7) as a function of the chemical potential u. Since each 
site is totally independent from others, one can study the single-site Hamiltonian, 


H = Uô, — ph, (4.2.2) 


instead. Since the Hamiltonian (4.2.2) commutes with the number density operator ù = Ay++ny, 
we can take the number basis to find the ground state. Let us defines the partition function 


Z=tr [exp (-2å)| : (4.2.3) 
then the number basis gives an explicit result: 
Z = 1 + 2e®” te UH) (4.2.4) 


The number density is given as 
10 2(eßb —(U—2u) 
nZ = (eP. +e ) . 
B Op 1 + 2e%# + e-8U—2n) 
and, in the zero-temperature limit 6 — oo, it shows a non-analytic behavior: 
0 for pO <0, 


n(8 =œ)=4 1 for 0< p/U <1, (4.2.6) 
2 for pu/U>l1. 


At the half-filling u = U/2, Eq. (4.2.5) shows that (nr) = 1 for any 8, which comes from the 
particle-hole symmetry iH 2 — ñ (or, Ao > 1 — ño) in (4.2.2). 


(4.2.5) 


4.2.2 Path integral formulation and the sign problem 


We have explicitly seen that the Hubbard model can be solved analytically using the number 
eigenstates in the strong coupling limit, or the static limit t = 0. However, if one changes the 
basis of the Hilbert space in taking trace, the sign problem emerges and the Hubbard model in 
the strong coupling provides us a good lesson. 

Let us first derive the path integral expression of the partition function. Using two- 
component complex Grassmannian variables y = (py, Y), the fermion coherent state is defined 
by 

I) = exp(—vel)|0), (4.2.7) 
which satisfies ¢|w) = y|). In the following, the summation over spin is taken without showing 
it explicitly. We take the convention in which all the Grassmannian variables anticommute with 
each other and with fermionic creation/annihilation operators. It is easy to check that 


(hila) = et, f dy*dyjyje "(y| = 1, (4.2.8) 
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and 


tr(O) = f dy*dye™®" * (—y Oly). (4.2.9) 
Let |n+,n,) be a fermionic Fock state, and the inner product with coherent states is 
(0,01%) = 1, (1, 0/0) = Yr, (0; 11%) = dy, (1, IY) = —d yyy. (4.2.10) 


Using these relation, one can compute the matrix element of the time-development operator in 
the following way: 


e PP ple A77 leh!) = exp [y (pb — eA") + (1 — €O7Y) (hte*™ Ha) (ype )] 
~ exp — |y* (y — eA!) + ATU (žep) (piee )] : (4.2.11) 
The first equality is exact thanks to the knowledge on the number eigenstates, and we take an 


approximation for small A7 in the last equality. In order to circumvent the quartic interaction, 
let us introduce an auxiliary field y vis the Hubbard-Stratonovich transformation, so that 


AT AT ; wishes ade j 
V xU [evo (y — iU (wper™*yy + vend). (4.2.12) 


Inserting (4.2.12) into (4.2.11), we obtain that 
edt 
AT p + 3 wr ( yen iare): 


a FY (abe) SN oa a7 fwe xp — 
(4.2.13) 


In order to take the naive continuum limit of the path integral expression, we must exponentiate 
the auxiliary field in the fermionic bilinear term: 


Ar? 
e^™H(1 +iATy) = exp (ar tip) 4 T 4 oars!) 


. Ar a 3/2 
exp | AT(u + iy) + FV Jo + O(AT” $) 


exp (a7 (u + £ + ie) + olar) ' (4.2.14) 


Here, ((y”)9 means the expectation value of y? with the Gaussian weight in (4.2.13), which gives 
U/Ar. We find the path integral expression of the partition function as 


k=1 k=1 
x f Tlia exp (-> > tia (vest -eN enna) (4.2.15) 
k=1 k=1 
with the antiperiodic boundary condition Wy, = —v,. Equation of motion in terms of y shows 


(â) = (Wire oh erritUPy,) = — 2 (oy) (4.2.16) 
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The expectation value of y is nothing but that of the total number density. This formula 
will be used later in order to compute the number density. Since Eq. (4.2.15) is quadratic in 
the fermionic field Yy and w*, we can perform the Grassmannian integration explicitly. The 
spectrum of the fermion bilinear operator is given by 


N 


—1)ri B : 
Aeg, u) = 1 — cP exp N? S "(ive +p+U/2). (4.2.17) 
k=1 


and then in the continuum limit 


sim. Mi x) l = (1 + exp [ dr (ip(T) + u + u/2)) (4.2.18) 


Since fermions have two flavors, the overall square is taken. Since the fermionic determinant 
contains only the Matsubara zero mode of y, the path integral of non-zero Matsubara modes 
of y gives a trivial Gaussian integration without coupling to u. Therefore, the path integral is 
now reduced to an integral of zero Matsubara mode Ypg, and we obtain that 


| B U B 
= sag | Pr (1 + exp p (ie +u+ a exp -pp Pte: (4.2.19) 


This integration can be performed analytically to find (4.2.4). Instead, we apply the Lefschetz- 
thimble method because (4.2.19) has the sign problem, and this simple model may lead us to 
the deeper understanding of the problem. 


Sign problem 


In order to explicitly show that the path integral (4.2.19) contains the sign problem, let us 
compare it with the sign-quenched partition function. In order to make the sign problem less 
severe, we first shift the integration variable yp, to 


Pog — 210 for u/U >} 
b=) Ye—i(ut 5) for =3<u/U <3, ea 20) 
Pog for p/U <—3. 


After shifting the integration variable, we take absolute values of the integrand, and the result 
is called the phase quenched partition function. For u/U < —1/2, 


ing = | =O flid+n+U/2) |’ e7882. (4.2.21) 
Jar a 
For —1/2 < u/U < 3/2, 
z= de [1 + e186)? @-80?/2U+B(u+U/2)?/20 (4.2.22) 
Pa \/2nU/B 


The ratio Z/Zpq is shown in Fig. 4.2 as a function of the chemical potential at GU = 30. For 
small (large) chemical potentials, u/U < —1/2 (u/U = 3/2), the ratio Z/Z,, is almost one 
because the dominant contribution comes from the unoccupied (occupied) state, and the sign 
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Z/ ‘Zquench 


Figure 4.2: Z/Z,, as a function of u at BU = 30. 


problem is sufficiently weak. In the later section, we can interpret this semiclassically because 
the important configurations accumulate around the origin Ypg ~ 0 (Ypg ~ 2iU). According to 
Fig. 4.2, there is another remarkable condition in which the sign problem disappears again by a 
simple shift of the integration variable: At the half-filling u = U/2, the shift of the integration 
variable solves the sign problem. By introducing a new integration variable by Ybg = Ybhg + WU, 
we find that 


| B bo B 
b 2 
Zhalf-filling = al Yoo (2 cos —; | exp - Fe -—U’)|, (4.2.23) 


and the integrand is positive definite. The sign problem can be evaded in these three special 
case even using the path integral formalism, but this does not happen in general. Indeed,when 
u/U = 0 or 1 at which number densities (4.2.6) jumps, the sign problem is severe because 


3 
Bg oo <1. (4.2.24) 


We explore the origin of this sign problem in the path integral formalism by applying the Picard— 
Lefschetz theory to the integration (4.2.19). If we naively use (4.2.21) instead of (4.2.22) for 
the whole region, the sign problem becomes more severe so that Z/Zpą x e-° for u = 0. 

When the sign problem becomes severe u/U = —0.5, the second-order phase transition 
happens in the phase quenched theory at T = 0. This situation is quite similar to what 
happens in the finite-density QCD at ug = m,/2. 


4.2.3 Lefschetz-thimble analysis at strong couplings 


Let us apply the Picard—Lefschetz theory to (4.2.19). The effective action of this system is 
obtained as 
B 


S(z) = Tis — 2In (1 + exp B (i: ++ 3) ; (4.2.25) 


which satisfies the real-valued condition of the partition function (4.1.4); S(z) = S(—Z). As 
we have shown in Sec. 4.1, this gives a strong constraint on the configuration of Lefschetz 
thimbles to ensure the real-valuedness of physical observables in the decomposition (4.1.10). 
The logarithmic function has branch singularities at fermionic Matsubara modes z = i(u + 
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Figure 4.3: Behaviors of Morse’s downward flow equation at GU = 30, for u/U = 2 (a) or 
u/U = 0 (b). Star-shape black points show singular points of logarithm, and red blobs show 
complex saddle points z,. We only show the region Re(z) > 0 because of the symmetry of the 
flow equation under z ++ —z. (Figures are taken from Ref. [128].) 


U/2) + (2l + 1)r/8 for £ € Z. The integrand converges to 0, i.e., S — +00, at these points. 
The flow equation of Eq. (4.2.25) reads 


dz Zz  2iexp6(-iz+p+ 2) 


- . (4.2.26) 
dt U 1+ exp 8 (-iz+ y+ 9) 


Semiclassical analysis and gradient flows Let us find the set of saddle points, which is an 
important step not only for the saddle-point approximation but also for the Lefschetz-thimble 
method. The saddle-point condition of the effective action (4.2.25) is 


2U 


— - f (4.2.27) 
1 + exp —8 (ize + u+ E) 


2y = 


In order to simplify our analysis, we consider a limiting case where T « U,|y|. We can 
approximately obtain the saddle points as 


3U — 
lm =i (u + z) +T (20m +iln3, T £) +O(T°) (4.2.28) 
oth 


2 2 


for m € Z by assuming that the second term is much smaller than the first one. The real-valued 
condition says that z_, and zm form a pair. If u > 3U/2 or u < —U/2, there is another solution 


_ f %U+o0(T) for p> 3U, 
i { 0+o(T) for u< —4U. eee) 


Behaviors of the downward flow equation (4.2.26) are shown in Fig. 4.3 with U = 1 and 
BU = 30. Let us discuss the case when the chemical potential is sufficiently large (u = 2U). 
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u/U 


-1/2 0 1 3/2 


Figure 4.4: Schematic behavior of the number density. The exact result at 6 = oo is shown 
with the solid blue curve. The dashed red line shows the naive expectation from the positions 
of the complex saddle points, Im(z,,) ~ u + U/2, for —1/2 < u/U < 3/2. 


In this case, Morse’s flow is shown in Fig. 4.3 (a). Only one Lefschetz thimble J, associated 
with z, contributes in the Lefschetz-thimble decomposition. Other dual thimbles shown with 
green dashed lines do not intersect with the original integration cycle R as in Fig. 4.3 (a). Thus 
intersection numbers in Eq. (4.1.3) vanish and we have 


z= | dzexp —S(z). (4.2.30) 


Within the saddle-point approximation using the complex saddle point z,, the number density 
is given as 


— iZ, 9 

U | 
and thus the saturation of fermions is well described. The exactly similar thing happens also 
for u/U < —1, and then the saddle-point approximation gives n ~ —iz,/U = 0. This analysis 
explicitly shows why the sign problem is significantly weakened with an appropriate shift of 
integration variables for u/U < —1 and u/U = 2, as we have show in Fig. 4.2. 

In the following, let us concentrate on the case where the sign problem is severe, —1/2 < 
u/U < 3/2, at low temperatures BU = 10. In this parameter region, the typical behavior of the 
flow is shown in Fig. 4.3 (b). All dual thimbles Km intersect with the original integration cycle 
R. This shows that all the saddle points zm contribute to the partition function®, so that the 
interference among them may not be negligible. This interference requires a careful treatment 
of the semiclassical analysis in order to solve the sign problem. Behaviors of the flow equation 
in Figs. 4.3 (a) and (b) are very different. This big difference can be explained by the Stokes 
phenomena [64], which occur around u/U œ 3/2 and also around u/U œ~ —1/2. However, 
this is irrelevant to the non-analytic behaviors of the number density, and then we stick to the 
analysis for —1/2 S u/U S< 3/2. 

There is a big question on the behavior of the number density for —1/2 < u/U < 3/2. Since 
the number density is given by n = Im(z)/U, we may naively expect that the number density 


ne (4.2.31) 


3This statement can be shown at half-filling due to the particle-hole symmetry, but in general cases we have 
checked it numerically only for |m| < BU/2r. For sufficiently large |m|(s> GU/27), the Stokes phenomena may 
happen and it may cause a peculiar behavior. However, the effect of that subtlety is quite small, and it does 
not change the result of the following argument. 
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is 

ali( 2) ~ $ a (4.2.32) 
by using (4.2.28). This is shown with the dashed red line in Fig. 4.4, which is totally different 
from the exact result (the solid blue curve). This behavior is, however, also expected from the 
phase quenched theory Z,, in (4.2.22). The similar problem happens in the finite-density QCD 
at Hq = m,,/2 (see Sec. 1.2), and thus this model is a prototype of the sign problem appearing 
in the finite-density QCD. This is the Silver Blaze problem in the one-site Hubbard model". 
This aspect will be discussed in detail in Sec. 4.3. 


Semiclassical partition function Let us denote the classical action at the saddle point z, 
as S,(:= S(z,)). Substituting the approximate expression (4.2.28) into Eq. (4.2.25), we have 


Qn? 


Re (Sm — So) > au” (4.2.33) 

Im Sm ~ 27m (5 + 5) ; (4.2.34) 
and 

So > "a (r + =) l (4.2.35) 


Equation (4.2.33) shows that subdominant thimbles Jm can be comparable with the dominant 
one Jo for BU > 1 so long as |m|(# 0) is not so large. According to Eq. (4.2.34), these 
different thimbles have different complex phases, and thus a contribution from one Lefschetz 
thimble is in general canceled by other ones. Exactly at the half-filling, u/U = 1/2, the 
complex phase is always an integer multiple of 27i, and such a cancellation is absent. This 
gives an interpretation on the reason why the sign problem disappears at the half-filling from 
the viewpoint of Lefschetz-thimble approach. 

Let us compute the partition function Za only using these semiclassical information. Ne- 
glecting the higher order correction in the Gaussian fluctuation, let us simply compute 


Za := ` exp — Sn. (4.2.36) 
This can be computed analytically in the following way using the elliptic theta function: 


= 1 
Zaz e 50 (: +2 P cos 2NT ($ + 5) pa 


n=l 


1 
= e9; (- (5 + 5) ise ) (4.2.37) 


“Recently, the one-site massive Thirring model is studied in Refs. [132-134] using the Monte Carlo simulation 
with one-thimble approximation. The one-thimble approximation reproduces the naive expectation, which does 
not explain the Silver Blaze problem correctly. 
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Here we used an approximate expression of the saddle points (4.2.33-4.2.35) to obtain (4.2.37) °. 
Using this result, in the limit 6 — oo, one can analytically show that 


1d 2 (1< u/U < 15), 
Na = 39, Poa >< 1 (0<pu/U <1), (4.2.38) 
Bou 0 (—0.5 < u/U <0). 
Of course, this good agreement is accidental to some extent, since we have neglected higher order 


contributions in the semiclassical analysis. This result, however, still indicates the usefulness 
of the semiclassical approximation even when the sign problem is quite severe. 


Lee—Yang zeros and fermion spectrum Let us discuss this non-analytic behavior from the 
viewpoint of Lee-Yang zeros [97, 98, 100-107] in this semiclassical expression (4.2.37) (For Lee- 
Yang zeros see Sec. 3.1.3). Using the infinite-product expression of the elliptic theta function, 


oO 


63(z,q) = I (1 — g^) (1 +287 cos(2z) + 4*7?) k (4.2.39) 
=1 
the zeros of the partition function in the complex u plane satisfy 
(2 a ) sal wn (4.2.40) 
cos | 2r= ) = — (2 — 2. 
U BU 


for some positive integer £. Since (4.2.37) is originally derived as an approximate expression 
for —1/2 < u/U < 3/2, we restrict the real part of the complex chemical potential to —U/2 < 
Reu < 3U/2. Then, the zeros of the partition function are 


u =inT (2k — 1), u =U +inT(2k — 1), (4.2.41) 


for k € Z. In the limit T — 0, these zeros form straight lines crossing to the real axis at 
u = 0,U, which causes the first-order phase transition as shown in (4.2.38). 

Let us apply the same analysis to the quenched partition function. Let us again define the 
phase quenched semiclassical partition function by 


Zp i= ` lexp — Spn]. (4.2.42) 


n=— o0 


In the phase quenched approximation, interference of complex phases expilmS,, disappears, 
and the chemical potential dependence on p is given by exp — So. This is consistent with (4.2.22). 
Lee-Yang zeros does not exists if —U/2 « Reu < 3U/2. Therefore, the only possibility of the 
phase transition in this phase quenched approximation originates from the Stokes phenomena 
around u/U ~ —1/2,3/2, which usually gives the second-order transition [91]. 

In order to investigate this failure of the phase quench approximation more deeply, we 
analyze the fermion spectrum (4.2.17) at complex saddle points. For simplicity, let us take the 
continuum limit at first with an appropriate renormalization, 


U 
Ap (Oh) = dim TNA, u) = — fe —1)niT +ip +u + 5} . (4.2.43) 


>The Gaussian integral gives an m-independent prefactor (1 + BU (u/U — 3/2) (u/U + 1/2) 71 at the lead- 
ing order. It does not change a qualitative behavior and thus it is neglected in the following. 
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The spectral representation of the number density is given by 


i0+e 2 


2 € 
n(y,u) =—2T > 


ng, u) e BG HU?) + (4.2.44) 


l=- 
If the sign problem is mild enough, u < —U/2, then the eigenvalue at the saddle point z, = 0 
is 


A(z, pw) & — (oe —1)riT +t >) (4.2.45) 


Since ReAj°"(z,, u) cannot be zero, the chemical potential dependence of the number density is 
exponentially suppressed as T + 0, 


2 


However, for u > —U/2, Reù can be zero in the vicinity of saddle points, and indeed one finds 


U 
2 


3U — u 
A (m u) = — | (2(€ +m) — 1riT -Tln 2 zs (4.2.47) 
u 


Within this approximation, the number density becomes 


2 H 1 
a ft) e~blzm+u+U/2) 4 1 = U 27 eee) 


The near-zero fermionic mode provokes the fake second-order transition at u = —U/2 again, 
as shown with the dashed red line in Fig. 4.4. From the viewpoint of fermion spectrum at 
complex saddle points, nothing special happens in the vicinity of the correct transition points 
u = 0,U. In the language of Lefschetz thimbles, this implies that their topological structure 
does not change at all around u = 0,U, as we can see in Fig. 4.3 (b). 


4.2.4 Importance of interference 


Numerical results Let us consider one-, three-, and five-thimble approximations in order to 
analyze the importance of the interference among multiple Lefschetz thimbles. We consider the 
case where the sign problem is severe, i.e., —1/2 S u/U X 3/2 and BU = 10. The (2m + 1)- 
thimble approximation takes into account only Jo, J+1,.--,; Jem in Fig. 4.3 (b). The partition 
function (4.2.19) in this approximation reads 


Z~Z™ := I dz e750) 4 Snare f dz e7503, (4.2.49) 
Jo k=1 Tk 


Each contribution of the summation in the second term comes from a pair Jy and Jp, and 
becomes manifestly real because of the reflection symmetry z + —Z (see Sec. 4.1 and Ref. [109]). 
Since the Lefschetz thimble Jp is homologically equivalent to the interval in the complex plane 
i(u + U/2) + [(2k — 1)nT, (2k + 1)nT], the above integration (4.2.49) is reduced to 


+(2m+1)rT 
Zim) — f dz exp -S(x +i(u + u/2)), (4.2.50) 


(2m41)rT 


72 Chapter 4. Lefschetz-thimble approach to the Silver Blaze problem 


TL 
Z0F 4 i > 
--- 3-thimble | +: \ 
— 5-thimble A i 
rsh! } 


R E kT 


mE a 


= 
| 
= 


0.5 1.0 1.5 


Figure 4.5: Behaviors of the number density n = —i(y)/U as a function of u at BU = 30. The 
thin solid black line shows the exact solution. Other lines, dotted green one, dashed blue one, 
and thick solid red one, show the result of one-, three-, and five-thimble approximations, which 
integrate over Jo, Jo U Jai, and JU S41 U J2, respectively. The five-thimble approximation 


gives an almost exact result, and most parts of the solid red line overlaps with the solid black 
one. (This figure is adapted from Ref. [128].) 


with x = z —i(u + U/2). The result on the number density at GU = 30 is shown in Fig. 4.5, in 
which one-, three-, and five-thimble approximations are compared with the exact computation. 
The one-thimble approximation, which is shown with the dashed green line, is an approx- 


imation to integrate only over the Lefschetz thimble J, and it almost gives the mean-field 
result, 


H 1 
NMF = U + z’ (4.2.51) 
for —1/2 < u/U < 3/2. Figure 4.5 shows that the one-thimble approximation is not sufficient 
to describe the plateaus of the number density in the region —1/2 < u/U < 3/2. One can 
easily check that this is also obtained from the sign-quench approximation (4.2.22). 

The result of the three-thimble approximation provides us a useful lesson for an application 
of the Lefschetz-thimble approach to the sign problem. The number density diverges at several 
chemical potentials around a rapid crossover, although the result is improved around each 
plateau. Let us analyze this divergence by using the semiclassical analysis. For that purpose, 
we introduce the semiclassical partition function with the (2m + 1)-thimble approximation by 


Zim) = > exp — Sk. (4.2.52) 


k=—m 


In the three-thimble approximation, the semiclassical partition function approximately behaves 


as 
E 2 
Z® ~ eS h +207 30 cos m (1 $ a) (4.2.53) 
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Since GU > 1 and thus e~2*°/8U ~ 1, the second term is not suppressed, compared with the 
first term. This eer partition function vanishes at some chemical potentials around 
p/U =0 and 1 ( vanishes at u/U ~ +0.04, which is roughly consistent with Fig. 4.5). This 
example ane em that picking up a part of the Lefschetz-thimble decomposition 
can violate the physical requirement, such as the positivity condition on the thermodynamic 
quantities. In this instance, the incompressibility (On/Oy)/n? should not be negative. 
According to Fig. 4.5, we can conclude that five Lefschetz thimbles Jo U J+ı U S42 are 
necessary in order to explain rapid jumps of the number density in terms of the chemical 
potential at a low temperature GU = 30. Indeed, the five-thimble approximation (solid red 
line) and the exact computation (solid black line) almost overlap with one another in Fig. 4.5. 


Necessary numbers of Lefschetz thimbles How many Lefschetz thimbles are necessary 
in general at a given lower temperature 8? Let us propose the criterion to neglect the Lefschetz 
thimbles Jm for large |m]: 

2 _ gil) 


cl 


<1. (4.2.54) 


Assuming that the number of Lefschetz thimbles, (2m + 1), is sufficiently large, one can replace 
Z\™ in the denominator by the semiclassical partition function Za. Since jg hag" | < 
2| exp — Sml, the contribution from Jm can be negligible if 


2| exp — Sml 2 


4.2.55 
p—s (4.2.55) 


Here € is a controlling parameter of an error, and € < 1. Using the approximate results (4.2.33) 
and (4.2.37), we can solve this inequality with respect to m: 


lm) > 4/ -Z n £63 (n(u/U + 1/2), e-27?/8U), (4.2.56) 


This criterion gives different results depending on u, and thus let us first derive the strongest 
restriction. If the sign problem is severe, i.e. u = 0, the elliptic theta function is exponentially 


small with respect to BU: 
03 (5e) 4) 2BU pu /s. (4.2.57) 
T 


Therefore, for reasonable ¢ such as € œ 0.1, the effect of € is negligible because of its mild 
logarithmic dependence. The criterion gives 


in the limit of GU > 1. It means that we need (2[8U/4r| + 1) thimbles in order to describe 
the rapid crossover of the number density for the one-site Hubbard model. 

Let us discuss this behavior from alternative point of view. Since GU becomes large, the 
relevant integration region becomes smaller as |Re z| < ,/U/@. However, Fig. 4.3 shows that 
the length of each Lefschetz thimble J, is proportional to 1/8, and thus we need more Lefschetz 
thimbles in order to cover the relevant integration region. That number is clearly proportional 


(4.2.58) 
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Figure 4.6: Behaviors of the number density n = —i(y)/U as a function of u at BU = 100. 
In order to explain the Silver Blaze phenomenon, the 7-thimble approximation is already rea- 
sonable. In order to describe the correct transition behavior, 17 Lefshetz thimbles become 
necessary. (This figure is adapted from Ref. [128].) 


to y BU, and is much smaller than Eq. (4.2.58). This is because interference of complex phases 
among different saddle points makes the semiclassical partition function much smaller than 
that of the naive expectation if u/U ~ 0 and 1. This indicates that the explanation of the 
plateau in the vicinity of the half-filling is much easier than that of the rapid crossover. The 
necessary number of Lefschetz thimbles for each phenomenon is also largely different. Indeed, 
for u/U = 1/2, one can show that 


A oe BU 
AS, Os (r,e 2 ie) ~ E (4.2.59) 


and, roughly speaking, we need 


Be In v8 
27? e./BU 


Im| z 


(4.2.60) 


for Jm being negligible in order to achieve a good accuracy only around the half filling. This is 
consistent with the result of the above heuristic argument. Now, we can understand that the 
large gap of these two estimates (4.2.58) and (4.2.60) comes from the sign problem in summing 
up Lefschetz thimbles, and the semiclassical analysis provides a reasonable rough indication. 
In Fig. 4.6, we plot the number density at GU = 100, and the exact result is compared with 
those of 7- and 17-thimble approximations. The result of the 17-thimble approximation (thick 
solid red line) completely overlaps with the exact one (thin solid black line). The 7-thimble 
approximation (dashed blue line) is insufficient to describe the phase transition at u/U ~ 0 
and 1, but it nicely describes the silence of the number density in the vicinity of u/U ~ —1/2, 
1/2, and 3/2. 
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Figure 4.7: Schematic illustration of the Silver Blaze phenomenon in finite-density QCD. This 
plots the baryon number density ng as a function of the quark chemical potential uq. In the 
phase quenched QCD, the second-order phase transition happens at Hq ~ m,/2 (the dashed red 
line), but the full QCD must experience nothing until the nuclear liquid-vapor phase transition 
at uq ~ my /3 (the solid blue line). It is similar to Fig. 4.4 for the one-site Hubbard model. 


4.3 Speculation on the Silver Blaze problem at finite- 
density QCD 


We have studied the one-site Hubbard model in Sec. 4.2 as a toy model of the sign problem. We 
elucidate that the interference of complex phases among multiple saddle points is important 
to explain n = 0 for —U/2 < u < 0 and the rapid jump of n at u = 0. In this section, we 
first discuss the similarity of the mathematical structure of the one-site Hubbard model to that 
of finite-density QCD, and mention the speculation on the baryon Silver Blaze problem. The 
discussion in this section is adapted from Ref. [128]. 

Let us recall the baryon Silver Blaze problem. We consider properties of finite-density QCD 
at T = 0. In general, the quark determinant, Det [y4( 4 + m) — u4], will have the (nontrivial) 
/lg-dependence. Naively thinking, this means that the baryon number density ng arises for any 
chemical potentials. However, since we know that the lightest baryon in the QCD spectrum is 
nucleon, the baryon number density ng must be zero for any 3ug < my — B ~ 923 MeV at 
T = 0. It is difficult to explain this empirically trivial fact by using the path-integral expression 
of finite-density QCD. It contains the essence of the sign problem, and this is called the baryon 
Silver Blaze problem [58]. Indeed, the unphysical early onset of ng has been observed in the 
reweighting method of lattice QCD [135, 136] (Fig. 4.7 schematically illustrates this situation). 

So far, we have understood the baryon Silver Blaze problem only for small quark chemical 
potentials wu, <_m,/2 [58, 106, 137]. In general, the eigenvalues of y4( +m) can be expressed 
as 

Agm) Ay Hg) = €)(A) — Ha — 6 )(A) + (Qn + aT (4.3.1) 
with -nT < j < mT, and they do not depend on uq. Using this expression, the quark 
determinant becomes 


Det AU +m) = Haq] 


= Ne Pe ™O)/2 TT (1 4 PGs) ett) , (4.3.2) 
j 
= N I] ef les +i¢;) (1 + e Ales-Ha—i#)) (1 + e@ Ales tHe ties) . (4.3.3) 
E€j>0 


where N is the normalization factor. Here, we implicitly use the zeta-function regularization, 
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and see Ref. [137] for its derivation. Its j4,-dependence at sufficiently low temperatures and at 
[lq > 0 becomes crystal clear according to the following expression [58, 137]: 


Det [aa +m) = A Ka [I (1 ae g B(es(A)— Ha) giB5(A)) (4.3.4) 
Depar yell 

The quark determinant becomes independent of u4 at T = 0 if uq is smaller than the minimum 
of ¢;(A)(> 0). For statistically significant gauge fields, it is confirmed that min(e;(A)) = m,/2 
at Ug = 0 [58, 106, 137], and thus ng = 0 for uq < m,/2. This also means that the Lefschetz- 
thimble decomposition of lattice QCD at uq < m,/2 and at sufficiently low temperatures is 
identical to the original integration cycle up to an exponentially small correction because the 
quark determinant at 4g = 0 is positive-definite. 

In contrast, as computed in Eq. (4.2.18), the fermion determinant in the one-site fermion 
model is 

Det fo- = (vo) FEF al = (1 + ePUAP) (4.3.5) 
We can ask the same question: How can we show n = 0 for u < 0 using the path-integral 
expression when the fermion determinant depends on u? This Silver Blaze problem for u < 
—U/2 can be solved in the same manner for the baryon Silver Blaze problem at [ugk < m,/2. 
At u < —U/2, the fermion determinant becomes independent of u as T — 0. This also explains 
why the Lefschetz thimble J, at ps < —U/2 is almost identical to the original integration cycle 
R (The similar discussion shows the Lefschetz thimble J, at u > 3U/2 is almost identical to 
the line with Im(z) = 2). 

For —U/2 < u < 0, the -dependence of the fermion determinant (4.3.5) becomes exponen- 
tially large, and thus the property n = 0 at T = 0 is still veiled. At u ~ —U/2, zeros of the 
fermion determinant (star-shape black points in Fig. 4.3) move close to the real axis, and let us 
recall that Ją = R for u < —U/2 up to an exponentially small correction. In the case of one-site 
fermion model, this triggers Stokes jumps and the original integration cycle is decomposed into 
multiple Lefschetz thimbles for u > —U/2. One can no longer replace the expectation value 
of the number density by the number density at a complex saddle point. The significance of 
interference among multiple Lefschetz thimbles is identified in order to explain not only the 
rapid jumps of n but also the property n = 0 for —U/2 < u < 0 at the zero temperature. 
Furthermore, if the number of Lefschetz thimbles is insufficient, the thermodynamic stability 
is violated in the one-site fermion model. 

We speculate that this interference of complex phases also play an important role in finite- 
density QCD to analyze the Silver Blaze problem and also the sign problem beyond half of 
the pion mass. For Hak > Mr/2, the quark determinant becomes dependent on Hak and highly 
oscillatory on the statistically significant domain of real gauge fields. This means that zeros of 
the quark determinant around jig, ~ m,/2 move closer to that domain. It would be natural to 
conclude that the statistically dominant region of the original integration cycle is decomposed 
into multiple Lefschetz thimbles around Hak ~ m,/2. The complex phase of the fermion deter- 
minant becomes ill-defined at its zeros. Thus, the constant-phase condition of the integrand on 
a Lefschetz thimble would be hard to be satisfied without such decomposition, when zeros of 
the fermion determinant are located in the vicinity of the Lefschetz thimble. 

Because of the complexity of QCD, we cannot show these statements so far. Future study of 
QCD-like models based on Lefschetz thimbles will be crucial to develop our understanding on 
the baryon Silver Blaze problem and on the finite-density QCD. Since the statistically significant 
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domain of real gauge fields can be partially mimicked by chiral random matrix models and by 
instanton liquid models, it would be desired to study those models from this viewpoint of 
Lefschetz thimbles. 


4.4 Brief summary 


In this chapter, we discussed the sign problem appearing in the path-integral expression of the 
partition function from the viewpoint of the Lefschetz-thimble approach. 

In Sec. 4.1, we first establish a theorem which ensures the real-valuedness of physical observ- 
ables. Since we have to complexify the integration variables to deform the integration contour, 
it is not obvious whether physical observables take real values in a manifest manner. Our anal- 
ysis shows that the antilinear symmetry of the gradient flow is its key. If the system accepts 
the charge conjugation, we can construct the antilinear charge conjugation operation, and it 
ensures that physical observables take real values. 

This theorem is helpful to understand the sign problem appearing in the mean-field ap- 
proximation and in the perturbation theory. Since the fermion chemical potential makes the 
effective action complex, the gap equation in the mean-field approximation has no solution 
when the sign problem exists. This problem can be solved by applying the background field 
method and the Lefschetz-thimble analysis, which gives the mathematical foundation for the 
results of previous studies [113, 114]. This knowledge is used to analyze the heavy-dense QCD 
for a lucid demonstration of this method. 

In Sec. 4.2, we apply the Lefschetz-thimble method to the one-site Fermi Hubbard model in 
order to understand the essence of difficulty hidden in the sign problem of finite-density QCD. 
The one-site Hubbard model is quite simple but its path-integral expression has the severe 
sign problem. If the phase-quenched approximation is applied, the fictitious phase transition 
appears at u = —U/2 although the number density must stay zero until u = 0. This reminds 
us what happens for QCD, because phase-quenched QCD experiences the second-order phase 
transition at H4 = m,/2 although there is nothing until the liquid-vapor phase transition 
H4 = (mn — B)/3. Therefore, the one-site Hubbard model has a severe Silver Blaze problem. 

Using the Lefschetz-thimble approach to the path integral, we succeeded to unveil the secret 
of this Silver Blaze problem for the one-site Hubbard model. For u < —U/2, we can apply 
the complex mean-field approximation described in Sec. 4.1. Around u = —U/2, the Stokes 
jumping happens because of fermionic zero modes, and the topological structure of Lefschetz 
thimbles drastically changes, and the situation totally changes. In order to explain the Silver 
Blaze phenomenon, the interference of complex phases among complex saddle points play a 
pivotal role. We discuss Lee-Yang zeros within the semiclassical analysis, and the correct 
phase transition can be explained in a clear way. 

Based on the semiclassical analysis, the necessary number of Lefschetz thimbles is estimated. 
It is proportional to the inverse temperature, 3U/27, and we confirmed its correctness by using 
the numerical computation. We hope this analysis suggests the general strategy to evaluate 
the difficulty of tackling the sign problem using the Lefschetz-thimble approach. 

In Sec. 4.3, we discuss the implication on the Silver Blaze problem of finite-density QCD 
based on our analysis of the one-site Hubbard model. These theories share the difficulty of 
the sign problem from various points of view. Outside of the pion condensation of the phase- 
quenched QCD, the contributing Lefschetz thimbles would not have different complex phases. 
This means that we can apply the complex mean-field approximation or also other conventional 
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approaches to the sign problem. If u, > m,/2, then the fermion determinant would divide the 
integration cycle into multiple Lefschetz thimbles with different complex phases. In summing 
up those phases, the baryon number density stays zero, which breaks our naive expectation 
obtained from the pion condensation in the phase-quenched theory. This speculation must be 
justified/denied in the future to deepen our understanding on the baryon Silver Blaze problem. 


Chapter 5 


Complex Langevin method revisited 
via Lefschetz-thimble integrals 


The complex Langevin method has been considered as a sign problem solver for a long time. 
Even how severe the sign problem is, this method solve it for some cases. However, it is also 
known that the complex Langevin method shows wrong convergence in other cases, and this 
situation is annoying when we apply it to the strongly-correlated systems. 

In Sec. 5.1, we show a simple criterion for incorrectness of the complex Langevin method 
by relating it with the Lefschetz-thimble path integral through the Dyson—Schwinger equa- 
tion. This clearly shows that the naive complex Langevin method cannot solve the Silver Blaze 
phenomenon at finite-density QCD for u > m,/2. In Sec. 5.2, we propose to revise the com- 
plex Langevin method to evade the proof of its incorrectness. We check our criterion for the 
incorrectness numerically for two zero-dimensional examples, and test the modified complex 
Langevin method. 

The discussion in this chapter is adapted from Ref. [138]. 


5.1 Complex Langevin method and its failure 


After a brief review of the complex Langevin method, we show that it gives wrong convergence 
in generic cases by relating it with the Lefschetz-thimble path integral. 


5.1.1 Brief review on basics of complex Langevin method 


For simplicity, we discuss an oscillatory integral of one variable x, which can be extended to 
multiple integrals in a straightforward way, 


Owe Z [ dreS@"O(). (5.1.1) 


The classical action S(x) is complex valued in general, which makes the Monte Carlo simulation 
of (5.1.1) difficult because of the sign problem. One proposal for solving the sign problem is 
to take the noise average of the solution of the classical equation of motion with random 
noises [139-142]. For this purpose, we need to solve the Langevin equations for complex values 
z along the fictitious time direction 


Oozn(8) = —S'(z,(0)) + VANO), (5.1.2) 
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with 0 being the continuous fictitious time. 7(@) is real Gaussian noises satisfying (7(@)), = 0 
and (7(8)n(@")), = 26(6—6'). The function z, shows the classical trajectory of the system when 
the random force ņ is fixed. Since the action S(x) is complex, the right-hand side of Eq. (5.1.2) 
is complex. Thus, complexification of the variable z is unavoidable, and this method is named 
the complex Langevin method. 

For that purpose, let us review the result of Ito calculus. We can restrict our consideration 
on the set of holomorphic operators O(z) for later purpose. The stochastic calculus shows that 


d 
5 (OC-a(6)) n= ROM (29(8)))y — (O" (8) S"(2n(0))))n (5.1.3) 
If we assume that the expectation values reach equilibrium, then the set of expectation values 
satisfy the Dyson-Schwinger equation: 
dim (O'(2n(9)S"(2n(9))) mn = lim A(O” (2n(8)))n- (5.1.4) 
— +00 — o0 
The real Langevin method is equivalent to the path-integral quantization; if S(x) were real, the 
ensemble average (O(x,(@))),, can be shown to converge into (5.1.1) as 0 — oo by rewriting this 
Dyson-Schwinger equation into the Fokker—Planck equation [143]. There is a folklore theorem 
stating that the complex Langevin method also gives correct expectation values, but let us first 
show that it is true only for very restricted cases in the next subsection. After that, we propose 
a new prescription to evade this breakdown. 


5.1.2 Semiclassical inconsistency 


In the following, we assume that expectation values of holomorphic operators O(z) converge 
as 0 — oo. As we have already seen, It6 calculus shows that they must satisfy the Dyson- 
Schwinger equation 

(O(2n)S"(Zn))n = RO" (Zn) n- (5.1.5) 
Here the argument 6 = oo of z is omitted. The Lefschetz-thimble method shows that the 
expectation value of the holomorphic operator O(z) at 6 — oo can be represented, by using a 


sum of complex contour integrals on steepest descent paths J, of ReS(z) (a € £ labels the 
critical points z, of S) [64, 65, 76], as 


(O(&))y = D [ dze-S@/4Q(z). (5.1.6) 


o€y 


Here, d, is a complex number in general. If and only if dẹ is an intersection number (Ke, R) 
between the steepest ascent path X, and the original integration path R, the original integration 
(5.1.1) is obtained: 


Ce te i dreS@V"Q(2) & dy = (Ko, R). (5.1.7) 


We already proved in Sec. 2.2 that each integration on Lefschetz thimble respect the Dyson- 
Schwinger equation. Another linearly independent solution does not exist since the Dyson- 
Schwinger equation requires that S’(z) = 0 in the limit A — 0, which is nothing but the saddle- 
point condition. Thus, all the possible integration cycles must be associated with complex 
saddle points. 
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Let us take the semiclassical limit A + +0 of (5.1.6). For semiclassical analysis, we expand 
S(z) around each complex saddle point z, as 


S22) = 8,4 ia + O(623). (5.1.8) 


Let us first analyze the left hand side of (5.1.6). If Rew, > 0, the solution of the equation of 
motion (5.1.2) can converge into z, as 6 + oo in the limit of A — 0. On the other hand, it is 
impossible for Rew, < 0. In the semiclassical approximation, we have 


(Ole S e0lz,); (5.1.9) 


where c, > 0, and co = 0 if Rew, < 0. Next let us analyze the right hand side of (5.1.6). In 
the semiclassical approximation, the integration along the thimble becomes 


2rh 

| dz e S@O(z) ~ fe Ol20)- (5.1.10) 
1 /Qth _g 

a5 2p) a Pelt .1.11 

o= zN u, do (5.1.11) 


Notice that c, > 0, however, d, must be an integer (Ko, R) in order to give the original in- 
tegration (5.1.1). These two statements clearly contradict with one another for generic cases, 
especially when there are several dominant saddle points with different complex phases. There- 
fore, complex Langevin methods cannot reproduce original integrals at least in the semiclassical 
limit. 

Let us give a notice about subdominant saddle points. Since the Borel resummation of 
higher-order perturbations may cancel the subdominant contributions of complex saddle points, 
we need a careful case-by-case study in order to conclude inconsistency from the above argu- 
ments. Nevertheless, we strongly believe that the above argument is still robust and applicable 
to judge the breakdown of complex Langevin method for wide class of classical actions S(x). 
We give a clear and simple criterion for the incorrectness of the complex Langevin method. 

Let us give a few comment on previous studies. There is a formal proof [144, 145] on the 
correctness of complex Langevin method, but it relies on several nontrivial assumptions. Recent 
study [146] elucidates that this assumption is questionable especially if some singularities exist 
in complex classical action S(x), such as poles or branch points. Our semiclassical analysis show 
that this breakdown happens in more generic cases. Furthermore, our criterion for incorrectness 
of the complex Langevin method can be applied before performing the complex Langevin 
simulations. 


Now, we reach 


5.2 Modification of complex Langevin method 


In the previous section, we showed that the complex Langevin method gives an unphysical 
solution of the Dyson—Schwinger equation in generic cases. We give a trial to modify the 
complex Langevin method in order to evade the semiclassical inconsistency. Although this trial 
is still incomplete, it may open a way to perform numerical multi-thimble lattice simulation, 
which has not yet been achieved in the current technology. 
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5.2.1 Proposal to modify complex Langevin method 


Let us assume that if Rew, < 0 then (Ko, R) = 0 for a complex classical action S(z). Otherwise, 
the number of contributing saddle points in the semiclassical limit is different between the 
complex Langevin method and the original integral. In such cases, I have no idea how to 
modify the complex Langevin method to evade this inconsistency. Therefore, we require that 
S(x) satisfy this restriction. 

Let us propose a way to circumvent the inconsistency of semiclassical analysis by adopting 
this conjecture as a working principle. We denote the equilibrium distribution of the complex 
Langevin method by P i.e., (O(z)), = f dxdy P(x, y)O(« + iy). In the semiclassical limit, P 
will be approximated as a sum of localized distributions P, around complex saddle points z,, 
as P = )°.coP;, which gives the expectation value (5.1.9). Here, each P, is a normalized 
probability density that localizes around z,. In order to evade the inconsistency, we need to 
introduce another distribution function, 


P(x,y) = Io R)4/ he SIP ay). (5.2.1) 


Here, Z is a normalization factor, as usual. As a result, the new average f dxdyP (x, y)O(x+iy) 
follows the correct semiclassical limit, and thus the Dyson—Schwinger equation is not violated at 
the semiclassical level. In order to represent this new average as an ensemble average, introduce 
the following phase function, 


(Ko, R)4/ te Se/h 
(x,y) =X Me 


Co 


Xo(2,Y), (5.2.2) 


oO 


where x, is almost the characteristic function of supp(P,) so that y,P, ~ 6,;P, because of 

locality of P,. When c, = 0, we must put x, = 0. Using this phase function, we can easily 

find that 

(®(2, Yn )O(Sn + Yn) 0 
(®(2y,Yn))n l 


This seems to be the only way to circumvent the inconsistency in the semiclassical limit for 
the complex Langevin method. However, we must give a notice that this prescription contains 
several unsolved problems to be tackled: (i) Phase factors are calculated in the semiclassical 
limit, and they can suffer from higher order corrections of A. (ii) In order to apply this pre- 
scription, we must compute c, in (5.1.9) and its higher order corrections in A. (iii) When h is 
not so small, P does not necessarily localize around saddle points and thus we may not be able 
to find good characteristic functions Xo. We must evaluate an error affected by this failure of 
the prescription in order to give a reliable prediction. 

The first problem (i) is not a big deal, because we can compute higher order corrections by 
usual perturbation theory. The third problem (iii) can be a fatal disease for strongly coupled 
systems if quantum fluctuations give huge corrections. This must be an important future task 
when the complex Langevin method is applied. The second one (ii) may look similar to the 
first one, but there is no standard way to compute c, to my best knowledge. Since the Fokker- 
Planck operator of the complex Langevin method is non-Hermitian, the developed techniques 
in quantum mechanics are not necessarily useful. 


J axayPe. noe + iy) = (5.2.3) 
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In order to demonstrate that this prescription may work well, we put a working hypothesis 
on cz. Let us assume that 
f dze 5) 


This solution means that the complex Langevin method gives an extension of the so-called 
phase quenched approximation so as to include complex saddle points: 


(O(2n))n & Do R), [eNO 20). (5.2.5) 


o 


Co = (Ko, R} 


— 
A 

2 
A 

~~ 


(5.2.4) 


Note that there are many other possibility to satisfy semiclassical consistency. This is just a 
one of possibilities and it has not been proven that this is indeed the solution of the complex 
Langevin method. Computation of c, directly from the complex Langevin method is of great 
importance if one would like to solve the sign problem using it. 

Once we adopt this working hypothesis on c,, we can write down the phase function ® in 


an explicit manner as 
[wo] —itm(s 
p= y (Ky, R)4{ Te T 5.2.6 
DT (5.2.6) 


The expectation values of holomorphic operators are computed by (5.2.3). 


5.2.2 Numerical test of the proposal 


In this section, we check that our criterion for incorrectness of the complex Langevin method, 
and test our proposal of modification by applying it to two examples with and without a 
singular drift term. As a nonsingular potential, we consider the double-well potential with a 
complex linear term, which is called the complex double-well potential here. As a singular one, 
we consider the one-site Fermi Hubbard model which is studied in Sec. 4.2 using Lefschetz 
thimbles. 


Complex double-well potential The first example is a nonsingular action 


rt r? 


S(x) = i a iaz, (5:27) 


with a > 0. This action has three saddle points on the complex plane, and let us see their 
properties in the limit œ — +00. In the case a > 1, the x? term in (5.2.7) becomes negligible 
and we can give simple expressions for three saddle points as 


1/3 57i/6 


z walde/6 2 ~a en is le (5.2.8) 


Around these points, we evaluate the classical action as 


3 ; , 
Bae oe). & gare + 3a2/3er/3 557, (5.2.9) 


K 


3 , : 
S (z2 + 62) = A + 3a? ez, (5.2.10) 


3 
S(z+6z) > = oe — 8a7/3§2?. (5.2.11) 
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Figure 5.1: Test of the complex Langevin method for the complex double-well potential S(x) = 
x*/4—x?/2+iax with a > 0. (a) Scatter plot obtained by the complex Langevin for a = 3.8. 
For large a, two local distributions around saddle points are well separated. (b) Expectation 
values of z?. The red boxes are obtained by the original complex Langevin method, and the 
blue circles by our new formula (5.2.3). (Figure (b) is taken from Ref. [138].) 


Only two of them, zı and 22, have positive Re|w,] as we can easily see, and contribute to the 
Langevin simulation in the semiclassical limit a — oo. In the present model, the same is true 
for the Lefschetz-thimble method, and only these two saddle points give contributions in the 
Lefschetz-thimble decomposition: 


rte) = ff does f tow 
Ji I2 
3v3 
T ad = E). (5.2.12) 


K 


32/3 


The last expression is the asymptotic expression which is valid only for a >> 1, but it clearly 
tells us the oscillatory behavior of the partition function. This oscillatory behavior comes from 
interference of complex phases between dominant saddle points 21,2, and thus the naive complex 
Langevin method is expected to give wrong results for this model. 

We numerically solved the complex Langevin equations with the fictitious time step € = 
5.0 x 10-7. We adopted a higher order algorithm [147], and errors were estimated by using the 
jackknife method. The result for (z*) is shown in Fig. 5.1 (b). The red squares show the result 
for the naive complex Langevin method, and they are clearly incorrect as we expected. 

Let us try our proposal on modification of the complex Langevin method. Since there are 
two contributing saddle points with Refz1] > 0 and Ref|z2] < 0, we divide the complex plane C 
into two regions for positive and negative real parts (See Fig. 5.1 (a)). We introduce the phase 
function by 


n 
=y |S”( some e ilm S(zc) 6((—1)°** 2), (5.2.13) 


i 
o=1,2 P (2o) 


where O(x) is the Heaviside step function. Here, saddle points 21,2 are computed exactly, 
although we only showed their approximate expressions for a > 1 in (5.2.8). The result 
for (®(x,y)z*)/(®) is shown with blue circles in Fig. 5.1 (b), and those points are roughly 


5.2. Modification of complex Langevin method 85 


consistent with the exact result (black dashed line). Especially for large a’s, our prescription 
on the modification works nicely for this model. 

We numerically checked whether the Dyson—Schwinger equation is satisfied for this model, 
but it turns out not to be satisfied for z” with n > 3. The equilibrium probability distribution 
seems to show a power-law decay in the Im(z) direction for this model, and this might be a 
reason for this illness. Our criterion for incorrectness still gives the sufficient condition for the 
complex Langevin method to break down, and we can clearly see it in this example. What is 
surprising is that our prescription of the modified complex Langevin works nicely for z” with 
n <4 in this case, and it is an open problem to understand its reason. 


One-site Fermi Hubbard model Next, we analyze the action with a logarithmic singular- 
ity. We analyze the one-site Fermi Hubbard model studied in Sec. 4.2, which is the simplest 
model and shows the severe sign problem like finite-density QCD (see also Refs. [128, 132—134]). 
For reader’s convenience, let us briefly describe properties of this model again. 

The Hamiltonian of this model (4.2.2) is given by 


A = Ufyiy — u(y + îy), (5.2.14) 


with ñ, = ĉłtê,. At the zero temperature, this system shows the first-order phase transitions at 
u/U =0, 1. A path-integral expression of the partition function Z = tr [exp (-så ) is given 
by the following one-variable integral (4.2.19), 


Z= yzy |% 1 +exp£ | ip» eg Er) (5.2.15) 
2nU i E 2 2U P8 


The expectation value of the number density can be computed as n = Im(y)/U. For u/U < 
—0.5 or u/U > 1.5, behaviors of the number density can be easily described just by solving the 
classical equation of motion, which comes from the fact that the u-dependence of the fermion 
determinant becomes trivial. What was difficult to understand in this path-integral expression 
is the step-function behaviors of n for —0.5 < u/U < 1.5, and the sign problem in this region 
is expected to have strong connection with that of finite-density QCD at m,/2 < u < my/3. 

Let us again stick to the case —0.5 < u/U < 1.5, since this is the best testing ground for the 
complex Langevin method. In this region, the contributing complex classical solutions (4.2.28) 
are labeled by integers m € Z as (BU > 1) 


2 


2 


U 307 — 
oni (u B) +T (2mm it Ty 3 (5.2.16) 
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At these saddle points, the classical action is evaluated as 


Qn? 


Re (Sin = So) = ao” (5.2.17) 
H 1 

Im Si, ~ 2 PS, 21 

mS. nm (4 + 3 (5.2.18) 


In this model, the most dominant contributions comes only from a saddle point zo, and then, 
strictly speaking, one cannot apply our proof of incorrectness of the complex Langevin method. 
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Figure 5.2: Test of the complex Langevin method for the one-site Fermi Hubbard model. 
(a) Scatter plot obtained by the complex Langevin method for BU = 30 and u/U = 0. This 
has a remarkable similarity with Fig. 4.3 (b). (b) The red boxes are obtained by the original 
complex Langevin method, and the blue circles by our new formula (5.2.3). (Figure (b) is taken 
from Ref. [138].) 


However, other saddle points zm (m #4 0) also give almost dominant contributions at GU — oo 
as long as |m| is not too large. This suggests that our criterion can still conclude that the 
complex Langevin method shows wrong convergence in this model. That is, because of the 
difference of complex phases among Sm, the complex Langevin method is expected to give a 
wrong answer for this model. 

Figure 5.2 (a) shows the scatter plot obtained by the complex Langevin method obtained by 
the complex Langevin method. Here, Im(z) is subtracted by the offset u +U/2 in order to look 
at details of the distribution. They does not show step-function behaviors and instead show the 
linear dependence on p. It is quite similar to the structure of Lefschetz thimbles in Fig. 4.3 (b), 
and the distribution is nicely separated by fermionic zero modes. The result for the number 
density of the naive complex Langevin method is shown with red squares in Fig. 5.2 (b). This is 
expected because the position of the complex saddle points is Im(z,,)/U œ u/U +1/2, and this 
indicates that the complex Langevin method at least loses the complex phases among saddle 
points. 

We introduce the phase function ® by 


y= X eim Smg ( (2m —1)r < Ba < (2m + yr). (5.2.19) 
The result of the modification, Im(® - z)/(®), is shown with blue circles in Fig. 5.2 (b). Al- 
though we cannot observe the perfect agreement between our modified computation of the 
complex Langevin method and the exact result, step-function behaviors can be observed after 
modification. It is an important future study to refine the modification procedure and explore 
the Silver Blaze phenomenon using the modified complex Langevin method. 


5.3 Brief summary 


In this chapter, we discussed the relation between the Lefschetz-thimble path integral and 
the complex Langevin method through the Dyson—Schwinger equation. Establishment of this 
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relation enables us to give feedbacks from our new findings on the Lefschetz-thimble approach 
for other approaches to the sign problem, and vice versa. 

In Sec. 5.1, we formally relate the complex Langevin method and the Lefschetz-thimble 
path integral by using Ito calculus. Since It6 calculus gives the Dyson—Schwinger equation if 
the Langevin dynamics equilibriates, the ensemble averages of the complex Langevin method 
must be given by a complex linear combination of Lefschetz-thimble integrals. Only when it is 
given by the linear combination with appropriate integers, does the complex Langevin method 
simulates the original quantum system. 

Using the semiclassical analysis, we derived the sufficient condition for the complex Langevin 
method to show the wrong convergence. This condition is quite simple: The complex Langevin 
method is always wrong if the classical action has several dominantly contributing saddle points 
with different complex phases. 

Since we have identified a missing part in the complex Langevin simulation, we can propose 
a way to change the complex Langevin method in Sec. 5.2. Since the phase factors appear- 
ing in summing up multiple Lefschetz thimbles are missed, we introduce them and perform 
the reweighting. This is the minimal modification in order to circumvent the semiclassical 
inconsistency in Sec. 5.1. 

Our proposal is checked numerically. If we apply the complex Langevin method to the 
one-site Hubbard model, it shows the fictitious phase transition which we observed with the 
phase-quenched approximation in Sec. 4.2. After performing the reweighting, the result on 
the number density is improved and shows a rapid transition at correct values of the chemical 
potential. 
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Chapter 6 


Summary and Outlook 


In Chapter 2, I first give a brief review on a foundation of the Picard—Lefschetz theory. This 
gives an extension of steepest descent methods by applying it to the multiple oscillatory integral, 
and the steepest descent integration cycles are called Lefschetz thimbles. In order to understand 
how it works, I compute the real-time Feynman kernels of free particles and the harmonic 
oscillator. 

In Chapter 3, spontaneous breaking of chiral symmetry is discussed from the viewpoint of the 
Lefschetz-thimble approach. This study tackles one of the important aspects of nonperturbative 
QCD using simple 0-dimensional fermionic models, Gross-Neveu and Nambu-—Jona-Lasinio 
models. In both cases, an auxiliary bosonic field for the chiral condensate is introduced, and its 
integration cycle is decomposed into Lefschetz thimbles. This allows ones to study the analytic 
property of the partition function for complex four-fermion interactions before performing the 
path integral. Thanks to this property, one can easily compute Lee-Yang zeros of these models 
in the large-N limit. 

In this chapter, I also develop the computational tool of Lefschetz thimbles with the contin- 
uous chiral symmetry. The Lefschetz thimble method can be naively applicable only when all 
the saddle points are non-degenerate, which is not the case if there exists an exact continuous 
symmetry. Even when the continuous symmetry is explicitly broken, the remnant of the sym- 
metry introduces a practical difficulty in its application because of the appearance of quasi-zero 
modes. I find that the quasi-zero mode of the gradient flow can be decoupled approximately, 
and Lefschetz thimbles become computable in an efficient way. This technique turns out to be 
quite useful, and it is applied for the 0-dimensional Nambu—Jona-Lasinio model. 

In Chapter 4, I consider two examples with the sign problem in order to understand the 
origin and structures of the sign problem. The first one is the Polyakov-loop effective model 
of the heavy-dense quantum chromodynamics (QCD). In order to study it, I developed the 
mean-field approximation based on the Lefschetz-thimble approach so that it is applicable even 
if the sign problem exists. One can judge whether this works nicely by computing topologies 
of Lefschetz thimbles. This gives the mathematical foundation to study the Polyakov-loop 
effective models at finite quark densities. 

The second one is the one-site Fermi Hubbard model, which is easily and exactly solvable but 
has the severe sign problem in the path-integral expression. In this case, the above mean-field 
approximation is not applicable, and I elucidate its reason by studying topological structures 
of Lefschetz thimbles. I show that interference of complex phases among complex classical 
solutions play a pivotal role to understand the sign problem and the correct phase structures 
using path integrals. This interference is shown to become also important in the finite-density 
QCD with light flavors if the baryon chemical potential exceeds the pion mass. 
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In Chapter 5, I relate this newly developing technique, the Lefschetz-thimble path integral, 
with a conventional approach to the sign problem, the complex Langevin method. I show that 
the original complex Langevin method cannot give correct results if interference of complex 
phases among classical solutions is important. This gives the simple criterion for incorrectness 
of the complex Langevin method without doing its numerical simulation, and suggests that the 
complex Langevin method is not applicable to the finite-density QCD. To resolve this failure, 
I propose to modify the complex Langevin method and it is numerically tested. Our proposal 
on modification improves the result for the one-site Hubbard model but it is not perfect yet. 
This study strongly suggests that the Lefschetz-thimble approach is quite powerful to invent 
or justify other sign problem solvers. 

In this thesis, I elucidate what is an essential difficulty of the sign problem when the path- 
integral formulation is applied. This requires some changes or modifications of conventional 
approaches to the sign problem, and this study may be able to infer how to do them. For 
that purpose, I must study further details and evaluate the amount of computational costs. 
Especially, it is important to understand the structures of Lefschetz thimbles when the system 
has infinitely many degrees of freedom. 

In order to explore the phase diagram of QCD using the path-integral approach, I believe 
that we must further develop the theoretical and technical foundations of Lefschetz-thimble 
methods. It is quite important to fully understand how the sign problem in finite-density QCD 
is correlated to the low-energy dynamics of QCD characterized by confinement and spontaneous 
chiral symmetry breaking. Each aspect of the problem is partly explored in this thesis, and 
I would like to know how these results play into the realistic problems. As a next step, it is 
interesting to revisit the structure of the sign problem in the chiral random matrix model and 
the instanton liquid model. These models simplify the dynamics of gauge fields in QCD, and 
thus one would be able to learn how the chiral symmetry breaking affects the difficulty of the 
sign problem in an explicit way. It is fascinating if these studies shed a new light on the cold 
and dense nuclear matters and a longstanding and new problems of neutron stars. 

So far, the path integral on Lefschetz thimbles relies on sophisticated mathematics. Its 
applications can be discussed based on solid mathematical foundations thanks to this aspect. 
However, if one could relate this analysis with our conventional physical intuition on quantum 
field theories, then this method must have wider applications and become practical tools. I 
believe that this is possible, because this method is suitable for the semiclassical study of 
quantum physics. Dynamical aspects of field theories, such as real-time phenomena, have vast 
unexplored areas compared with static ones, and I would like to tackle those problems in future 
by extending the applicability of this technique. 


Appendix A 


Homology 


A.1 General property 


Homology can be defined in a quite abstract way, and it is sometimes useful to compute the 
homology since we can deal them without considering details of its construction. Therefore, let 
us just remark the Eilenberg-Steenrod axioms [148] of the homology at first. 


Kilenberg—Steenrod axiom 

For every pair (X, A) of topological spaces, i.e. A C X, Abelian groups H;(X,A) (i > 0) 
and boundary maps 0: H;(X, A) > H;_1(A,)(=: H;_1(A)) are assigned, and a map of pairs 
f : (X, A) > (Y, B) induces a group homomorphism f, : H,(X,A) > H,(Y, B). They further 
satisfy the following five axioms: 


1. [Naturalness] The boundary map 0 is natural, i.e. for a map f : (X, A) > (Y, B) the 
induced maps must commute with 0; f, 00 =00(fla).«. 


2. [Homotopy] If two maps f,g : (X, A) > (Y, B) is homotopic, the induced maps fs, gx 
to the homology are the same. 


3. [Excision] Let (X, A) be a pair and U C A. If the closure of U is inside the interior of 
A, then the inclusion (X \ U, A \ U) > (X, A) induces H;(X \ U, A \U) = H;(X, A). 


4. [Dimension] For the one-point set x, H;(«) = 0 for i > 0. In the singular homology, 

5. [Additivity] Let X is given by a disjoint union of X,’s, then H;(X) = @, Hi(Xa)- 

6. [Exactness] For each pair (X, A), two inclusions i: A > X and j : (X,0) => (X, A) 
induces the long exact sequence 


E (Ay Hx) Ss Be Ay BA) See, (A.1.1) 


With these axioms, the homology can be uniquely determined. For some topological spaces 
like spheres, one can compute their homologies directly from these axioms. 

If one were familiar with the homology group H,,(X) but not with the relative homology 
group H,,(X, A), one can think of H,,(X,A) as H (X/A) for n > 0 (for singular homologies). 
For n = 0, Ho(X, A) = Ho(X/A, *) and it is the free module of one rank less than Ho(X/A). 
These properties follow from the excision axiom. 
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Mayer-—Vietoris exact sequence 


The Mayer-—Vietoris sequence provides a helpful way to compute the homologies, and this 
theorem follows from the above Eilenberg—Steenrod axioms. Let A,B C X be topological 
subspaces of X, and X be covered by interiors of A and B. We can consider natural inclusions 
i: AAB A, j: ANBOB,k: AG X, andl: B— X, then the Mayer—Vietories sequence 
states that the following is the long exact sequence: 

= Hi Ane) BAe.) SS aos ene, Ca 
The map ô, needs a little explanation. Let us set U = B \ A, then the closure of U is contained 
in B because X is covered by interiors of A and B by the assumption. Using the excision 
property, one finds 


A, (X, B) ~ H,(X \ U,B\ U) = H(A, ANB). (A.1.3) 
Then, ô, is given by the composition, 
ð; : Hy(X) > H,(X,B) ~ H,(A, AN B) 2 Hy1(ANB). (A.1.4) 


For spheres 9*, this states that H,,(9") = H,-1(S*~') by dividing it into two hemispheres. 
In the case of singular homologies, one can iteratively show that Ho(S*) = H,(S*) = Z and 
HS") = 0 for n Æ 0,k: 


A.2 Smale—Witten theory 


Morse homology provides a useful way to understand homologies of compact smooth mani- 
folds. Semiclassical analysis of a supersymmetric quantum mechanics gives a physical idea of 
homologies [149, 150]. For its relevant extension in our context, see Refs. [151-154]. 


Morse inequalities Let (M,g) be a smooth Riemannian manifold with finite dimension, 
and consider h : M — R. For a moment, let us assume that M is compact. The function h 
is called a Morse function on M if all the critical points p € & = {x € M|dh(x) = 0} is 


non-degenerate, i.e. 
di O*h 
© | arið 


For each critical point p, the Morse index A(p) is defined by the number of negative eigenvalues 
of the Hessian. 
The strong Morse inequality gives 


NOPO — X (dim, (M) = (1+ t) X Ont” (A.2.2) 


ped n>0 n>0 


#0. (A.2.1) 


P 


for some nonnegative integers Q,. This gives #{A(p) = n} > dim H, (M). This inequality 
can be shown for any homologies satisfying the Eilenberg-Steenrod axioms. Especially if all 
Qn = 0, h is called a perfect Morse function. 
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Gradient flow and stable/unstable manifolds Let us still stick to the smooth compact 
Riemannian manifold M. Taking a local coordinate x‘, Morse’s downward flow equation is 
defined as 

dat ay dh 

dt Ors 
There are two types of solutions in general; one is the constant solution at a critical point, 
x(t) = p, and another one is a one-dimensional curve connecting two critical points. Let’s 
call this one-dimensional manifold as an “instanton”, because this is indeed an instanton of a 
supersymmetric quantum mechanics [149]. For each critical point p, we define the descending 
and ascending manifolds by 


(A.2.3) 


lim z(t) = p} ; (A.2.4) 


i= {0 E€ m | „lim z(t) = p}, Kp := {0 eM jim 
respectively. The stable/unstable manifold theorem states that J, and K, are smooth mani- 
folds, which are homeomorphic to open disks, and dim J, = A(p) and dim kK, = dim(M) — X(p). 

The pair (h, g) of the Morse function and the metric is called Morse-Smale if J, and K, 
are transverse to each other for any p,q € X: For any r € Jp N Ka, Trp OT Ky = T,M. This 
condition gives an exact sequence, 


0> T: hN K > TI) 8 T-a) > T.M = 0, (A.2.5) 


for r € Jp N Ka. This is a generic condition, in other words, small perturbations can always 
satisfy this even when a Morse function h is not Morse-Smale. We assume that (h, g) is Morse- 
Smale in the following. 

Since dim(.J, N Ka) = A(p) — A(q) and 


Jp 1Kq = {x(0) | z(—00) = p, x(+00) = q}, (A.2.6) 


the dimension of the moduli space of instantons is given by A(p) — A(q). Therefore, if A(q) = 
A(p) — 1, the number of instantons connecting two critical points p and q is finite. The number 
of such instantons is precisely necessary information for computing homologies of M [149]. 
Consider the free Z-module generated by descending manifolds J, graded by Morse indices, 


C(M,h) = @®Ca(M, h), Cr(M,h)= X ZY, (A.2.7) 


ped, A(p)=n 


The boundary operator 0: Cn > Cn-1 is given by 


= >. needs (A.2.8) 


ged, A(q)=A(p)—1 


where n(p,q) is a number of instantons connecting p and q. Here, an instanton is counted as 
+1 whether it preserves the orientation of the exact sequence (A.2.5) or not. One can check 
that 0? = 0, and the homology of the complex C(M, h) is nothing but H,(M). This also shows 
the strong Morse inequality. 
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Sketch of Witten’s proof Let us recall the de Rham complex at first. Denoting the set of 
n-forms on M by 2"(M), the de Rham complex is 


OF PoS aM) S ee S 4m). (A.2.9) 
Cohomology of this complex is called the de Rham cohomology Hjp(M), and the Poincaré 
duality states that H;,(M) ~ H. a ae ) as R-vector spaces. One can introduce an L? inner 
product on Q” by 


(a|B) = f ane, (A.2.10) 


where x is the Hodge star operator and @ is the complex conjugation of a. If degrees of a and 
b are different, we set (a|3) = 0. Let us regard Q* as a “Hilbert space” of quantum mechanics. 
As a “Hamiltonian”, we use the positive-definite operator 


2 
Ho = ddt + dtd (- -> zar) , (A.2.11) 


where (d'a, 8) = (a, d8) (The last equality is true only on a flat local coordinate). This is 
nothing but the Laplace operator. The Hodge theory tells us that H7,(/) is nothing but the 
space of vacua of degree n, i.e., Hip(M) = {a € Q” | Hoa = 0}. 

Instead of the exterior derivative d, one can consider another operator d; = e~“"de"” for 
t > 0. Here h is a Morse function, and assume that (h, g) is Morse-Smale. Since this is just a 
conjugation, the de Rham cohomology is invariant for any t. Let us consider a Hamiltonian 


H, = dd! + did, (A.2.12) 


then we would like to compute the space of vacua in order to obtain (co)homologies. In 
Ref. [149], the spectrum of H; is computed in the limit t > oo semiclassically. 

Because of the anticommuting property of differential forms dz’ A da? = —dz/ A dz’, the 
space of n-forms Q” can be regarded as the Hilbert space including n fermions, and a> dz’? Aa 
gives the creation of ith fermion ®t. The Hamiltonian H; is aad as 


H =X (-® +? (Anh) +t a = Gace (A.2.13) 


n 


In the limit t + oo, low-energy wave functions must localize around critical points p € X 
because of the potential term (0,h)?. However, the bosonic part of Hamiltonians still have 
large energy due to zero-point oscillation. In order to cancel it, we need to make the fermionic 
energy maximally negative, and thus fermions must be occupied for negative eigenvalues of the 
Hessian of h and other fermions unoccupied. Let us denote this wave function by |p). 

So far, by approximating the system as harmonic oscillators around critical points, we 
obtained the space of classical vacua X"(M,h) with n fermions, and dimX” is the number of 
critical points p with Morse index X(p) = n. Because of instantons (A.2.3), some elements of 
X” may not be quantum vacua, which must be eliminated. The derivative ô : X” = X™*1, 


6|p) = > n(r,p)|r), (A.2.14) 


red, A(r)=A(p)+1 


satisfies 6? = 0. The cohomology of the complex (X*, ô) gives the de Rham cohomology [149]. 


Appendix B 


Path integral on Lefschetz thimbles 
and supersymmetry 


Lefschetz thimbles can be interpreted as ground states of M = 2 supersymmetric quantum 
systems [152-154] (see also Ref. [65]). It is possible to use this reformulation of the Lef- 
schetz thimble in order to compute expectation values of observables. This work is taken from 
Ref. [155]. 


B.1 Lefschetz thimble and SUSY quantum mechanics. 


Let us consider an N-dimensional real integral as a “quantum field theory” defined by a (com- 
plex) classical action S(x). In this theory our goal is to compute an expectation value of an 
“observable” O(a) defined by 


(0) = nf dg & O(a) , (B.1.1) 


where z = (a), a), ...,2)) € RY and the normalization M is chosen such that (1) = 1. 
The starting point in our discussion is to reformulate this theory in an equivalent and more 
treatable way using a complexified representation: 


(0) = fev daz P(z,Z) O(z) . (B.1.2) 
Here 2 = ol + ia? and Z = ol — ia? with 2 gO € R and f dzdzZ represents the 
integration over the whole complex plane; i.e. ee dx, lien dx. The choice of the gener- 
alized weight function P(z,Z) may not be unique. Indeed, a trivial example is P(z,Z) = 
NeS) [ó (z® — 2). At the cost of complexifying the variables, nevertheless, it is often the 
case that P(z,Z) could be endowed with more desirable properties for analytical and numerical 
computation than the original e®). 

A clear criterion to simplify the integral is to find P(z,Z) such that the phase oscillation 
can be as much suppressed along integration paths as possible, while in the complex Langevin 
method P(z,Z) is optimized to become a real probability. To suppress the phase oscillation, let 
us pick up a saddle point z, satisfying S’(z,) = 0. The steepest descent cycle or the Lefschetz 
thimble J- of the saddle point z, is defined with a fictitious time t as 
dx(t)  AReS 


pT par’ mle ti) (B13) 
j 


J= f-0) = zı (0) + ix(0) 
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This is a multi-dimensional generalization of the steepest descent path in complex analysis, 
which we will refer to as the downward path. The original integration path on the real axis 
in Eq. (B.1.1) can be deformed as a sum of contributions on J, weighted with an integer ms; 
ie, had = X, Mo Iz dz. In mathematics it is established how to determine m, from 
the intersection pattern between the steepest ascent (upward) path from z, and the original 
integration path [62, 64, 65]. It is important to note that ImS is a constant on each Lefschetz 
thimble for the application to the sign problem [67, 70]. 

In the following, let us restrict ourselves to N = 1 for simplicity, because the generalization 
is straightforward. So far, the Lefschetz thimble is constructed as a line, and let us find a two- 
dimensional smooth distribution P(z,Z) according to Ref. [65]. For that purpose, we define the 
“delta-functional one-form” (J) supported on the Lefschetz thimble so that 


f O(z)eFdz = f lT) A Olze Pdz. (B.1.4) 
A c 


For instance, ô(R) = 6(y)dy. Such delta-functional forms (Js) (on a Kähler manifold) have 
a path-integral expression from the supersymmetric quantum mechanics [152-154] (see also 
Secs. 2.8 and 4 of Ref. [65] for more details in this context). Integration (B.1.4) can be repre- 
sented as 


0 dx; Res 
(O) = N J Diep, mY] apli f an A e) 


co 


0 2 
x exp |- fa m( Sy $ a Je O(2(0)) eS) (y +i) (0) . (B.1.5) 
Here x,p are bosonic fields and z,~ are fermonic ghost fields, and z(t) > z, as t => —oo. 
We should note that an integration in terms of z is promoted to the path integral on z(t) for 
t < 0, while the observable and the weight O(z(0)) exp S(z(0)) are functions of z(0) only. Let 
us outline how these two expressions (B.1.4) and (B.1.5) are equivalent [65, 152-154]. We first 
integrate out p(t) to get the Dirac delta function, 


_f dz; Res dz; Res 
fo epli f arp ( E) =(4 + Aa, ) ‘ (B.1.6) 


This delta function constrains the path integral on x(t) to a gradient-flow line defining Lefschetz 
thimbles. Since z(—o0) > z,, this path integral for t < 0 gives a delta-functional support on 
Jo. However, the delta function produces an unwanted determinant factor. As is well-known, 
the path integral on ghost fields z(t), w(t) for t < 0 can eliminate that factor as 


d Res d Res 
for Dy exp - f dt mi( 585 + aaa) | = pet( $3 j E) . (B.1.7) 


oO 


Now, we obtain an integration over surface variables x(0), (0), and denote them by z, w. 
Locally, the Lefschetz thimble % can be expressed as zeros of a certain function f, then we 
can find that the path integral (B.1.5) eventually gives 


feryat Eunoe O (itis) = f DAEA (B18) 


which is nothing but the local expression of the original integration (B.1.4). Going back to 
(B.1.5), this shows that the so-called residual sign problem comes from the fermionic surface 
term Y1(0) + iv2(0) because one can identify y;(0) = dz; as above. 
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B.2 Hamilton dynamics 


Importantly, with these added fields, p;, 7, Yi, the action is BRST exact under a transformation; 
6x; = Wi, du; = = 0, ÔT; = —ip;, bp; = = 0 [65]. By definition the nilpotency 62 = 0 is obvious. 
Thanks to the Boundary fermionic operator in (B.1.5), the surface term is BRST-closed so long 
as the observables are holomorphic. This makes a sharp contrast to the complex Langevin 
method that could also acquire the BRST symmetry but it is violated by the surface term. 
Because of the BRST symmetry we can add any BRST exact terms without changing the 
original integral, and it is useful to insert = af dt p?. In summary, the effective Lagrangian that 
describes the fictitious time evolution is ven by the following topological theory: 


ee ee da; _ OReS i (45 , Res wb, 
Gg 9 Pi Pi dt oxi dt iJ OLO a 


-ô ; or - Oy ORES 
E MIDP T ag Ox; 


which is nothing but a Legendre transform of an effective Hamiltonian: 


= Ging Í Res, „ OReS 7 1?ReS die 


(B.2.1) 


with [x;,p;] = id,; and {ai bj} = 6;;. The fermion number F = fri + Froth is a conserved 
quantity of this Hamiltonian. After the time evolution from t = —oo only the ground state with 
the lowest energy eigenvalue remains, so that the generalized weight is given by P(z,Z)dzdz = 
W(z,z) A edz, where Y(z,Z) is the ground state wave-function and converges to 6(J,) in 
the limit e; > +0. Note that the weight factor exp S(z) is necessary in this formula, since the 
wave function designates only the integration cycle Js. We can further simplify this Hamilton 
problem by choosing ¢ = £1 = €2. Performing the conjugate transformation UV = e~R°5/*W’, the 
first derivative terms are eliminated as 


T E 9 1 /ðReS 2 B 1@ReS N 
“= 3 z i 2e ( Ox; 2 2 0x,0z; [ii v] : (B.2.3) 


This describes supersymmetric quantum mechanics with the superpotential ReS [65]. 
By restricting ourselves to the F = 1 sector, we can define the effective potential in a form 
of 2 x 2 matrix-valued function that amounts to 


ve — 1 [(2Res P OReS\*) / @ReS/dxz2  8Res/ðx1ð1 (B.2.4) 
of De | Oxy Oxy OReS/Ax,0r2 —OReS/Ox? j’ = 


and then the Hamiltonian is 


ef o o? 
i | + Vog. B.2. 
Hes 2 & Z) Ver (B.2.5) 


This reformulation provides one way to compute Lefschetz thimbles for numerical purpose by 
solving the ground state of the above H/, at some e [155]. 
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